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1 SUMMARY 


An innovative computational technique, Nonconical Relaxation (NCOREL), 
has been established for the treatment of three dimensional supersonic 
flows. The method is nonlinear in that it solves the nonconservative finite 
difference analog of the full potential equation and can predict the formation 
of supercritical cross flow regions, embedded and bow shocks. 

The method implicitly computes a conical flow at the apex (R = 0) of a 
spherical coordinate system and uses a fully implicit marching technique to 
obtain three dimensional cross flow solutions. This implies that the radial 
Mach number must remain supersonic. The cross flow solutions are obtained by 
using type dependent transonic relaxation techniques with the type dependency 
linked to the character of the cross flow velocity (i.e. , subsonic/ 
supersonic). 

The spherical coordinate system and marching on spherical surfaces is 
ideally suited to the computation of wing flows at low supersonic Mach numbers 
due to the elimination of the subsonic axial Mach number problems that exist 
in other marching codes that utilize Cartesian transverse marching planes. 

NCOREL has been applied successfully to a wide variety of body and wing 
configurations in the Mach number range from 1.5 to 3. The method represents 
a significant advancement in accuracy over linear methods. This report 
describes the practical application and use of the computer program NCOREL. 
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2 INTRODUCTION 


A new computational technique. Nonconi cal Relaxation (NCOREL), utilizing 
the full potential flow equation has been established for the prediction of 
inviscid nonlinear supersonic aerodynamics. This method has proven to be 
remarkably accurate for predicting flows in the low to moderate supersonic 
Mach number regime (i.e., 1.5 to 2.5). Unlike linear potential methods (e.g., 
panel methods), NCOREL can predict the formation of supercritical cross flow 
regions, bow and embedded shocks. 

Nonlinear effects are important to the aerodynamicist in the supersonic 
regime because he cannot intelligently or reliably tailor a wing to optimize 
lift and drag without knowledge of the formation of “embedded crossflow 
shocks" which can significantly alter the wing characteristics. Embedded 
crossflow shocks not only induce drag but can also lead to boundary layer 
separation. The high lift characteristics of supercritical crossflow regions 
can be beneficial to the aerodynamicist if he can exercise control over the 
formation and strength of crossflow shocks which are typically synonymous with 
these supercritical crossflow regions. Strong embedded cross flow shocks are 
most prevalent on wings that have subsonic leading edges. The strength of the 
crossflow shock is proportional to the angle of attack or lift coefficient, 
cross sectional thickness, and leading edge radius of curvature for an 
uncambered wing. Nonlinear effects are also important on wings that have 
supersonic leading edges or bodies whose slope is greater than the Mach 
angle. For this class of configurations, significant errors become prominent 
in linear theory and a nonlinear theory becomes necessary. Linear panel 
methods must use approximate techniques to compute bodies that have panels 
inclined greater than the Mach angle. Hence, significant departures from 
linear methods can be expected locally in regions of supercritical crossflow 
and more globally in regions where the bow shock position differs 
significantly from the Mach cone as is the case for wings or bodies with 
supersonic boundaries (i.e., normal Mach number greater than one). 

The primary intent of this report is to document and illustrate the use 
of the computer program NCOREL. The name NCOREL is an anacronym for Non- 
COnical RELaxation which is the basic technique used in NCOREL to obtain 
numerical solutions to the three dimensional nonconservative full potential 
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equation. NCOREL has in its ancestry tne COREL code which was developed for 
supersonic conical flows (Ref. 1,2). A computer program called SUNPOT for 
supersonic Nonlinear POTential flow was developed prior to NCOREL for the Air 
Force Wright Aeronautical Laboratory (see Ref. 3,4). This early version of 
NCOREL captured all shocks and was lacking in numerical efficiency compared to 
the present code. 

The success of COREL in its application to the design of a SC3 Conceptual 
Conical wing (i.e. , Supercritical Conical Camber) detailed in Ref. 5, paved 
the way for its extension to NCOREL. A landmark application of NCOREL to a 
three dimensional high lift maneuver wing (Ref. 6,7) demonstrated that for the 
first time aerodynami cists could reliably predict the detailed spanwise 
behavior of realistic supersonic wings in a high lift maneuver regime. 
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3 ANALYTICAL APPROACH 


The approach used in NCOREL consists of casting the nonconservative full 
potential equation in a spherical coordinate system. Conformal mappings are 
then used, if necessary, to generate a grid in the spherical cross flow plane 
capable of resolving typically large gradients in the vicinity of the wing 
leading edge. The three dimensional full potential equation written in a 
spherical coordinate system reduces identically to the conical two dimensional 
equation at r = 0. The primary advantage to this formulation is that a three 
dimensional computation can be automatically initiated with a conical flow at 
the apex using the same set of equations. Other supersonic marching codes 
either require a separate conical computation or another code to generate 
starting plane data. Unlike other supersonic methods utilizing nonlinear 
equations, NCOREL is a fully implicit marching technique which is not hindered 
by CFL constraints. The only constraint in the marching direction is that due 
to geometric accuracy. Hence, a three dimensional wing computation can be 
achieved in 10 to 50 steps depending upon accuracy requirements and geometric 
variation in the marching direction. 

The radial spherical coordinate is used as the marching direction. This 
is an absolute necessity for a method oriented toward wing flow fields at low 
supersonic Mach numbers. Other marching methods that utilize the axial 
cartesian coordinate as the marching direction very quickly encounter subsonic 
axial Mach numbers below Mach 2 for wings that are not highly swept. Explicit 
marching codes that use conformal mappings also typically encounter difficulty 
for wing geometries at low supersonic Mach numbers because the mesh size in 
the vicinity of the leading edge tends to diminish greatly. The vanishing 
mesh size and CFL condition in combination greatly restrict the step size in 
the marching direction, making the computation impractical. 

In NCOREL, finite difference solutions to the nonconservative full poten- 
tial equation are obtained at each spherical crossflow station commencing at 
the apex, R = 0, with a conical solution for the geometry obtained at R = AR. 
Transonic type dependent relaxation techniques form the basis for the 
numerical finite difference crossflow plane solutions. The current finite 
difference equations are second order accurate in elliptic crossflow regions 
and first order accurate in the hyperbolic crossflow region and spherical 
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marching direction, A rotated difference scheme (Ref. 8) is used in the 
supersonic crossflow region. The type dependency is coupled to the character 
of the crossflow Mach number on the spherical surface. The character of the 
crossflow terms switch from elliptic to hyperbolic as the crossflow Mach 
number exceeds sonic conditions. The numerical technique allows for the 
capture of both embedded crossflow and bow type shocks. Recently, in Ref. 9, 
the bow shock was fitted isentropically as a boundary. It was found that 
fitting the bow shock further increased the accuracy of the technique. As 
will be demonstrated in a later section, NCOREL has exhibited a high level of 
success for flows over isolated bodies and wings. 
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4 GENERAL METHODOLOGY 


4.1 TRANSFORMATIONS 

The method uses a series of three transformations from the physical 
cartesian (x,y,z) or spherical (o),4,r) coordinates to the computational evenly 
spaced rectangular crossflow mesh. The coordinate transformations are: 

I Spherical Coordinates (o),4,r) 

X = r cos 0 ) sin 

y = r sin w sin 4> (1) 

z = r cos 

Note that the Cartesian coordinate sytem used in NCOREL is atypical of the 
standard aircraft coordinate system in that the Z direction is axial and (X,Y) 
are the transverse plane coordinates. 

II Stereographic Projection (p,q,t) 

p = cos CO tan 4>/2 

a sin CO tan <\>/2 (2) 

t = r 

III Conformal Joukowski Mapping (p,0,R) 


S-S 
^ ^0 

S+I 

0 


r-s^/2 

= (^-) 

r+s^/2 


(3) 


where S = p + iq, r = pe^® and SQ(r) is the singularity location. 
IV Shearing Transformation (X,Y,Z) 


X = 0 

Y - 

^ ■ C-B 
Z = R = r 


(4) 
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where p = B(0,R) = Body Surface 

p = C(0,R) = Outer Boundary or Bow Shock Surface 
4.2 RELAXATION SCHEMES (SLOR) 

The three dimensional full potential equation in a spherical coordinate 
system can be represented by, 

Lc(4>) + RLnc^**^) " 0 (5) 

where is a conical operator and is a nonconi cal operator. At R = 0, 
the above equation reduces to the conical quasi two dimensional full potential 
equation. A perturbation potential is introduced as 

5 = q + ^ = V* + ■^ (6) 

where 

5 = uTg + tip + wTp 

This potential is further reduced to 

4. - RF(p,0,R) (7) 

The mapped space velocities are then defined as. 


U 


= -^ + U 
pH ° 


V = 



V 


W = F + RF + W 

P GO 

where U„, V„, and W„ are the freestream velocity components in the mapped 0, 

P, and r directions, respectively. H is the metric of the two conformal 
mapping transformations. 

If the Joukowski mapping singularity Sq is a function of r, the following 
mesh derivatives are generally nonzero and defined by, 

h^ = 5p/5r 
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h« ^ ae/ar 

2 P 

These mesh derivatives are used in the following transformation, 

ar aR ae ^2 ap 

The finite difference analog of the nonconservative full potential 
equation in the computational domain can be written as, 

- Bj) + (A2 - B2) + (A3 - B3) Fyy + (S) 

+ ®s''zz + ®9’"xz ®10’'yz + ®11 ''z = ° 

where the coefficients B^- = r¥^. vanish for conical flows. The X and Y terms 
are crossflow plane derivatives evaluated according to the type of crossflow, 
(i.e. , subsonic/elliptic or supersonic/hyperbol ic). The Z terms in the 
marching direction are always assumed to be hyperbolic. This implies the 
condition that the radial velocity must remain supersonic. The Z derivatives 
are evaluated using first order accurate upwind differences. The first order 
accurate second derivative F^^ requires that the potential at two previous 
spherical planes be stored. Hence, the potential at a total of three cross 
flow stations must be stored. 

Figure 1 shows the coordinate meshes in the mapped crossflow plane and in 
the sheared computational domain. The conformal mapping transforms the cross 
sectional projected geometry into a near circle which clusters grid points in 
the vicinity of the leading edge in the physical space. Two basic types of 
shocks can occur, the bow shock or Y-shock and the crossflow shock or X- 
shock. Since the freest ream must always be supersonic, a bow shock will 
always be present. If the bow shock is captured within the mesh (BSC method), 
the bow shock is defined by a crossflow sonic line. The crossflow velocity 
beyond the bow shock position will be supersonic and inside the bow shock only 
subsonic crossflow will exist if a crossflow shock or X-shock is not 
present. If the flow around the leading edge accelerates to supersonic 
velocities an embedded crossflow shock will also exist. For most flows, the Y 
or bow shock is due to the V component of velocity being negative (towards the 
body) and supersonic and for the X-shock, or crossflow shock, the U component 
of velocity becomes positive and supersonic if a > 0. 
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Successive Line Over Relaxation (SLOR) is used as the basic numerical 
procedure in the computer program NCOREL. In SLOR, the values of the 
potential on any given computational coordinate line X = constant or Y = 
constant are over-relaxed if the flow is subsonic and the equation is elliptic 
in nature. To determine the value of the correction to the potentials on a 
given line, a tridiagonal matrix must be solved using the appropriate boundary 
conditions as end conditions in a Gaussian elimination procedure. The 
efficient convergence of a SLOR scheme is best achieved when the sweep 
direction is not opposite to any supersonic zones. Hence, two sweep 
directions are then possible for the present problem as illustrated in Fig, 

1. In sweep A, values for the correction are determined on a Y = constant 
line (Ring) and the sweep direction is towards the surface or decreasing J. 

For sweep B, values for the correction are determined on an X = constant line 
(column) and the sweep direction is around the body towards increasing I or 
increasing 0. Both sweeps are in the direction of the supersonic zones if a > 
0. Each of these methods have been applied to the current problem. Sweep A is 
the original COREL method of Ref. 1 and 2. Sweep B was later implemented in 
Ref. 9. The test of which method is the best is how well they perform on 
problems of practical interest. Figure 2 shows a comparison of the two 
methods on some sample conical flow problems. The column method turns out to 
be two to five times faster than the ring relaxation method. The supersonic 
flow problem is largely determined by the formation of the bow shock. The 
column relaxation method inverts the tridiagonal matrix from the body to the 
outer boundary. Hence, the body boundary condition immediately affects or is 
felt by the bow shock. In the ring relaxation method, it takes approximately 
N iterations, where N is the number of grid points in the Y direction, for the 
outer boundary to feel the effect of the surface boundary condition. This is 
the major reason why the column sweep method is more efficient. There may be 
a region for a particular class of shapes or for a very high angle of attack 
where the ring relaxation method may be more efficient. For most practical 
cases, the Column Relaxation Method has been found to be more efficient for 
both conical and three dimensional flows and is used in the present NCOREL 
coFTputer program. 
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Figure 2 Convergence Rates for Conical Flow Solutions, 
Ring vs Column SLOR 
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4.3 BOUNDARY CONDITIONS 


For the symmetric half plane problem, periodic end conditions are used at 
the boundaries 1=2 and I = IC. This symmetry condition is enforced by 
storing the reflected potential at I = 1 and I = IC + 1. The values of the 
potential are set as: 

F(1,J) = F(3,J) 

F(IC+1,J) = F(IC-1,J) (9) 

If quarter plane symmetry exists at zero angle of attack, a reflected 
boundary condition is enforced at I = 1 + IC/2 and the sweep is conducted from 
1=1+ IC/2 to I = IC. 

The body or surface boundary condition of flow tangency or the vanishing 
of normal velocity is implemented by defining a dummy mesh ring inside the 
body at J = 1 (see Fig. 1). The boundary condition of zero normal velocity is 
then used to determine values of F (1,1) as a function of F (1,2) and F (1,3) 
and the body geometry by solving for the Y derivative of the potential. In 
this way, the surface boundary condition is automatically implemented when 
centered values of the derivative Fy are determined at the body surface J = 

2. The surface boundary condition is implemented for each iteration of the 
relaxation scheme until convergence. For the bow shock capture method (BSC), 
the value of the potential at the two outermost rings is set equal to zero. 
This assumes that the two outermost rings are at the freestream conditions, 
and an outer end condition of zero correction will be imposed in the 
tridiagonal solution. Hence, care must be taken so that the bow shock is 
captured within these coordinate rings. Generally, the bow shock should be 
captured with at least two rings beyond its positon where freestream 
supersonic crossflow velocities exist, not including the two rings held fixed; 
or J of the bow shock should be at least JC-3 to insure the proper solution. 
This occurs because the bow shock is captured across three mesh lines. One 
point will occur on the lower pressure or freestream side and another on the 
high pressure side of the bow shock with one point in between at approximately 
sonic conditions. If the outer boundary is chosen incorrectly, no supersonic 
crossflow velocities will occur within the two fixed rings at JC and JC-1, and 
erroneous results can occur. 
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In the Bow Shock Fit (BSF) Method, the boundary condition of zero 
potential is imposed only at the outermost ring. In this method, the 
outermost ring will be the bow shock, and in the absence of a crossflow shock, 
only subsonic crossflow velocities will exist for the entire mesh. The 
restrictions and requirements on the outer boundary that occur in the BSC 
Method are absent in the BSF method. The first derivative of the potential Fy 
is second order and taken one sided into the mesh or 






This derivative is necessary to determine the bow shock jump conditions. The 
different bow shock methods and their implementation will be explained more 
fully in a subsequent section. 


4.4 MESH REFINEMENT 

Mesh refinement is currently used in NCOREL to obtain the conical 
starting solution. A coarse mesh is used to initiate the solution with zero 
potential (i.e. , freestream conditions) used as starting potential values. 

This mesh is then halved to obtain a finer mesh. The potential solution from 
the coarser mesh is interpolated onto the finer mesh as a starting solution 
for the potential. Mesh refinement can speed up the computation by as much as 
a factor of four or five. Figure 3 shows a sample computation with and 
without mesh refinement. At little computational cost, a coarse mesh will 
very quickly establish the salient features of the flow field, such as the bow 
shock and embedded supersonic crossflow region. This flow field then becomes 
the starting guess for the next mesh. This is the primary reason for the 
enhanced computational efficiency. A total of two mesh refinements (i.e., or 
three grids) can be used in NCOREL for the conical solution. The nonconical 
solutions are then obtained by marching on the finest grid using the potential 
solution at the previous station as a starting solution for the new station. 


4.5 BOW SHOCK METHODS 

4. 5. 1 Bow Shock Capture (BSC) 

The capture of the bow shock and its position within the computational 
grid is largely determined by the prescription of the crossflow plane outer 
boundary. If the shape of the outer boundary is much different from that of 
the captured bow shock, the captured shock position in the crossflow plane may 
look somewhat like that depicted in Fig. 4a. The outer crossflow sonic line 
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Figure 3 Effect of Mesh Refinement 
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(i.e. , captured bow shock) can cut through several J rings or J lines in the 
computational grid. This is computationally inefficient in that many grid 
points of freestream conditions outside the bow shock will be carried along in 
the computation. To compensate for this and to better control the position of 
the bow shock, a method is used to align the bow shock with the computational 
grid. This procedure is beneficial not only in that the number of freestream 
mesh points beyond the bow shock can be controlled, but also in that it aligns 
the grid with the bow shock, which should yield a better captured shock 
solution and increased accuracy. 

Three bow shock capturing options have been included in NCOREL designated 
by the parameter lOPTS, 

Figure 4 illustrates each of these options. The first option (lOPTS = 0) 
is the default option. In this option the outer crossflow boundary is held 
constant and equal to a factor (EP - 1) times the rotated Mach cone for both 
conical and nonconical stations. Hence, the bow shock will typically cut 
through the cross flow computational mesh, as illustrated in Fig. 4a. For 
nonconical stations, the location of the bow shock will generally move within 
this mesh. Axially expanding bodies will cause the captured bow shock to move 
toward the surface (J = 2) and compression surfaces will cause the shock to 
move toward the outer boundary (J = JC). Hence, holding the outer boundary 
fixed can cause problems in that the movement of the bow shock has to be 
anticipated in positioning the outer boundary for the conical solution. If 
this is not done, the bow shock can hit the OLKter boundary causing erroneous 
results. In theory, the outer boundary could be tailored a priori to the 
particular body shape as a prevention technique. Instead, two other options 
were built into NCOREL. 

The first option (lOPTS = 1) is illustrated in Fig. 4b. This option will 
be initiated only at the conical station with a mesh refinement (KREF > 1). 

The first mesh will capture the bow shock as illustrated in Fig. 4a. The 
position and shape of the bow shock will be interpolated from the coarse 
grid. A new outer boundary will be automatically determined with the same 
crossflow shape as the interpolated bow shock. This shape will then be 
imposed NS grid lines or NS x AY beyond the interpolated bow shock sonic line 
to allow for proper capture. Hence, the parameter NS will approximately 
control the number of freestream mesh lines beyond the bow shock position. 
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Figure 4 Bow Shock Capture Methods (BSC) 


17 




The coarse mesh solution will then be interpolated onto this new finer mesh as 
a starting solution. The resulting solution on the finer mesh will look 
similar to that shown in Fig. 4b. The captured bow shock will now be aligned 
approximately with a Y grid line in the computational space. The number of Y 
grid lines beyond the shock is roughly determined by the parameter NS which 
shall be set > 4. 

The option lOPTS - 2 is similar to lOPTS = 1 but will occur at nonconi cal 
stations. This option, shown in Fig. 4c, not only aligns the captured bow 
shock with the cross flow grid, but also controls the vertical position of the 
bow shock by deforming the outer boundary as the marching solutions are 
obtained. Since mesh refinement is not used at nonconical stations, the 
previous station's interpolated bow shock shape and position will be imposed 
at the current station. This will prevent the captured bow shock from moving 
significantly through the mesh. 

For lOPTS = 0 and lOPTS = 1, the outer boundary will be held fixed in 
the nonconical computation. In lOPTS = 2, the outer boundary will move with 
the bow shock as the marching solutions are obtained. This causes R 
derivatives to be nonzero for the outer boundary. lOPTS = 0 is the fastest 
scheme and lOPTS = 2 requires the most computational time. The extra 
computational time is caused by the imposition of the interpolated solution 
onto a new mesh causing a somewhat slower convergence rate. 

A detailed study of the effect of each of these options on the solution 
has not been carried out. 

4.5.2 Bow Shock Fit (BSF) 

The bow shock fitting eliminates the need for capturing the bow shock 
within a prescribed outer boundary. Instead the outer boundary at J = JC 
becomes the bow shock. Its shape is governed by an isentropic shock jump 
condition, 



= Q, ) 

NORM p = c 


( 10 ) 
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where ^ is the density and insures continuity of mass flow across its 
boundary. In potential flow, the potential is continuous across a shock 
wave. Its normal derivative or normal velocity is discontinuous causing a 
pressure discontinuity or jump. The tangential derivative or velocity is also 
continuous across the shock wave. The reduced potential everywhere outside 
the bow shock in the cross flow plane is zero due to the freestream outer 
condition. Hence, its value at the bow shock boundary must also be zero. The 
value of the potential at the bow shock is then known, and its normal 
derivative is a function of the internal flow solution and the shape and 
location of the bow shock boundary. The normal derivative of the potential in 
Eq (10) is computed at the bow shock by a second order, one sided difference 
formula into the internal flow field. 

In the BSF method, the conical starting solution is by far the most 
difficult to determine because the bow shock boundary is not known a priori. 

An initial guess for its shape must be prescribed and iterated upon. At a 
nonconical station, the bow shock boundary is known and held fixed, and the 
slope of the bow shock in the marching direction is iterated upon. Iterating 
and perturbing the bow shock as a boundary for conical flow is also time 
consuming in a computational sense in that the entire cross flow mesh changes 
with the iteration process. As a result, the mapping metrics and freestream 
velocities must be recomputed after each iteration. In addition, the new bow 
shock shape must be underrelaxed so as not to disturb the convergence of the 
internal flow field. In a successful calculation, the internal flow residuals 
or errors associated with the potential flow Eq (8) and those associated with 
the shock jump condition Eq (10) will both diminish as the iteration procedure 
progresses. 

In conical flow, the isentropic shock jump condition Eq (10), can be 
written in a recursive tridiagonal form. 
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where the RHS of Eq (11) is determined using old values of the shock boundary 
C, and the LHS is then solved by Gaussian elimination for updated values of 
the shock boundary c\ Equation (11) was found to be the only stable method 
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to determine the bow shock boundary because it couples all the bow shock 
points. Periodic end conditions are used in the Gaussian elimination. The 
updated values of the bow shock position or shape are then under-relaxed. 

For nonconi cal flow, Eq (11) becomes 

C h C 

-pWRH = / (pU - U„)-pV + V„ + RH(h^ - -?^)(W„-pW)-W„RHC^ (12) 

In Eq (12), the slope of the bow shock location Cj^ at each point is determined 
as the internal flow is iterated upon. Because the crossflow mesh is held 
fixed, the nonconical calculation is not as sensitive, and the updated values 
of the shock slope in the marching direction can be over- relaxed. When 
both the internal flow converges and Eq (12) is satisfied, the new shock 
slopes are then used to determine a new shock boundary and mesh at the next 
nonconical station using the first order accurate expression: 

C.(R + AR) = C.(R) + Co AR (13) 

1 1 K. 

For the conical starting flow, an initial guess for the outer boundary 
must be prescribed that meets certain requirements. For subsonic leading edge 
crossflow shapes, the rotated Mach cone is usually a reasonable initial 
shape. Figure 5 illustrates some of the possible relationships between the 
rotated Mach cone and body shape in the cross flow plane. For other flows, 
the prescription of an initial guess is not as simple. The rotated Mach cone 
can not be used for near sonic or supersonic leading edge shapes since the 
Mach cone will touch or be within the boundary of the body surface (i.e. , 
usually around the leading edge). For body shapes, this condition can also be 
encountered at high angles of attack on the leeward surface. For other body 
shapes and freestream conditions where M^ sin is everywhere greater than 
one for the body, the Mach cone may lie entirely within the body. A crossflow 
mesh and a solution can only be generated when the initial guess is everywhere 
outside the body shape. Hence, the rotated Mach cone must be modified in such 
a way that an initial guess is prescribed that is everywhere within the 
boundary of the true bow shock and outside or coincident with the Mach wave. 

A solution may not be computed if any part of the initial guess lies within 
the rotated Mach cone or outside the true shock location. If a large part or 
all of the initial boundary lies outside the true shock location, the method 
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Figure 5 Mach cone Relationship to Type of Surface for Bow Shock Fit Method (BSF) 
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will attempt to capture the bow shock and is designed to fail to prevent this 
from happening. Since freest ream conditions or zero potential is used as 
initial values for the potential, the initial information generated at the 
prescribed bow shock boundary may be erroneous. This is another reason that 
under-relaxation must be used so as not to update the bow shock boundary too 
quickly causing an overshoot of the true bow shock position until the internal 
flow has somewhat converged. In a conical computation, supersonic crossflow 
velocities may occur within the boundary at a certain number of initial 
iterations. These supersonic crossflow points will vanish as the internal 
flow converges if a good shock guess has been prescribed, meaning that a 
compression is caused by the body which propagates toward the initial bow 
shock boundary and merges with it. A better way to initiate the conical BSF 
calculation might be to hold the boundary fixed (wg - 0.0) until the internal 
flow converges. This would cause the internal flow to converge, the 
compression to reach the boundary and supersonic crossflow points to 
disappear, and then allow for a higher bow shock relaxation value to be used 
because better information is obtained at the boundary. Since this is in 
general difficult to automate, NCOREL updates the boundary simultaneously with 
the internal flow iterations, and slightly smaller under-relaxation values 
must be used. The initial imposed shape for the conical bow shock on the 
crude grid can be checked by running NCOREL with zero for the bow shock 
relaxation parameter and then inspecting the crossflow velocities at the bow 
shock guess after a certain number of iterations (e. g. , 25-40). If the 
crossflow velocities are sonic or subsonic, then a good guess has been 
achieved for the bow shock. 

Low relaxation values (o)g < .10) are used for difficult conical BSF 
solutions. Difficult BSF solutions are those that fall into the Type II or 
Type III categories illustrated in Fig. 5. A crude grid (16 x 16) is used 
initially to determine the bow shock and solution. If mesh refinement is 
used, it is not essential for the crude grid bow shock jump conditions to 
converge totally, especially since a very low bow shock relaxation parameter 
may have to be used initially. 

Hence, given a reasonable initial guess, the bow shock fit method will 
iterate on the internal flow field and the bow shock shape simultaneously. 

The new shape of the bow shock at each iteration required to fulfill the shock 
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jump condition at each shock point will be under- relaxed. The solution is 
considered converged when both internal flow field residuals vanish or become 
small and the shock jump conditions are met at the shock boundary. 

Figure 6 illustrates the basic procedure for generating an initial guess 
for a bow shock boundary. The bow shock guess is generated by a combination 
of the rotated Mach cone and a normal displacement of the body surface 
boundary. The outermost boundary of these two shapes is then used as the 
initial shape. The displaced body boundary is controlled by an input 
parameter EP > 1. A factor (EP - 1.0) times the centerline half angle of the 
body is used as the basic displacement AS|_£ at the leading edge. So as not to 
yield a too discontinuous shape, the centerline displacement is taken to 
be the greater of the leading edge displacement or the displacement due to the 
Mach cone. The displacement at the leading edge and at the upper and lower 
centerlines are then faired linearly to determine the displaced body shape. 

If the Mach cone lies outside the displaced body, then that portion of the 
displaced body shape is replaced by the Mach cone. Other more clever bow 
shock guess generators may be found, but this procedure was found to be 
easiest to manipulate and the safest. 

Figure 7 illustrates several examples of actual conical flow BSF 
computations on a (16 x 16) grid for each of the type flow situations depicted 
in Fig. 5. Figure 7 shows the rotated Mach cone, initial shock guess and 
converged shock boundary. 

Mesh refinement is almost always used in conjunction with the conical BSF 
method since the first mesh solution is the most difficult to achieve because 
of the need for an initial shock guess. Refined meshes use the previous mesh 
solution for the bow shock boundary shape. 


23 



^ S, , = TAN ifJL + a)- TAN \I/. 

M D/ , _ IT 

CO -j 


As= DISPLACEMENT 

NORMAL TO SURFACE 


2 246-0 low 
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5 GEOMETRIC AND MAPPING REQUIREMENTS 


The full potential equation is solved in a spherical coordinate system. 
Hence, geometric cross sectional coordinates must be generated on successive 
spheres r = constant. A conical apex or nose cap is always assumed by 
defining the geometry between r = 0 and r = Ar (i.e. , Ar is the marching step) 
as a constant and equal to the cross sectional geometry at r = Ar. The 
geometry is defined one step ahead of the calculation. This is necessary to 
define the conical apex geometry and also to generate numerical second order 
accurate geometric derivatives (i.e., slopes, etc.). This technique 
inherently assumes that the geometry is smooth. Any discontinuities in the 
geometry will be smoothed to the order of the marching step size. At each 
spherical station, NG cross sectional points will be determined from the 
geometry routines. At a given r station, the leading edge coordinates are 
first determined in order to scale the points. Typically the NG geometry 
points will be clustered in the vicinity of the leading edge for better 
accuracy by a cosine or sine variation. This set of geometric points usually 
exceeds the number of points in the actual cross sectional final computational 
mesh for accuracy. The geometry points will then be transformed to the mapped 
space given a singularity location for the particular cross section. Periodic 
cubic spline coefficients are then generated for the NG mapped space 
coordinates. The actual IC mesh points, at equally spaced 9, are then found 
by interpolation from the NG input geometry points in the mapped space. A 
transformation of these interpolated mapped body points at equal 9 back to the 
physical space will then generate a mesh where points will be clustered in the 
vicinity of the wing leading edge. All geometric or surface derivatives are 
generated numerically in the mapped space. 

5.1 MAPPING REQUIREMENTS 

The mapping transformation will, in general, be a function of r if the 
geometric cross sections are nonconical and change in character as a function 
of r. The functional dependence on r comes about from the r variation of the 
singularity location. This variation of the mapping transformation with r 
results in mesh derivatives defined as h]^ = Pp and h 2 = pOr. If the geometry 
is conical, or if the singularity location is fixed with r for nonconical 
flows, these mesh derivatives will vanish or h^^ = h 2 = 0. 
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In general, the best choice for the singularity location is the center of 

radius of curvature of the wing leading edge. In practice, using a fine mesh, 

the solution is somewhat sensitive to the precise placement of the singularity 

location. The exception to this would be the special case of sharp (e. g. , 

small leading edge radius of curvature wings.) For this special case, the 

singularity must be placed at the center of radius of curvature. Hence, the 

singularity location is determined by the following approximate procedure. 

The leading edge coordinates are first determined (x^^-p, yfip» ^tip^‘ 

centerline half thickness at x = 0 is determined. The lateral or x location 

of the singularity is then determined as the focus of an ellipse or 
2 2 2 

X . = Y^. - Y The y location (ycin) oT the singularity is further 

sin tip x==0 '''Sin' 

determined to be on the mean line at the location Given the singularity 

location (x^^-^, yg-jn). the mapping transformation is then defined and the set 
of NG input geometric points can be transformed to the mapped space. At each 
r-station, the cross sectional geometry and singularity location are also 
determined at r + Ar. This is necessary because the r derivative of the 
singularity location is required to determine the mesh derivatives h^ and 
h 2 * Numerical second order accurate singularity derivatives are computed at r 
+ Ar, requiring the location of the singularity at the upstream station 
(r + 2Ar). 

The premise for geometry generation in formulating NCOREL was that users, 
in general, will not have spherical geometry at their disposal. Hence, the 
conputer code has been tailored to the availability of transverse (z = 
constant) cross sectional geometry. The code then internally iterates for the 
spherical geometry given transverse geometry data. The procedure is 
illustrated in Fig. 8a. An initial iteration takes place to determine the 
leading edge coordinates at a given r-station. Since the lateral singularity 
location is determined approximately as the focus of an equivalent ellipse, an 
additional approximation is made by determining the centerline thickness at z 
= z^^p and not on the spherical surface. The singularity location is then 
determined solely by the transverse data at z = z^-jp. The premise for this 
approximation is that the cross section at r = constant will not differ 
dramatically from the cross sectional characteristic shape at z = z^ip* 
may not always be the case and a more exact placement of the singularity may 
be necessary. This approximation, though, greatly simplifies the generation 
of the singularity coordinates. 
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Figure 8 Singularity Placement and CON USE Logic Used in NCOREL 
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5.2 GEOMETRY OPTIONS 

Several geometry options have been built into NCOREL based on the 
availability of transverse geometric cross section data being as follows. 

IQPTG - 1: (Axi symmetric bodies of revolution) 

A routine is supplied that will generate cross sectional coordinates 
given the radius of the body as a function of z. 

IQPTG = 2: (Analytical Geometry Package) 

A geometry definition based on simple conical parameters was developed 
and was found to be very useful in defining both conical and three dimensional 
wings, bodies, and wing-body configurations. The cross section was simply 
defined as the superposition of a spanwise thickness and camber distribution. 
The basic wing cross section is restricted to spanwise elliptical thickness 
and circular arc camber. 

In this analytical representation of the wing geometry, the geometric 
parameters are specified as parametric angles. Thus, a wing cross section can 
be entirely specified by five angles: 6p, 6^, 6j, 6p. The symbols 6p, 

and 6j represent planform angles defining the body lines of the wing and 6p is 
an angle controlling the centerline thickness variation, and 6 jq the leading 
edge camber angle In the cross sectional plane (see Fig. 9). Cubic equations 
are used to axially curve fit the angles in the general form: 

6(z) = a(z - z^*)^ + - z-j)^ + c(z - z^.) + d 

The spatial location of the body lines can then be simply defined as: 

X = z tan (6(z) ) 

The thickness of the wing was assumed elliptical but was allowed to have 
an inboard flat section. The camber line, defined separately, was also 
allowed to have a flat inboard section with a circular arc camber line to the 
leading edge. To define the maximum camber and twist of the cross section, a 
leading edge camber angle is specified (6-j-q). The thickness and camber flats 
are defined separately and, thus, do not have to coincide. The wing cross 
sectional geometry was then simply defined by adding and subtracting the 
elliptical thickness to the camber line. 
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The present capability of this geometry package allows for a transition 
from a cone to a conically cambered wing section. A further deformation from 
a conical wing section to a conical wing-body cross section can also be 
generated. Figure 10 illustrates the basic capability of the geometry 
package. This package can be used to generate conical wing or body or wing- 
body cross sections for conical solutions. It can also be tailored to 
generate three dimensional deforming geometries of the above types by a 
judicious choice of input parameters. Variable sweep wings can also be 
generated with this definition. Figure 11 shows some sample geometries that 
can be generated. Figure 11a shows a typical body line deformation to a 
conical delta wing. Figure 11b shows a cross sectional deformation from a 
circular cone to an elliptical cross section. Figure 11c shows the 
deformation from an elliptical cross section to a cambered cross section. The 
last section of the geometry can deform into a wing-body cross section defined 
as shown in Fig. 12a. The last conical section of the geometry requires one 
additional angle for the body definition. If the body angle is made equal to 
the wing thickness angle the cross section will degenerate into a wing cross 
section. Figure 12b shows an example of the geometry that can be generated. 

LOPTG - 3: (Analytical User Supplied Geometry) 

This geometry package is supplied by the user in terms of a subroutine. 
The basic requirements of the subroutine are illustrated in Fig. 8 and solely 
requires that the transverse geometry be defined (i.e z = constant). Figure 
8a illustrates the use of the subroutine in finding the leading edge coordi- 
nates and singularity location for a call to the subroutine with IFLAG = 0. 
Figure 8b illustrates the general usage of the subroutine for defining the 
cross sectional points. The code will then use this information to generate 
the spherical cross sectional geometry. Note that the user is not restricted 
as to the placement of the singularity as the focus of an ellipse. 

lOPTG = 4: (Bi -cubic Surface Patch Geometry) 

This is a somewhat more flexible geometry package although a much more 
conputational ly time consuming package. This geometry package uses the Harris 
Wave Drag Input where a wing is defined by airfoil sections whose thickness 
and camber are specified as a set of input points. A body or fuselage is 
defined as a set of transverse cross sections where the coordinates of each 
cross section are specified as input. The geometry package is an adaptation 
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Figure 12 Conical Wing-Body Geometry Generation 
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of the Craidon bi-cubic surface patch geometry program of Ref. 10. It was 
chosen because the Harris Wave Drag Input format is widely used by aerodynam- 
i cists and by NASA in defining geometries and also is compatible with many 
linear panel codes. The bi-cubic surface patch definition was chosen because 
a continuously defined and smooth geometry must be available to NCOREL. The 
geometry package as adapted to NCOREL can generate arbitrary wings, bodies, or 
wing-body geometries based on point wise input. It treats wing and body 
components individually and superimposes these cross sections to determine a 
wing-body cross section. This geometry package also works well with the 
constraint that the individual wing or body geometries are smooth. A detailed 
description of the basic geometry package and mathematics of the bi-cubic 
surface patches and spline fitting can be found in Ref. 10. Essentially, a 
set of input points are converted to quadrilaterals. Each quadrilateral is 
then defined to be a bi-cubic surface patch defined by a 4 x 3 boundary 
matrix. Since a large number of patches (e. g. , 2000) can be defined for a 
vehicle, all of the patch matrices can not be stored. The entire set of patch 
matrices are written on disk and read when needed. In the adaptation of the 
patch geometry to NCOREL, it was found necessary to store the patches required 
to define the geometry for a particular spherical cut. The I/O time required 
to read and sort through the entire set of patch matrices for one cross 
sectional point was found to be prohibitive. Hence, for each spherical cut 
the patch matrices are read and sorted only once and the required patches 
stored. This increases the core requirements of NCOREL but reduces the I/O 
related CPU time enormously. Since the patches are defined parametrically by 
two cubics in U and W coordinates, a Newton iteration in two variables has to 
be used in order to determine an arbitrary surface coordinate. This is also 
costly in computational time and causes the patch geometry to be the slowest 
geometry option in NCOREL but the most user oriented. 

It should be noted that the Craidon patch geometry routines use the 
standard aircraft coordinate system denoted by (X^, Zy^) where Xy^ is axial, 

Yy^ is spanwise and Zy^ is vertical, whereas NCOREL uses a nonstandard 

coordinate system (Xq, Y^^, Zq) where Z^^ is axial, Xq is spanwise and Yq is 
vertical. These coordinate systems are referenced to and used interchangeably 
in NCOREL. 
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In general, any geometry package that can conform to the simple format 
illustrated in Fig. 8 can be coupled to NCOREL, 
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6 AERODYNAMIC FORCES AND MOMENTS 


To tailor NCOREL to the specific needs of aerodynami cists, a detailed 
knowledge of forces and moments generated by a configuration is essential. 
Since NCOREL computes spherical solutions, an adaptation is required to 
generate forces and moments that are useful to the engineer. To do this, an 
interpolation scheme is used that first recomputes all of the AR spherical 
surface data to Az spanwise (z = constant) data. The interpolated data is 
then used to compute axial incremental forces and moments on each spanwise 
segment Az. These incremental spanwise forces and moments are then integrated 
to yield a running account of the axial force and moment generation. To 
accomodate wings with variable or constant sweep trailing edges, a special 
interpolation scheme was developed to determine the forces and moments aft of 
the centerline trailing edge station. Figure 13 illustrates the basic 
technique. The z-station interpolation proceeds with the same step size as 
the spherical surface data until the trailing edge station. At the centerline 
trailing edge, z = ZTEC, spanwise data is interpolated. 

For all of these stations, each radial grid line is used for inter- 
polation to z = constant stations. Hence, IC-1 mesh lines are interpolated to 
yield IC-1 interpolated points for a z = constant station. Aft of the 
centerline trailing edge station, z = ZTEC, the z location of the interpolated 
stations is governed by the intersection of a radial grid line with the 
trailing edge. The actual interpolated z values and number of trailing edge 
stations is governed by the number of grid lines used in the computaton. 

Hence, for each of the interpolated z stations aft of the trailing edge the 
number of interpolated values will diminish by two. Finally, at z = ZTIP, 
only one value will be interpolated and the interpolation scheme will 
terminate. If the computation to R = ZTIP has not been completed the 
interpolation scheme will terminate when the interpolation from x = x-p^ to x = 
^LE completed. The interpolated z-station data for the surface 

coordinates and pressure coefficients is then written onto TAPE2. 

The trailing edge computation is contingent upon a user supplied 
subroutine XTES which yields a ZTE value given XTE. If a delta wing or body 
is being computed, the last interpolated station occurs at a z = constant 
station specified in the input. 
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Figure 13 General Interpolation Scheme for Arbitrary Planforms 
Used for Aerodynamic Computation 
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After the interpolation has been completed, a force subroutine is called 
which reads the interpolated z = constant data from TAPE2. The wing or body 
is now divided into a set of quadri 1 ateral panels where the four pressure 
coefficients are known at each corner point. For variable sweep trailing 
edges, the spanwise sections contain quadrilaterials and two triangular panels 
which include the trailing edge points. To compute the force and moment on 
each quadrilateral, the quadrilateral is divided into four different 
triangular planes or two pairs of triangles governed by the two quadrilateral 
diagonals. The normal, axial force, moment and areas are then computed for 
each pair of triangular planes and averaged. This is done because if only one 
set of triangular planes was used, the projection of an upper surface 
quadri 1 ateral onto a lower surface quadrilateral will result in different 
triangular planes for the two elements. Hence, for a symmetrical wing, exact 
concellation of normal forces might not occur, resulting in a net small but 
finite normal force on the vehicle. By taking both sets of triangular planes, 
an exact cancellation of these normal forces will occur, if quarter plane 
symmetry exists in the geometry and pressure coefficents, yielding a more 
exact computation. The pressure acting at the centroid of each panel is 
assumed to be one third of the pressure at the corner points. 

The forces, moments and areas are summed spanwise at each z - station. 
Projected planform area, surface area and cross sectional areas are all 
computed. The user can monitor the forces, moments and areas generated by 
each z increment. These incremental values are then integrated yielding an 
axial account of the total forces, moments and areas. 

Normal and axial forces are computed as well as normal and axial force 
coefficients, lift and drag coefficients, and moment coefficient. 

Figure 14 illustrates several examples of both incremental and integrated 
aerodynamic output by NCOREL. The first portion of the figure illustrates the 
planform and surface area of the configuration nondimensional ized by the 
maximum planform area; the surface area is always slightly greater than twice 
the planform area depending upon the cross sectional character. The second 
portion of the figure shows the cross sectional area nondimensional ized by its 
maximum value. The remaining two portions of each figure show lift- 
coefficient and drag-coefficient information. Figure 14a shows results for a 
simple symmetrical delta wing with a subsonic leading edge, elliptical cross 
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section and a parabolic centerline thickness. Figure 14b shows results for a 
symmetrical, supersonic-leading-edge, arrow wing with a 4% thick, four digit 
NACA airfoil section and a highly swept trailing edge. Downstream of the 
centerline trailing edge of the arrow wing, the sectional lift coefficients 
increase significantly, causing an overall increase in integrated lift 
coefficient. The axial variation of drag coefficient is distinctly different 
for the two wings. Figure 14c shows a typical output for a wing with variable 
sweep leading and trailing edges. The sectional spanwise lift coefficients 
also increase significantly aft of the centerline trailing edge. Figure 14d 
shows the aerodynamic output for a lifting missile body with an elliptical 3 
to 1 cross section and a Haack-Adams area distribution at high a. The forward 
half of the body has high lift characteristics. The decreasing lifting 
characteristics are due to the expansion of the body. Figure 15 shows the 
moment coefficients for the four different configurations. The aerodynamic 
output represented in Fig. 14 and 15 can be very useful to the designer in 
tailoring the behavior of a particular configuration. 
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Figure 15 NCOREL Aerodynamic Output, Moment Coefficients 
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7 GEOMETRIC RESTRICTIONS 


Several geometric restrictions currently exist in NCOREL. NCOREL was 
originally developed for the computation of three dimensional flow over smooth 
bodies and wings. Any abrupt or discontinous charge in geometry will result 
in a discontinuous expansion or embedded shock. Currently NCOREL captures all 
embedded shocks. Theoretically, it can compute embedded supersonic shocks. 
This shock is different from the embedded cross flow shock because the flow 
across an embedded cross flow shock changes from supersonic/hyperbolic to 
subsonic/elliptic. This characteristic change does not occur across a 
supersonic shock. The interaction of embedded supersonic shocks with the 
fitted or captured bow shock has not been tested in NCOREL. A supersonic 
shock can be generated by an abrupt change in leading edge sweep of a wing or 
body angle. The occurrence of this shock in combination with a supercri tical 
cross flow region and shock also has not been tested in NCOREL. NCOREL will 
capture these supersonic shocks but will have a tendency to smear them across 
several mesh points, making their accuracy questionable. 

Conical flow computations are limited to attached flow conditions where 
the bow shock remains attached to the apex of the surface. Detachment can 
occur for a given Mach number at a high angle of attack or at a given Mach 
number as the local surface deflection angle is increased. Typically, just 
prior to detachment, subsonic total Mach numbers will be found in the conical 
solution. These solutions can not be used as a starting solution for a 
marching computation. 

High angle of attack computations over t hin wing cross sections may fail 

under certain conditions. If the cross section is thin and as the angle of 

attack is increased, a small region of near limiting suction pressure 
2 

(i.e., C ~ ~~~ 9 ) appear on the leeward surface. Further increase in 

*^1 i m 

the angle of attack may cause the computation to diverge or fail. Two 
explanations are possible; either not enough resolution or mesh points are 
available to resolve the large gradients in the vicinity of the leading edge 
or an inviscid solution does not exist. Sometimes very thin conical wings 
must be started on finer grids in order to converge and a large coefficient of 
artificial viscosity must be used. If these two devices fail, then the above 
explanations probably hold. 
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NCOREL computations are also limited to detached leading edge bow 
shocks. The bow shock will become attached to the leading edge for supersonic 
leading edges with zero radius of curvature. This condition can not be 
computed with NCOREL because a computational grid between body and bow shock 
can not be generated at the leading edge. Although this exact condition can 
not be computed, wings with small leading edge radii can be computed which 
approach this condition but still have the bow shock at a small but finite 
standoff distance. 

In its application to wings, only wings with supersonic trailing edges 
can be treated at the present time. An embedded supersonic shock will be 
generated at the trailing edge of this wing (i.e. , Mach number normal to 
trailing edge remains supersonic). An approximation for the wake is made by 
extending the camber line of the wing into the centerline inboard of the wing 
trailing edge. The wake surface is then treated as a flat plate with the 
boundary condition of flow tangency. This results in erroneous wake pressures 
that are not matched but should not affect the wing pressures if the wing has 
a supersonic trailing edge. This approximation becomes exact for a 
symmetrical wing at zero angle of attack. Because the actual wing trailing 
edge point is not demarcated in the cross flow plane, some smearing of this 
shock and a slight upstream effect on the wing will occur in the solution. 

This behavior will be corrected and an extension to subsonic trailing edges 
will be included in a future version of NCOREL which treats the wake properly. 

Another area of difficulty in wing applications is the planform of a 
clipped delta or trapezoidal wing where the abrupt change in leading edge 
angle has been smoothed. The difficulty occurs when the flow expands around 
the leading edge. At zero angle of attack, the flow must expand through an 
angle equal to the planform angle of the wing. This usually does not present 
a problem to NCOREL. A realistic wing will typically approach zero thickness 
in this region of the wing. At angle of attack, the flow approaching this 
region may already have a low pressure supercritical region on its leeward 
surface. The flow then has to be expanded further through the half planform 
angle of the wing. NCOREL does not always yield a solution for this type of 
geometry in the clipped region of the wing. Limiting pressures may be 
obtained associated with large crossflow Mach numbers. It is suspected that in 
some cases an inviscid solution to this flow does not exist and the only way 
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that a solution can be obtained by NCOREL is to modify the leading edge 
geometry making the expansion less severe. This adverse behavior is a 
function of the planform angle or sweep, thickness and camber in the vicinity 
of the expansion and the angle of attack and may not always occur. 

Wing-body cross sections can be treated with NCOREL with no apparent 
difficulties as long as the necessary mesh resolution is available. 

Difficulty will be encountered in treating a wing body planform where the wing 
emerges abruptly and at a position significantly downstream of the apex of the 
vehicle. In this situation, a strong embedded shock will be generated by the 
wing emergence. NCOREL will capture this shock but will have difficulty 
resolving the flow field if the initial cross section of the wing is thin and 
a supercritical region has developed. An insufficient number of mesh points 
may lie between the captured wing shock and wing surface, causing convergence 
difficulties. A blended wing-body configuration with strake blending is most 
suitable to be treated by NCOREL. 

Since it is impossible to evaluate every application that NCOREL could 
possibly be applied to, it is left up to the user to evaluate its capability 
when going beyond the limit of isolated smooth bodies and wings. In applying 
NCOREL to a particular configuration, the user should keep in mind the 
possible areas of difficulty mentioned above. 
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8 APPLICATIONS 


NCOREL has been applied successfully to a variety of conical flows and 
three dimensional isolated bodies and wings. Several simple wing-body 
configurations have also been treated successfully. The wing-body problem 
still needs further investigation, and the user should keep this in mind in 
attempting this type of computation. 

8.1 CONICAL FLOWS 

NCOREL can treat almost any kind of conical body, wing or wing-body 
configuration. Conical flows can be treated using any of the geometry 
packages or a user supplied routine where the cross section is invariant with 
the axial coordinate. This can be achieved analytically or with a simple 
nondimensional point wise input. The first station solution at R = 0 is 
always treated as a conical solution. The only inherent limitation in NCOREL 
is that the bow shock remain attached at the apex. The typical solution 
preceding detachment will contain subsonic total Mach numbers. NCOREL will 
fail at some point when the angle of attack, Mach number, or surface angle 
exceeds the condition of an attached bow shock. 

NCOREL will also fail for a thin wing at some high angle of attack. 
Preceding the failure or divergence of a solution, a region of limiting 
suction pressure or near vacuum pressures will typically occur on the leeward 
surface near the leading edge. The divergence of the solution indicates one 
of two things: either a lack of resolution (i.e. , need for a finer mesh) or 

that an inviscid solution can not be achieved under those flow conditions. 

Figure 16a shows an example of some typical conical surface pressure 
coefficient distributions obtained by NCOREL on a very fine mesh (86 x 58) at 
M(a = 2,0 and angles of attack up to 15° for a thin elliptic cone. A cross 
flow shock has formed at a = 5° and gets stronger and moves inboard with 
increasing a. At a = 15°, a small portion of the suction plateau is close to 
limiting or vacuum pressure. Higher angles of attack could not be computed 
because of the presence of these limiting pressures that correspond to very 
high cross flow Mach numbers. Figure 16b shows a similar set of surface 
pressure distributions for a wedge wing cross section with a small leading 
edge radius of curvature. At a = 5°, a large cross flow shock appears on the 
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Figure 16 Conical Surface Pressure Distributions for Two Dissimilar Cross Sections = 2.0, 0^ = 20 ) 




leeward surface. At a = 10°, a small leeward region of near limiting 
pressures has already occurred. Higher angles of attack could not be computed 
even with this fine mesh. Figure 16b shows clearly why leading edge 
separation occurs early on wings with small leading edge radii. 

Figures 17 and 18 show a variety of interesting bow shock fit conical 
flows that can be computed with NCOREL. These figures show the cross section, 
cross flow sonic line, and cross flow streamline pattern. Supersonic cross 
flow velocities exist everywhere within the sonic line which is usually 
terminated by a cross flow shock. The cross flow streamlines emanate from the 
bow shock (i.e. , the outer boundary line) and terminate at the body surface 
either at a vortical singularity or cross flow saddle point all of which are 
cross flow stagnation points (i.e., U = V = 0), Winglike cross sections will 
have three stagnation points, two of which will be vortical singularities 
located in the leeward and windward symmetry planes. The third cross flow 
stagnation point will occur near the leading edge of the wing. Body like 
cross sections will not have the leading edge stagnation point. 

Figure 17a shows the solution for a thin supersonic elliptic cone at zero 
incidence showing the formation of a supersonic cross flow zone on both 
leeward and windward surfaces. The leading edge bow shock is also approaching 
an attachment condition. Hence, supersonic cross flow regions can appear on 
both upper and lower surfaces depending upon the freestream and geometric 
conditions. An interesting phenomenon can occur on cambered wings at low 
angle of attack. Figure 17b shows a highly cambered wing at zero incidence. 

In this flow a supersonic bubble and cross flow shock occurs on the lower 
surface of the wing. 

Figure 17c shows a subsonic cone at high incidence. In this flow the 

leeward vortical singularity lifts off the body surface. Figure 17d shows an 

interesting flow over a subsonic elliptic cone where a large supersonic cross 

flow bubble exists on the leeward surface terminating in a very strong cross 
flow shock. Another sonic line exists on the windward surface. This is a 
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Figure 17 Conical Bow Shock Fit Crossflow Streamline Patterns 


Figures 18 a, b and c show the solutions for a thin elliptic cone at 
incidence and Mach numbers corresponding to subsonic, near sonic and 
supersonic leading edges. The cross flow region becomes larger and the cross 
flow shock moves inboard with increasing Mach number until finally the cross 
flow region extends all the way to the bow shock. 

Figure 19 shows two conical wing-body flow solutions. Figure 19a shows a 
very blended wing-body cross section, and Fig. 19b shows a more abrupt wing- 
body juncture. It is interesting to note that an additional stagnation point 
occurs in Fig. 19b at the windward wing-body juncture. The windward vortical 
singularity moves from the symmetry plane to the wing-body juncture. 

Figure 20a shows spanwise experimental pressure distributions for a 
conical cambered wing designed with the use of COREL (see Ref. 2). Figure 20a 
also shows the computed NCOREL solutions for both the BSC and BSF methods. 

Both methods compare well at the near sonic condition M^ = 1.66. At the 
higher Mach number, supersonic leading edge condition, the BSF method yields 
better windward correlation with the pressure data. Figure 20b shows the 
corresponding differences between the computed bow shock locations for the two 
methods. 

In addition to the cambered wing that was designed for a shockless 
recompression at M„ = 1.62 and C^ - 0.40, the equivalent flat or uncambered 
wing was tested. Figure 21 shows the cambered and flat wing cross sections at 
the design lift condition of C|_ ~ 0.40. Excellent correlation and validation 
of the method was achieved at the design condition for the cambered wing, 
indicating a virtually shockless recompression. In comparison, the flat wing 
shows a higher suction pressure plateau and stronger cross flow shock. The 
comparison with NCOREL is good everywhere except in the vicinity of the cross 
flow shock. Some boundary layer separation may be occurring that is not 
accounted for in the inviscid NCOREL solutions. 

8.2 NONCONICAL FLOWS 

8. 2. 1 Analytic Body Shapes 

A series of analytically defined forebodies was tested at supersonic 
speeds ranging from M^^ = 1.70 to 4.50, and the pressure data are available in 
Ref. 11. Two of these body shapes were chosen to determine the accuracy of 
NCOREL with respect to surface pressure data at M<„ = 1.70. All the 
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Figure 18 Conical Bow Shock Fit Crossflow Streamline Patterns for a Thin Elliptic Cone, 6 ^ = 20 °, 5^ = 2 °. and a = io*^ 
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Figure 19 Conical Bow Shock Fit Crossflow Streamline Patterns for a Wing-Body Cross Section, Mqo= 2.0, a = io°. 6 ^ = 20 ° 
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(A) SPANWISE PRESSURE DISTRIBUTIONS (B) COMPUTED BOW SHOCK LOCATIONS 

2246-042W 

Figure 20 Conically Cambered Wing 
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conputations were carried out on a (38 x 30 x 50) mesh where 38 points were 
distributed around the body and 50 radial marching steps or Z/L = 0.02 were 
used. The conformal mapping was not implemented for these body shapes. The 
total run times ranged from 4 to 7 min on an IBM 370/168 with the spherical 
geometry generation requiring approximately 2 - 3 min of the total run time. 

Figure 22a and 22b show the computed windward and leeward symmetry plane 
pressures compared to the test data at = 1.70 for a parabolic arc body with 
a lobed cross section at a = 5.35° and 10.35°, respectively. The test data do 
not exhibit a smooth behavior due to either probe or geometric inaccuracies. 
Excellent overall comparison is achieved at these angles of attack. Figure 23 
shows a comparison of computed azimuthal pressure distributions with data for 
the same body at = 1.70, a = 10.35° for several axial locations. Excellent 
agreement in the azimuthal shape of the pressure distribution is achieved with 
a slight difference in the magnitude of the pressure levels. 

Figure 24 shows another set of comparisons for a fuselage forebody shape 
at = 1.70, a = 1.00° and a = 5.01°. According to Ref. 11, the first five 
pressure taps were incorrectly located and hence do not necessarily lie on the 
symmetry planes and should be disregarded. Figure 25 shows the azimuthal 
pressure distributions for the same body at = 1.70, a = 1.00° and several 
axial locations. 

Overall excellent agreement with the test data was achieved. 

8.2.2 M issile Bodies 

Experimental force and moment coefficients for two missile bodies are 
reported in Ref. 12. Circular and 3 to 1 elliptical cross section body shapes 
were developed from a minimum drag Haack-Adams area distribution. NCOREL was 
used to compute the body alone characteristics of these two shapes in Ref. 

13. Figure 26a shows the moment and normal force coefficients as a function 
of angle of attack for the ciruclar cross section body. A linear Woodward 
panel method calculation and experimental data points are also shown. NCOREL 
slightly overpredicts the moment and underpredicts the normal force. It is 
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Figure 22 Axial Pressure Distribution on a Parabolic Arc Body 
with Lobed Cross Sections at = 1.70 
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These curves dramatically demonstrate the need for a nonlinear method, 
especially at high angles of attack. The discrepancy between the experimental 
data and NCOREL is probably due to viscous effects that are entirely neglected 
in NCOREL. 

8. 2. 3 Analytic Arrow Wing 

Reference 14 reports on a series of wind tunnel tests for four 
analytically defined arrow wings. This experiment was selected because of the 
availability of detailed spanwise data at several stations for a three- 
dimensional configuration. The basic symmetrical wing configuration referred 
to as Model 1 with a sweep angle A = 63.43° was used for the present study 
(see Fig. 27). The analytical geometry was kindly supplied by G. Moretti of 
PINY in the form of a computer subroutine. 

The geometry consists of a wing with a straight leading edge and a 
circular centerbody emerging from the wing downstream of the apex. The 
variable sweep trailing edge of the wing remains a supersonic trailing edge 
for all Mach numbers tested. The geometry was run in NCOREL at the lowest 
test Mach number of 2.36. A (57 x 41) cross flow plane grid and 40 radial 
steps were used for the computation. 

Figure 28 shows the correlation achieved for a = 6° and span stations 
Z/L = 0.30, 0.50, 0.65, and 0.80, The computed results at a spanwise station 
were obtained by interpolation of the results from the spherical surfaces. 

The body has already emerged from the wing at the first span station (Fig. 
28a), where Z/L = 0.30. The grid definition of the body is adequate at this 
station. The computed pressures do not indicate any localized spanwise effect 
due to the wing-body juncture, although the test data show a slight variation 
in the vicinity of the juncture. A cross flow shock is not indicated in 
either the experimental data or the conputed results at this station. Figures 
28b and 28c show two further downstream stations. The formation of a weak 

rrncc: f 1 nw chnrif ic nnw PviHpnt* At thpcp «t;itinnc thp nrirl Hpfinition nf thp 

body is rather poor although the correlation in pressures remains excellent. 
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Figure 27 Arrow-wing Wind-tunnel Model Showing Pressure Orifice Locations <Ref. 14) 
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Figure 28 Spanwi$e Pressure Distribution on an Analyticaiiy Defined Arrow Wing, Moo = 2.36, a ~ 6"" 





by the wake approximation. Remarkably good comparison is achieved on the wing 
surface using this approximation. The computed results indicate trailing edge 
shocks that occur due to the supersonic nature of the trailing edge. The 
local solution for the shocks in the vicinity of the trailing edge is not good 
because of the erroneous flat plate approximation as well as boundary layer 
effects and artificial smearing of the oblique trailing edge shocks. Smearing 
of the trailing edge shocks occur because grid points are not necessarily 
located on the exact trailing edge of the wing and also due to the use of 
numerical central difference formulae for the computation of body spanwise 
slopes. The computed body pressures do not correlate with the test data. 

This was expected because the wake of the wing has a definite influence on the 
body, unlike on the wing. Some difficulty in the geometry of the flat plate 
wake occurred as indicated in Fig. 28d by the wake crossover. 

The computed results have a tendency to slightly underpredict the 
supercritical pressures and indicate a slightly downstream position of the 
cross flow shock in comparison with the test data. This discrepancy may in 
part be due to the rotational effects associated with the high Mach number of 
2.36 which are neglected in the present potential theory. 

Higher angles of attack were not computed because the experimental 
pressure distributions indicated a flat leeward pressure distribution with the 
absence of a cross flow shock. This is generally indicative of leeward 
surface boundary layer separation, and correlation with inviscid solutions 
would not be expected, 

8.2.4 Supersonic Cruise Aircraft 

A more stringent test of the code in computing wing-body geometries is 
the type of configuration where the body origin lies upstream of the wing 
origin, which is typical of any realistic aircraft. This type of geometry 
leads to the formation of an embedded oblique wing shock. Reference 15 
contains spanwise surface pressure data for a supersonic cruise aircraft 
configuration at Mach numbers 2.30, 2.96, and 3.30. Reference 15 was chosen 
because it also contains detailed measured geometric data of the complete 
aircraft in the form of an input data set consistent with the patch geometry 
requirement of Ref. 10. 
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Figure 29 shows the panel definition of the aircraft and the location of 
some of the pressure stations used for comparison defined in aircraft model 
coordinates (in.)* The geometry data set of Ref. 15 was modified for use in 
NCOREL. The engine nacelles and vertical tails were removed. The wing model 
contains a total of 650 panels and the fuselage model, which was made simpler 
than that of Ref. 15, contains a total of 328 panels. The Joukowski mapping 
singularity was then tied to the evolution of the wing. Hence, the mapping 
singularity remains on the axis of the aircraft until the wing origin. Then 
the singularity moves continuously outward with the wing as the wing emerges 
from the fuselage. This procedure was found to be necessary because abrupt 
changes in singularity location cause abrupt changes in the mesh, leading to 
erroneous results. The camber of the wing centerline section was also 
modified to make the apex of the wing coincide with the origin of the fuselage 
so that the mapping singularity smoothly departs from the axis of the 
fuselage. 

Some preliminary computations were carried out at = 2.30. A (58 x 58) 
cross sectional mesh was used with a marching step size of AZ = 2 which 
yielded approximately 23 steps to the centerline trailing edge of the wing. 
Figure 30 shows an isometric plot of the spanwise cross sections and pressure 
distributions at = 2.30, a = 4.82°. As noted earlier, the numerical 
program marches on spherical surfaces using the spherical cross sectional 
geometry. Figure 30 was obtained from the spherical data by linear 
interpolation. A large, low pressure spike occurs as the wing emerges from 
the fuselage (x ~ 14). The spike remains throughout and a cross flow shock 
quickly develops. Interesting to note is the pressure depression that occurs 
due to the fuselage. Just after the wing emerges, the pressure distribution 
exhibits large oscillations on the leeward surface, causing convergence 
problems. This may be related to oscillations in the geometry or to the 
simultaneous development of the oblique embedded wing shock due to the wing 
emergence from the fuselage and the appearance of the supercritical cross flow 
shock. 

If a thin wing (i.e. , in comparison to the fuselage) emerges from the 
fuselage, a new scale to the problem is created. In the wing-body 
calculation, only a few points will be distributed between the newly emerged 
wing surface and oblique wing shock embedded in the fuselage flow field. As 
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the oblique wing shock moves away from the wing surface towards the bow shock, 
the resolution becomes better. It is suspected that this is the primary 
difficulty being encountered in the vicinity of the wing emergence from the 
fuselage. The lack of resolution can create oscillations and convergence 
problems especially if a supercritical region has developed. A stretching of 
the coordinates creating a cluster near the newly emerged wing surface may be 
necessary to alleviate this problem. 

The computed results show a large suction spike at the leading edge with 
a steep adverse pressure gradient. The computed solution indicates that the 
suction spike induces a cross flow shock. The flow then expands slightly to 
further recompress through a second cross flow shock. This unusual behavior 
is most probably caused by the small nose radius and rather large camber. The 
suction spike would most probably cause the leading edge to separate. Figure 
31 shows a series of comparisons at = 2.30 and a = 4.82° between a wing- 
body and wing alone calculation with the corresponding test data. The 
presence of the body has a minor effect on the windward pressures and a 
significant but local effect on the leeward surface, yielding a compression 
from the wing-body juncture to the centerline. Both wing and wing-body 
calculations exhibit leading edge suction spikes, with the wing-body 
exhibiting a more inward location and a slightly stronger cross flow shock. 
Overall the two calculations compare favorably with only minor and expected 
differences. The leeward surface experimental pressure data do not indicate a 
suction spike or a cross flow shock. The absence of a cross flow shock in the 
experimental data as the wing becomes thinner must be due to leading edge 
boundary layer separation. 

The lack of correlation in the leeward surface pressures led to a closer 
inspection of the test data. Two intermediate stations where only leeward 
pressure data points were taken are plotted for comparison in Fig. 32a and 32b 
for a = 2.81° and a - 4.82° without the wing alone solutions. Except for the 
leading edge pressure spikes remarkably good correlation was achieved with the 
leeward experimental and theoretical pressures at these two intermediate 
stations. Interesting to note is that at the two intermediate stations the 
experimental pressures exhibit cross flow shocks at the higher angle of attack 
(a = 4.82°). Also, if one compares the stations of Fig. 32 with the stations 
shown in Fig. 31, the experimental leeward data of the two intermediate 
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stations are distinctly different from their neighboring stations which raises 
some concern about inconsistencies in the experimental data. 

8. 2. 5 Zero Lift Wave Drag of a Delta Wing 

Figure 33 shows a series of numerical computations at several Mach 
numbers for a three dimensional delta wing (A = 71.57°) at zero angle of 
attack. The Weber version (see Ref. 16) of the Squire delta wing has a 
circular arc centerline thickness distribution with elliptic spanwise cross 
sections. To compute the entire wing in a spherical coordinate system, a flat 
plate was inserted for the wake geometry. This is a valid approximation since 
the wing is symmetric and at zero angle of attack. The computed pressures on 
the wing were integrated to obtain the zero lift wave drag values at several 
different Mach numbers and compared to experimental values and the Squire 
linearized analysis. This computation was suggested by Mason because it had 
been noted earlier in Ref 17 that linearized theory has a tendency to 
overpredict the drag. Hence, the opportunity exists of achieving lower levels 
of drag using a more accurate nonlinear analysis tool. The nonlinear code 
shows a marked improvement above Mach 1.8 for the difficult prediction of zero 
lift wave drag. The discrepancies at the lower Mach numbers are being 
investigated. At the Mach numbers 1.5 and 1.6, the linearized results show 
slightly better agreement. The accuracy of the nonlinear analysis depends 
upon the accuracy of capturing the oblique trailing edge shocks. To better 
resolve these shocks, a windward differencing was used to determine the 
numerical surface slopes in the radial marching direction. Less smearing of 
these shocks occurred when this was done and slightly better correlation was 
achieved. No attempts were made in the analysis to actually delineate the 
trailing edge of the wing. Aft of the trailing edge, the spherical cross 
sections take on a lobed appearance outboard of the wake region. A grid point 
on the body surface did not necessarily fall exactly on the trailing edge of 
the wing. Also, centered differences were used to determine the 
circumferential body slopes. All of this tends to artificially smear the 
oblique trailing edge shock. The surface pressures on the wing should not be 
influenced by the trailing edge shocks, but, unfortunately, the above 
numerical approximations that tend to smear the trailing edge affects the wing 
pressures slightly. 
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Figure 33 Computed and Measured Zero Lift Wave Drag on the 
Squire Delta Wing 
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8. 2. 6 NASA/Grumman Demonstration Wing 


A three dimensional supersonic wing was designed by Grumman Engineering 
for NASA Langley and tested in the Langley Unitary Plan Wind Tunnel. This 
wing was designed with the aid of NCOREL at a design condition of = 1,62 
and a = 12° with a lift coefficient Cj_ ~ 0.40. One of the objectives of this 
program was to use a nonlinear theory to design a realistic wing and achieve 
attached flow at the design condition by minimizing the strength of the cross 
flow shock. A post test evaluation of the wind tunnel results against the 
pretest predictions of NCOREL also was used as a validation of the accuracy of 
the code (see Ref, 6 and 7). 

Four spanwise stations, x = 10.6, 15.50, 19.90 and 24.40, were instrumen- 
ted for pressures. The last station occurs roughly at the centerline trailing 
edge of the wing. The results obtained with NCOREL (BSF) are for a (57 x 57) 
crossflow mesh and a marching step size of one inch. This yields approxi- 
mately 30 steps to the tip section of the wing. The NCOREL solutions are 
obtained on spherical surfaces. The code computes, by linear interpolation, 
both spanwise and chordwise pressure distributions. Figure 34 shows a com- 
parison between the predictions of NCOREL and wind tunnel data at = 1.62, a 
= 9.92°, At the first station, x = 10.60, a slight shift in pressure levels 
is apparent between theory and experiment. The tunnel data indicates a higher 
pressure level on both windward and leeward surfaces. At x = 15,50, the 
windward and the leeward supercritical pressures are in excellent agreement. 
The experiment indicates a stronger compression and higher subcritical 
pressures on the leeward surface. This trend persists at the next two 
stations. At this condition, a distinct cross flow shock has not yet 
developed except possibly at the last station, x = 24,40, It is interesting 
to note that NCOREL has no difficulties computing down to the last station 
where the wing is very thin and almost a flat plate in the vicinity of the 
centerbody. With this grid the body definition is not very good at the last 
two stations. The centerbody causes an expansion to occur. Unfortunately, 
the wind tunnel model was not instrumented in the vicinity of the center- 
body. Figure 35 shows the design condition of M„ = 1.62, a = 11.93°. The 
same trends between NCOREL and experiment persist, but a cross flow shock 
develops at x - 19.90. The experiment seems to indicate a slightly stronger 
cross flow shock but excellent agreement in shock position is achieved. At 
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Figure 34 NCOREL Comparisons With DEMO Wing Surface Pressure Data, Moo 1.62, a = 9.92° 
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Figure 3^5. NCOREL Comparisons with DEMO Wing Surface Pressure Data, M„ = 1.62, a = 11.93 
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the last station, x = 24.40, the experiment does not indicate the plateau in 
supercritical pressure that NCOREL predicts. The last plot. Fig. 36 shows the 
conparison at = 1.62, and a = 13.92 which corresponds to a ~ 0.47. In 
this figure, a crossflow shock develops at x = 10.60. A distinct shock is 
apparent in the test data at x = 19.90. NCOREL overpredicts the supercritical 
plateau at x = 24.40. 

A comparison between NCOREL and test data at = 2.00, a = 11.80° is 

shown in Fig. 37. The test data represent the model with the alternate 
leading edge, and the NCOREL predictions are for the basic leading edge. The 
difference beween the two leading edges are minor, with the alternate leading 
edge having slightly less camber very close to the leading edge. The 
agreement is excellent at all four stations. Overall, the agreement between 
the tunnel data and NCOREL is quite good. 

Figure 38 shows a series of isobar plots for the demonstration wing. The 
planform plotted is that computed by NCOREL where a small portion of the tip 
could not be computed. All the plots were generated by interpolating for 100 
pressure coefficient levels between the limit of -0.5 and + 0.5. The dark 
regions indicate large pressure gradients. This occurs typically near the 
apex, leading edges and in the vicinity of crossflow shocks. On the leeward 
surface at a = 7.92°, a small compression occurs as the wing gets thinner 
toward the tip section. The crossflow compression develops with angle of 
attack and moves inboard. Finally, at a = 13.92*^, a distinct cross flow shock 
is indicated in the isobar pattern. The leeward isobars become more parallel 
to the freestream with angle of attack. The windward surface isobar pattern 
does not change that dramatically with angle of attack. The compression due 
to the centerbody is quite apparent. An interesting feature is that the 
windward isobars tend to become parallel to the leading edge with angle of 
attack. The compression that occurs in the vicinity of the trailing edge is 
due to the trailing edge shock smearing. 

As an interesting comparison, the equivalent flat wing was computed at 
approximately the same lift, i.e. , = 1.62, a = 9.40°. The flat wing has 

the same planform and thickness distribution of the DEMO wing. Figure 39 
shows a comparison of the two isobar patterns. A very strong crossflow shock 
appears on the flat wing in comparison to the DEMO wing. 




78 




2246-039W 

Figure 36 NCOREL Comparisons With DEMO Wing Surface Pressure Data, Moo - 1.62, a = 13.92° 
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Figure 37 NCOREL Comparisons With DEMO Wing Surface Pressure Data, Moo = 2.00, a 
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9 COMPUTER PROGRAM DESCRIPTION 


9.1 ORGANIZATIONAL STRUCTURE 

The NCOREL computer program in its segmented form requires about 215 Kg, 
on a CDC computer. Segmentation was chosen over Overlay to reduce core 
requirements so that a more portable computer code would be created that could 
be run on either CDC or IBM machines. Segmentation, unlike Overlay, does not 
require changes in the source code to reduce core storage. Instead, a 
segmentation directives file is created that controls the sequential loading 
of subroutines in order to reduce core. A segmentation tree based on these 
directives is created that controls the sequential loading of subroutines in 
order to reduce core. The basic core requirements of the code is governed by 
the main level common blocks and subroutines and the branch segment or 
sublevel of the tree that requires the most core. Any number of other 
branches can be included in the logic without increasing core as long as the 
core of these branches are less than the core of the largest branch. 

Figure 40 illustrates the basic subroutine dependence used in NCOREL and 
its branches or sublevels. Each level may have other subroutines included or 
secondary branches using other subroutines. The name of a branch or sublevel 
coincides with a subroutine name. The INCLUDE statement identifies other 
subroutines that are to be loaded simultaneously with the routine that the 
branch name identifies. Figure 41 further illustrates the breakdown and 
calling sequence of the somewhat complicated INITN branch. Figure 42 shows 
the segmentation directives for NCOREL. The GLOBAL statements include the 
labeled COMMON BLOCKS within a particular level. For example, the labeled 
common block named BLKl is GLOBAL to the main level NCOREL. Hence, BLKl is 
always stored in the main level and this information is available to any of 
the sublevels. In addition, the loading of any sublevel or branch does not 
destroy this information. This is true of all the labeled commmon blocks 
GLOBAL to the main level NCOREL. Any other common block information that is 
only specific to a sublevel is made GLOBAL to that sublevel. When that 
sublevel or branch is returned to main and another sublevel loaded, the 
information is lost or written over by the new sublevel being loaded. The 
CONUSE sublevel is a user supplied geometry package that needs a CONUSE 
subroutine. Any other subroutines that are required in CONUSE must be named 
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(lOPTG = 2) ■ IXWI > 0 
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Figure 41 Subroutine Dependence in INITN Branch 
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Figure 42 Segmentation Directives for NCOREL 



by the user in the INCLUDE statement. Figure 43 shows the segmentation tree 
diagram that is generated for a typical run. The only variation to the 
segmentation will occur in the INCLUDE statement for the user supplied CONUSE 
geometry. Segmentation reduces the core requirement by a factor of two to 
three. If large core is available, NCOREL can be run without segmentation as, 
for example, on a IBM 3033. Any small subroutines that are called throughout 
NCOREL have been included in the main level. 

Figure 43 is also annotated to show the core required for each branch of 
the segmentation. The main program requires about 111 Kg. The relaxation 
routines that generate the cross flow solutions require about 164 Kg. The 
next largest branch, 174 Kg, involves the aerodynamic subroutines that 
interpolate and compute forces and moments. The largest branch of 212 Kg 
comes from the patch geometry subroutines. Hence, the patch geometry branch 
GEOM sets the core requirements for NCOREL. 

9.2 SUBROUTINE DESCRIPTIONS 

Table I groups the subroutines as to their function in NCOREL. 

MAIN PROGRAM NCOREL 

This is the main program that reads and writes an echo of the input 
control data from TAPE 1. In addition, it controls the calling sequence of 
the sublevels. Any geometry initialization and reading of geometric data from 
TAPE 5 and storing of input parameters occurs within the main program. 
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Figure 43 Segmentation Tree Diagram for NCOREL 
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Table 1. Subroutine Groups 


GROUP I 

(MAPPING, RELAXATION SOLUTION. AND SPHERICAL OUTPUT) 


GROUP lA 


GROUP IB 


NCOREL (MAIN) 


(AUXILIARY PROGRAMS OF 


GROUP lA) 


INITN - GEOMETRY MAPPING 




FLSOLV 


PCS 

PERIODIC CUBIC 

STRECN - METRICS 


DELTA 

INTERP 

■ SPLINE 

REL CApI RELAXATION 



RELFIT r SOLUTION 


CONF 
RCONF _ 

^ WING MAPPING 

OUTPN - SPHERICAL OUTPUT 
RENEW- REFINES 


SHOCKI - 

-SHOCK GUESS 

MARCHING STEP 





GROUP tl 

(ANALYTICAL OR USER 
SUPPLIED GEOMETRY) 

GEOMIA- READ INPUT 
GEOMA - SPHERICAL CROSS 
CROSS SECTION 
SINGA - SINGULARITY 
LOCATION 

CONBOD - AXISYMMETRtC 
BODY 

CONWIN - ANALYTIC WINGS 
CONUSE - USER SUPPLIED 


GROUP m 

(HARRIS INPUT AND PATCH GEOMETRY) 


GROUP IMA 
GEOMIN 

GROUP IIIB 

GROUP MIC 

START - READS INPUT 





GEOM - SPHERICAL CROSS 

XCUT 

TRANSVERSE 

VSOLV 

BI-CUBIC 

SECTION 

XCUTB 

- SURFACE 

VSOLVB 

.PATCH EQUATION 

SINGUL - SINGULARITY 

ZCUTB J 

COORDINATES 

TSOLV 

UWSOLV, 

SOLVER 

SURF ”1 

PACH _r GENERATION 

TIPLE 

SPHERICAL 



TIPB 

’ LEADING EDGE 



SPFIT - SPLINE FITTING 
PACHXY - PATCH CORNER 

YMBOD 




POINTS 





XPATCH - SPHERICAL CUT 





PATCHES 






2246-003W 
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Purpose 


Subroutine Name 


GEOMIA 


GEOMIN 


(START) 


Reads geometric parameters from TAPES and initializes and 
stores data used in the supplied analytical geometry 
package (lOPTG = 2). 

Reads Harris Wave Drag Input from TAPES and initializes 
data for Craidon Bi -Cubic Patch Geometry Option (lOPTG = 
4). 

Reads in Harris Wave Drag Input from TAPES and writes 
surface coordinate data on TAPEIO. 


(SURF ) 


Reads surface coordinate data from TAPEIO and controls 
writing of patches on TAPE 14. 


(PACH ) 

(SPFIT) 

(PACHXY) 


Computes and writes Patch coefficient matrix on TAPE14. 

Spline fits geometric surface data for patch definition. 

Reads PATCH matrices from TAPE14 and computes and stores 
maximum and minimum coordinates of BODY and WING patches 


INITN 


Computes spherical geometry and mapped space body and shock 
information for relaxation solution 


(GEQM ) 

(XPATCH ) 

(SINGUL) 

(GEOIA) 

(SINGA) 

(SHOCKI ) 

(CONF ) 

(PCS ) 

(DELTA) 

(INTERP ) 


Computes spherical geometry and singularity location at 
r = r« + Ar where r is the current solution station (lOPTG 
= 4).° 

Reads and sorts entire set of patches and stores patches 
required for spherical cut (lOPTG = 4). 

Computes leading edge and singularity location at 
r = r^ + 2Ar (lOPTG = 4). 

Computes spherical geometry and singularity at 
r = r^ + Ar. (lOPTG = 1, 2 or 3) 

Computes leading edge coordinates and singularity location 
at r = r^ + 2Ar. (lOPTG = 1, 2 or 3) 

Computes spherical coordinates of initial shock shape guess 
(BSF) or outer boundary (BSC) 

Maps input body or shock physical spherical coordinates to 
mapped space. 

Computes derivatives used in finding spline coefficients 
Determines spline fit parameters 

Interpolates from periodic cubic splines to find body or 
shock mapped space radii at mesh points. 
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(RENEW) Recomputes body geometry when marching step size is 

refined. 

FLSOLV Controls the calling of subroutines STRECN, RELCAP, RELFIT, 

and OUTPN 


(STRECN) Computes mapping metrics, nonconical mesh derivatives, and 

freest ream velocity components. Inverts mapping to obtain 
physical spherical mesh coordinates. 

(RELCAP) Computes spherical cross flow plane bow shock captured 

(BSC) relaxation solution for the reduced potential 
F(I,J) when IFIT = 0. 

(RELFIT) Computes spherical cross flow plane bow shock fit (BSF) 

relaxation solution for the reduced potential F(I,J) when 
IFIT = 1. Also recoRputes mapped space bow shock 
coordinates C(I), its derivatives, and mesh metrics. 

(OUTPN) Outputs spherical cross flow plane solutions on TAPE6 and 

interpolates for cross flow sonic line and bow shock if 
IFIT = 0. Also writes spherical solutions on TAPE7. 

(RCONF) Inverts mapping to obtain physical spherical coordinates 

given mapped space coordinates 


NOTE : Two Cartesian coordinate systems are referenced and treated 

interchangeably in the geometry routines. Zq = Yq - Zq = Xa 
where the subscript C refers to the cartesian coordinates associated 
with the spherical coordinate system and the subscript A refers to the 
standard aircraft coordinate system. 


SPANI 


(INTERZ) 

CHORDI 


(STREAM) 


AEROF 

(INTERN) 

(FORCE) 


Controls interpolation for transverse cross sectional data 
(Zq = = constant) by reading arbitrary transverse 

locations specified as input on TAPEl. 

Interpolates for a given ^ = constant station from 

spherical solutions read from TAPE?. Interpolated 
transverse coordinates and pressure data written on TAPEIO. 

Controls interpolation for chordwise data (Xq = Y^ = 
constant) by reading for arbitrary chordwise locations 
specified as input on TAPEl. 

Interpolates for Xq = Ya = constant station from spherical 
data read from TAPE?, writes interpolated coordinates and 
pressure data on TAPE14. 

Controls the computation of forces, moments, and areas. 

Interpolates and writes on TAPE2 transverse plane data used 
in force calculations from spherical data read from TAPE?. 

Computes forces, moments, and areas from set of 
quadrilateral panels determined by INTERN 
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(PLANE) Computes area, normal and axial force for a given triangular plane 


XTES 

GEQMA 

SINGA 

lOPTG =1 

lOPTG = 2 
lOPTG = 3 

GEQM 

(TIPB ) 

(TIPLE) 

(XCUTB) 

(XCUT ) 

(TSOLV ) 


A subroutine used by CONUSE, INTERZ, STREAM, and INTERN to 
yield the or given X^ or 

ANALYTICAL GEOMETRY SUBROUTINES 
(TOPTG = 1, 2, or 3) 

Iterates to find spherical geometry at r = r^ + Ar where r^ 
is current solution station (lOPTG = 1, 2,or^3) 

Computes leading edge coordinates and singularity location 
at r = r^ + 2Ar. (lOPTG = 1, 2, or 3) 

(CONBOD) Given axial station or Xa computes radius of 
axi symmetric body of revolution, (NOTE: equation of body 

radius R = R(Zc) is supplied by user) 

(CONWIN) A nalytical geometry package for conical and 
nonconical wings, bodies or wing-bodies. Given axial 
station Zq or X^, and spanwise coordinate Xq or Y^, 
computes surface coordinate Yq or Z^. 

(CONUSE)User supplied geometry package simulates use of 
CONWIN 

CRAIDON BI-CUBIC PATCH SUBROUTINES 

(iTptg = 4) 

Computes spherical surface geometry and singularity 
location from Harris Wave Drag Input and Bi-Cubic Surface 
Patches at r = r^ + Ar. 

Used in iteration procedure to find spherical maximum half- 
width of BODY cross section at given axial location. Given 
(Z, X) or (X,Y)a it computes Y or Zn of the body. Uses 
XCUTB, ZCUTB and \MB0D. 

Used in iteration procedure to find spherical wing leading 
edge coordinates. Given (Z, X) or (X,Y)^ it computes the 
surface coordinate Y^ or Z« of the leading edge. Uses 
TSOLV and XCUT. 

Used in iteration procedure to determine spherical BODY 
surface coordinates. Given (Z,X)^ or (X,Y)^^ it computes Y^ 
or Z^. 

Used in iteration procedure to determine spherical WING 
surface coordinates. Given (Z,X)^ or (X,Y)^ it computes Y^ 
or Z^. 

Used by TIPLE to compute WING leading edge coordinate from 
leading edge patches. Given U^^ and (Z,X)^ or (X,Y)^ it 
determines W by iteration to find Y^ or Z^. 
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(VSOLV) Used by XCUT and XCUTB, a Newton iteration method for 

finding U and W on a patch corresponding to given values of 
(Z,X)j, or (X,Y)y\ yielding the surface coordinate or Z^. 

(ZCUTB) Used by TIPB to find X^ or Yn surface coordinate of a BODY 

patch given the coordinates U,Y)^. or (X,Z)^. 

(VSOLVB) Newton iteration method for finding U and W on a BODY patch 

given (Z,Y)^. or (X,Z)^ yielding the surface coordinate Y^. 
or Z^. 

(YMBOD) Determines maximum and minimum coordinates of a set of BODY 

patches. Used by TIPB 

9.3 COMMON BLOCKS 


The following is a partial description of the labeled common blocks used 
in the basic NCOREL subroutines. Labeled common blocks associated with the 
various geometry options are excluded. 


Label Contents 

/BLKl/ LIN, LOUT, IC, JC , KT, DX, DY, WW, GA, EMINF, A02, 

GAMMA, RELW, DELMIN, PI, PI02, ALP, IPLOT, EST, NOUT, JOUT, 
IRUN, IPUN, REND, ISYM 

Variable Definition 


LIN = 1 

LOUT = 6 

IC 

JC 

KT 

DX 

DY 

WW 

GA 

EMINF 

A02 

GAMMA 

RELW 

DELMIN 

PI 

PI02 

ALP 

IPLOT 

EST 

NOUT 

IRUN 

IPUN 

REND 

ISYM 


Input device for namelist data 
Output print device 

Initial number of mesh points in X or 9 direction 
Initial number of mesh points in Y or p direction 
Maximum number of iterations in the relaxation routines 
Mesh size in computational X or 9 direction 
Mesh size in computational Y or p direction 
Over-relaxation parameter 
Specific heat ratio parameter 
Freestream Mach number 
Stagnation speed of sound squared. 

Ratio of specific heats 
Bow shock fit relaxation parameter 
Convergence tolerance on potential solution 
% 

nJ2 

Angle of attack (radians) 

Plotting parameter 

Coefficient of temporal artificial viscosity 

Print control parameter 

Restart control parameter 

Punch control parameter for restarting 

Spherical termination station 

Symmetry control parameter 
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Label 


Contents 


/GEO/ 

NG, NS, AZ, BZ. NBOD 


Variable 

Definition 


NG 

NS 

AZ 

BZ 

NBOD 

Number of spherical cross sectional geometry points 

Mesh control parameter for bow shock capture 

Tangent of leading edge planform half angle used in patch 

geometry 

Not used 

Number of spherical body cross sectional geometry points 

Label 

Contents 


/SURVEY/ 

KIT, IDIV 


Variable 

Definition 


KIT 

IDIV 

Number of iterations for convergence at an r- 
IDIV = 0, spherical station has converged 
IDIV = 1, spherical station has diverged 

-station 

Label 

Contents 


/TE/ 

IZW. ITE, ZTE, XTEBK, ZTIP 


Variable 

Definition 


IZW 

ITE 

ZTE 

XTEBK 

ZTIP 

Print control for interpolated Z - stations 
Triggers call of trailing edge subroutine XTES 
Default z-station for terminating force calculation 
Spanwise location of break in trailing edge 
Axial terminaton of wing tip 

Label 

Contents 


/AEROC/ 

CBAR, YCG, ZCG 


Variable 

Definition 


CBAR 

YCG 

ZCG 

Aerodynamic chord used in moment coefficient 
Vertical location of moment center 
Axial location of moment center 

definition 

Label 

Contents 


/RAD/ 

R, DZ, KR, KRAD, lOPTG, lOPTS, EP, IFIT 
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Variable 


Definition 


R 

Spherical r - station 

DZ 

Step size in marching direction 

KR 

Radial station integer 

KR = 1, conical station (R = 0) 

KR > 1, nonconical stations (R > 0) 

KRAD 

Marching step size refinement parameter 

lOPTG 

Geometry option parameter 

I OPTS 

Bow shock capture option parameter 

EP 

Bow shock guess parameter 

IFIT 

BSC or BSF option parameter 

Label 

Contents 


/MAP/ SO, SI, S2, PII, CSCAL, SOPR, SIPR, SON, SIN, S2N, 

SOPRN SIPRN, S2PRN 


Note: Complex numbers except for PII and CSCAL 


Variable 


Definition 


SO 

Singularity location 

SI 

Complex conjugate of 

S2 

Singularity shift in 

PII 

11 , constant 

CSCAL = 1.0 

Not used 

SOPR 

Radial derivative of 

SIPR 

Radial derivative of 

S2PR 

Radial derivative of 

SON 

Singularity location 

SIN 

£omplex conjugate of 

S2N 

S , at r = r + Ar 
Singularity shift in 

SOPRN 

Radial derivative of 

SIPRN 

Radial derivative of 

S2PRN 

Radial derivative of 

Label 

/BOWSHK/ 

THEB(201), RH0B(201) 


in mapped space, S^, ajt r 
singularity location, S , at r 
mapped space, at r 


mapped singularity, S^p at r 
conjugate, S at r 
singularity shift, S 2 p at r 

in mapped space, S^, at r + Ar 
singularity location, 

mapped space, at r + Ar 
mapped singj^l arity, S^p at r + Ar 
conjugate, S at r + Ar 
singularity shift, S£p at r + Ar 

Contents 

THES(201), RH0S(201) 


Variable 


Definition 


THEB 

RHOB 

THES 

RHOS 


Body theta locations in mapped space from input 
geometry cross sectional points 

Body radii in mapped space from input geometry cross 
sectional points 

Bow shock theta locations in mapped space 
Bow shock radii in mapped space 
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Label 


Contents 


/PLO/ 

XXB (60), YYB (60) 


Variable 

Definition 


XXB 

Physical body spherical coordinates 
points 

(X/Z) at mesh 

YYB 

Physical body spherical coordinates 
points 

(Y/Z) at mesh 

Label 

Contents 


/INT/ 

L(201), D(201) 


Variable 

Definition 


L,D 

Used in periodic spline fit routines 



Label 

/TIP/ 

Variable 

XTIP 

YTIP 

XSIN 

YSIN 


Contents 


XTIP, YTIP, XSIN, YSIN 

Definition 

Leading edge physical coordinate x 
Leading edge physical coordinate y 
Singularity physical coordinate x 
Singularity physical coordinate y 


Note: Labeled blocks /BLK2/, 

common blocks. 


Label 


/BLK2/ 


Variable 

F(I,J) 

H(I.O) 

X(I) 

Y(I) 

B(I) 

BPR(I) 

BSEC(I) 

C(I) 

CPR(I) 

CSEC(l) 

Hl(I.O) 


/BLK2P/, /BLK3/, and /BLK4/ are large core 


Contents 


F(60,60), H(60,60), X(60), B(60), BPR{60) BSEC(60) 
C(60), CPR(60), CSEC (60), Hl(60,60), H2(60,60) 

Definition 


Reduced potential at R 
Mapping metrics for mesh at R 
Sheared computational coordinate 
Sheared computational coordinate 
Body mapped radius 
Mapped body derivative, B 0 
Mapped body derivative, Bgg 

Bow shock (BSF) or Outer Boundary (BSC) mapped radii 
Mapped bow shock/boundary derivative, Cg 
Mapped bow shock/boundary derivative, ^09 
Noncomcal mapped mesh derivative, p» 

Nonconical mapped mesh derivative, p6r 
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Label Contents 

/BLK2P/ FP(60,60), FPP(60,60), BR(60), BRR(60), BRTH(60), 

CR(60), CRR(60), CRTH(60), BN(60), CN(60), CP(60), 
CPP(60) 


Variable 


Definition 


FP(I,J) 

FPP(I,J) 

BR(I) 

BRR(I) 

BRTH(I) 

CR(I) 

CRR(I) 

CRTH(I) 

BN(I) 

CN(I) 

CP(I) 

CPP(I) 


Reduced potential at R - AR 
Reduced potential at R - 2AR 
Nonconi cal mapped body derivative Bj^ 

Nonconical mapped body derivative B|^p 

Nonconical mapped body derivative 

Nonconical mapped bow shock/boundary derivative 

Nonconical mapped bow shock/boundary derivative Cj^|^ 

Nonconical mapped bow shock/boundary derivative ^R0 

Body mapped radius at R + AR 

Bow shock/boundary radius at R + AR 

Bow shock/boundary radius at R - AR 

Bow shock/boundary radius at R - 2AR 


Label Contents 

/BLK3/ VI(60,60), UI(60,60), WI(60,60) 


Variable 


Definition 


VI(I,J) Freestream velocity component in mapped p direction 

UI(I,J) Freestream velocity component in mapped 0 direction 

WI(I,J) Freestream velocity component in spherical r direction 


Label Contents 

/BLK4/ HP(60,60,), H1P(60,60), H2P(60,60) 


Vari abl e Def i ni ti on 

HP(I,J) Mapping metrics at R - AR 

H1P(I,J) Nonconical mapped mesh derivative at R - AR 

H2P(I,J) Nonconical mapped mesh derivative at R - AR 


9.4 I/O UNIT DEVICES 


The following is a description of the I/O units utilized in NCOREL and 
their purpose. 


Unit 


Usage 


TAPEl Namelist input device 

Input for arbitrary spanwise and chordwise 
interpolation 

TAPE2 Used in force computations (contains z 

station data) 
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TAPE5 

TAPE6 

TAPE7 


Input device for geometry data 
Print device for output 
Spherical data output device (surface 
coordinates, pressure and bow shock data) 
tapes Plow field data output device for restarting NCOREL at 

R > 0. Also, flow field data input device for restart. 
TAPEIO Used in patch geometry definition; Chordwise 

interpolated data output device 
TAPE14 Used in patch geometry definition; Spanwise 

interpolated data output device 

Table 2 shows the I/O unit devices and the subroutines which utilize them. 
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Table 2. External Device Usage 



TAPE1 

TAPE2 

TAPES 

TAPE? 

TAPES 

TAPE10 

TAPE14 









MAIN 

X 







INITN 





X 



OUTPN 




X 

X 



GEOMIA 



X 





START 



X 



X 


SURF 






X 

X 

PACH 






X 

X 

PACHXY 






X 

X 

XPATCH 






X 

X 

AEROF 


X 


X 




INTERN 


X 


X 




FORCE 


X 






SPANI 

X 



X 




INTERZ 




X 


X 


CHORDI 

X 



X 




STREAM 

2240-004W 




X 



X 
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10 NCOREL USER MANUAL 


10.1 INPUT PARAMETERS 

Card No. 1 Description Title of Geometry 

Card No. 2 Title Describing Freestream Conditions 

Note : Cards No. 1 and 2 are alphanumeric title cards and must be less than 80 

characters. 


Note : The remainder of the input is in namelist format. Recommended default 

values of certain parameters are suggested. 


Namelist Variable 

EM INF 

ALP 

ISYM 


Description 

Mach Number (M > 1) 

Angle of Attack (degrees) 

If ISYM = 0, half plane symmetry will be 
assumed. 

If ISYM = 1, quarter plane symmetry is 
assumed. 


GAMMA =1.4 


Ratio of Specific Heat 


EST = -0.1 


IC 

JC 


Coefficient of Temporal Artificial Viscosity 
(-5.0 < EST < 0.0) 

Ci rcumferential Mesh Points (IC-1) 

(Maximum IC = 58), must be even number 

Radial Crossflow Mesh Points (OC-1) 

(Maximum JC = 58), must be even number 


Note : IC and JC apply only for crude conical grid subject to conical mesh 

refinement. IC and JC must be even numbers. 


KMAX (1) 
KMAX (2) 
KMAX (3) 

KREF 


KMAXNC 


Maximum number of iterations for each 
conical mesh 


Conical grid refinements. Maximum two 
refinements. Refinement given by 
2(IC-1) and 2(JC-1). Final mesh not to 
exceed IC = JC = 58. (1 < KREF < 3) 

Maximum nonconical iterations. 
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DMIN (1) 
DMIN (2) 
DMIN (3) 


Convergence tolerances for each conical 
grid. Nonconi cal tolerance given by 
DMIN (1). 


Note: Currently, 10"^ to 10"^ is used for the bow shock capturing (IFIT = 0) 

tolerance on the correction to potential. 10“^ to 10"^ is used for bow 
shock fitting (IFIT = 1) tolerance on maximum residual. 


Namelist Variable 

Description 

W(l) = 1.5 
W(2) - " 

W(3) = " 

Relaxation parameter for SLOR 
scheme. 

RELSHK(l) 

RELSHK(2) 

RELSHK(3) 

Conical Bow Shock relaxation 
parameters. Bow shock fitting only. 
0 < RELSHK <1.0 

RELNC =1.5 

Nonconi cal bow shock fit relaxation 
parameter 

IFIT 

IFIT = 0, Bow Shock is captured 
IFIT = 1, Bow Shock is fit. 

lOPTS 

Bow Shock capture options 

lOPTS = 0, Outer boundary not conformed to 

shock shape. 


lOPTS = 1, after first conical mesh bow 
shock is conformed to captured shock shape 
for further conical mesh refinements. 


lOPTS = 2, Nonconi cal meshes are conformed 
to previous captured bow shock shape. 

NS = 4 (lOPTS = 0) 

= 4-8 (lOPTS = 1,2) 

Used in bow shock capturing scheme to 
position mesh beyond initial shock shape 
guess. 

EP > 1.0 (IFIT = 0) 
> 1.0 (IFIT = 1) 

If EP = 1.0, initial conical guess for bow 
shock shape is the rotated mach cone. 

Note: For Bow Shock Fit, IFIT = 

lOPTS = 0, NS = 0. 

1, the following parameters must be set: 

NOUT 

Controls output print, if NOUT = 0, no 
mapped plane variables and derivatives are 
output. Minimum output. If NOUT = 1, 
maximum mapped plane output. 

JOUT 

If JOUT = -1, minimum output for solution 
(only on body surface). 
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Namelist Variable 


IPUN 


DZ 

KRR (1) 
KRR (2) 
KRR (3) 
KRR (4) 
KRR (5) 


RZNEW (1) 
RZNEW (2) 
RZNEW (3) 
RZNEW (4) 
RZNEW (5) 

lOPTG 


If JOUT = 0, body surface and shock 
surface (IFIT = 1) output and three 
coordinate lines are output from body to bow 
shock. 

Description 

If JOUT = 1, maximum output of entire 
crossflow mesh. If JOUT > 1, every JOUT 
coordinate line from body surface to bow 
shock will be output. 

If IPUN = 0, the solution will not be 
punched out to restart. If IPUN = 1, the 
solution will be punched out on TAPES. 

If IRUN = 0, new run 

If IRUN = 1, restart run, reads data from 
TAPES 

Marching Step size 

If KRR = 0, marching step size will 
not be refined. 

If KRR > 0, marching step size will be 
decreased. New step size will be 
DZ = DZ/KRR 

If KRR < 0, marching step size will be 
increased. New step size will be DZ = 

DZ * IKRR]. KRR constrained to -2. 

R stations at which marching step size 
is changed. 


Geometry Options 


If lOPTG = 1, axisymmetric body subroutine 
CONBOD is used to generate geometry (user 
supplied or altered). 


If lOPTG = 2, subroutine CONWIN is used for 
geometry. This subroutine can generate 
conical and 3D wings or wing body with 
circular arc camber and elliptical cross 
section. Requires input on Tape 5. 


If lOPTG - 3, subroutine CONUSE will be used 
to generate geometry. User supplied. Any 
input data should be read from TAPES. 
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NG 

AZ 

BZ = 1.0 

REND 

NT 

IPLOT = 2 
lAERO 


ZCG 

YCG 

CBAR 

ITE 


ZTE 


IZW 


If lOPTG = 4, Craidon patch geometry 
routines are used with Harris wave drag 
input required on Tape 5. 

Number of geometry points generated in cross 
flow plane. Generally, NG > IC after mesh 
refinement. (NG must be odd) 

Maximum half angle in degrees of 
configuration (for lOPTG = 4 only). 

Not used. 

Solution terminates at R = REND. 

Alternate termination of program after NT 
steps. 

Inoperabl e. 

If lAERO = 0, aerodynamic coefficients are 
not computed. 

If lAERO = 1, aerodynamic coefficients are 
computed at end of output. 

Axial moment center. 

Vertical moment center. 

Aerodynamic chord used for nondimensional 
moment coefficient. 

If ITE = 0, trailing edge subroutine XTES 
will not be called. Force calculation will 
terminate at z = ZTE 

If ITE = 1, trailing edge subroutine XTES 
will be called and used by force 
subroutines. 

Default location of trailing edge. If ITE = 
0, a straight trailing edge with zero sweep 
or a body will be assumed where forces are 
computed to z - ZTE. 

If IZW = 0, interpolated z - station data 
used in force calcualtion will not be 
pri nted. 

If IZW = 1, interpolated z - station data 
will be printed. 
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IXWI If IXWI = 0, no chordwise stations are read 

for interpolation output. 

If IXWI > 0, IXWI chordwise interpolation 
stations will be read from TAPEl and x- 
station interpolated output will be printed. 

IZWI If IZWI = 0, no spanwise stations are read 

for interpolation output. 

If IZWI > 0, IZWI spanwise interpolation 
stations will be read from TAPEl and z- 
station interpolated output will be printed. 

Auxiliary Input on TAPEl 

This input uses a format of 10F7.4 and only occurs if either IXWI or IZWI 
are not zero. 

IZWI Input Set: 

Up to 25 arbitrary spanwise z stations can be input for interpolation, 
IXWI Input Set: 

Up to 25 arbitrary chordwise x-stations can be input for interpolation. 
Geometric Input (TAPE 5) 
lOPTG = 2 

A total of six input cards corresponding to five geometric axial sections 
are required for this geometry option with Format 7F10.5, 

Each card has the following format: 

z, 6p, 6j, 6-j-Qj Sp, Sg 

where z refers to z-station and the remaining six numbers are angles in 
degrees. 

Card No. Variables 

1 ZWl, DELFWl, DELEWl, DELTWl, DELCWl, DELPWl 

2 ZW2, DELFW2, DELEW2, DELTW2, DELCW2, DELPW2 

3 ZW3, DELFW3, DELEW3, DELTW3, DELCW3, DELPW3 

4 ZW4, DELFW4, DELEW4, DELTW4, DELCW4, DELPW4 
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5 

6 

Variable 

z 

OELFW 

OELEW 


DELTW 

DELCW 

DELPW 

DELPB 


ZW5, DELFW5, DELEW5, DELTW5, DELCW5, DELPW5, 
DELPB5 

ZW6, DELFW6, DELEW6, DELTW6, DELCW6, DELPW6, 
0ELPB6 

Definition 


z- station 

Spanwise angular location for camber 
If DELFW = 0, entire spanwise section will 
have circular arc camber, 

Spanwise angular location for elliptical 
thickness distribution 

If DELEW = 0, cross section will be 
ell iptical. 

If 0 < DELEW < DELTW, section inboard of 
DELEW has constant thickness, outboard of 
DELEW has elliptical thickness. 

Leading edge half angle (90° - A) 

Maximum circular arc leading edge camber 
angl e 

Centerline half thickness angle of wing 

Centerline half thickness angle of wing-body 
cross section. 


General Rules for Analytical Geometry Package (lOPTG = 2) 

1) The geometry must be initiated with a circular cone at z Z1 with 
DELFWl = DELTWl = DELPWl and DELEWl = DELCWl = 0. 

2) Card No. 2 for z = ZW2 is identical to card No. 1 except DELPW2 = 0, 

3) The geometry will deform from a circular cone at z = ZWl to a conical 
geometry section between z = ZW3 and z = ZW4. 

4) The angular definitions of cards No. 3 and 4 must be identical, 

since a conical section is assumed between z = ZW3 and z = ZW4, 

5) Wing-body cross section can only be generated for z > ZW5. Body 

angular definition only required for card no. 5 and 6. 
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6) Conical wing-body or wing section will be generated between z = ZW5 
and z = ZW6. Angular definitions on cards 5 and 6 must be 
identical. If DELP = DELPB on cards 5 and 6, wing-body cross section 
will not be generated. If DELPB > DELP, a wing body cross section is 
generated. 

7) DELE must be less than DELT. 

8) DELF must be less than or equal to DELT. 


Harris Wave Drag/Craidon Patch Geometry (IQPTG = 4) 

This input is also read from TAPE5. A detailed description of this 
geometry input is located in several references and will not be repeated. 

(see Ref. 10) 

Restrictions: Only wing, fuselage, or wing-body options apply to NCOREL. 

Fuselage must be treated as individual component and should be 
smooth (i.e. , no corners). Wing defnition must extend to 
centerline and origin of geometry must be at (0,0,0). For 
wing-body configurations, apex of wing should be on axis of 
coordinate system (Y^ ^ 

Note: Centerline wing or body twist should be removed from geometry and 

included in angle of attack. Severe centerline camber on bodies or 
wings that cause origin of coordinate or apex to lie outside of cross 
section will not run in NCOREL. 
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10.2 PRINTED OUTPUT 


Note: The first set of 

Output Block 
Spherical Geometry 

Variable 

IG 

R 

X 

y 

z 

OMEG 

PSI 

Output Block 
Map Singularity 

Output Block 

Body Mesh Data 
(Not printed if NOUT 

Block #1 
Variable 
X 

B 

BPR 

BSEC 

BR 

BRR 

BRTH 


printed output echoes the input parameters. 

Description 

NG Input Geometry points are printed at new 
radial station (R = R + DR). 

Description 

IG = 1, lower symmetry plane 
IG = NG, upper symmetry plane 

Spherical Radius 

Spanwise Coordinate 

Vertical Coordinate 

Axial Coordinate 

Spherical angle (degrees) 

Spherical angle (degrees) 

Description 

The location of the mapping singularity and 
derivatives at the current R and new R + A R 
station are printed for both the physical 
and mapped space. 

Description 

Blocks of variables are printed for 
= 0) the body surface (J = 2) used in solution at 
current R station. 

Description 

Mapped space ci rcumferential coordinate 

(X = e) 

Mapped space body radius (B = p) 

Mapped Space body derivatives B 0 

®00 

^R 

^RR 

Br 
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Block #2 


Variabl e 

X 

C 


CPR 

CSEC 

CR 

CRR 

CRTH 

Block #3 

X/Z 

Y/Z 

H 

UI 

VI 

WI 

HI 

H2 

Block #4 


HX 

HY 

HR 

HIX 

HIY 

HIR 

H2X 


Description 

Mapped space circumferential coordinate 
(X = 0) 

Mapped space shock (IFIT = 1) or outer 
boundary (IFIT = 0) radius 

Mapped space shock or outer boundary 
derivatives 

Ce 

^09 

^R 

^RR 

^R 


Physical Space Body Mesh 
Coordinates 

Mapping metric 

Mapped space freest ream velocities 


Radial mesh derivatives pr 
Radial mesh derivatives p0r 


Metric derivatives in computational space 
'’x 


'’lx 

hjy 

hip 

'^^x 
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H2Y 


H2R 

Output B1 ock 

Solution Convergence History 
(Not printed if NOUT = 0) 

Variable 

Iter 

DELMX 

I 

J 

DELAVG 

RESMX 

I 

J 

RESAVG 

KSUP 

DELS 

IS 

Output Block 

Bow Shock Mesh Location 


h2y 

*^2R 


Description 

Iteration 

Maximum correction to potential 

Mesh point where maximum 
Correction occurs 

Average correction 

Maximum residual 

Mesh point where maximum residual 
occurs 

Average residual 

If IFIT = 0, number of supersonic points 
outside captured bow shock 
IFIT = 1, number of supersonic crossflow 
poi nts 

(Only printed if IFIT = 1) 

Maximum residual in shock jump equation 

Shock Mesh point where maximum residual 
occurs. 


Only printed if IFIT = 0. The JS = J mesh location of the last 
supersonic point that occurs starting with J = JC-2 and decreasing J. The JS- 
1 mesh point corresonds to the subsonic side of the captured bow shock. 

Note: The J mesh location must not exceed JC-3 or erroneous results may be 
obtained. 


Output Block 
Mesh and Solution 
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Block #1 


Body coordinates and surface data (J = 2, I = 2 to IC) 


Variable 


RHO 

F 

X/Z 

Y/Z 

Z 

M 

Mr 

MC 

CP 

Block #2 


Description 

Cross flow mapped radius 
Reduced potential 

Physical space spherical coordinates 

Total Mach number 
Radial Mach number 
Cross flow Mach number 
Pressure coefficient 


Bow shock coordinates and surface data (J = JC, I = 2 to IC) 
Printed only if IFIT = 1 

Above output is printed for the bow shock surface. 


Block #3 

RADIAL OUTPUT DATA IN CROSSFLOW PLANE 

This output is printed only if JOUT > 0. Above printout occurs for I = 
constant and J = 2 to JC. If JOUT = 0, Printout occurs atl = 2, I=l+-|^ 
and I = IC. 


Block #4 

INTERPOLATED CROSSFLOW SONIC LINE 

If JOUT > 0, the interpolated physical coordinates (X/Z, Y/Z) of the 
embedded cross flow sonic line are printed. 

INTERPOLATED BOW SH OCK 

Occurs only if IFIT = 0 and JOUT > 0, interpolated physical coordinates 
of the outer sonic line. 

Note: The above printout occurs for each R station computed. 


Ill 


Note: The code will terminate at R = REND or after NT steps. If the message: 

SOLUTION HAS DIVERGED AT R = 

the solution diverges and run will terminate and proceed to next set of 
output. 

Output Block 

END OF RUN ITERATION SURVEY 

This set of output summarizes the number of iterations required for 
convergence at each R Station. 

TOTAL NUMBER OF ITERATIONS = 

This number reflects the total number of iterations required for the 
computation. 

Output B1 ock 

If IZW = 1 and lAERO = 1, spanwise data is interpolated and printed. 
Output Block 

AERODYNAMIC OUTPUT DATA 

IF lAERO = 1, aerodynamic coefficients are computed. 

Block #1 

HALF CROSS SECTIONAL AREA 

This output prints the cross sectional area for each of the interpolated 
z stations 

Block #2 

SECTIONAL OR INCREMENTAL AREAS, FORCES, AND COEFFICIENTS 

This output prints the forces and moments on each spanwise increment for 
half of the symmetric configuration. 

Variable Description 

Z Axial Station 

Incremental : 

SREF(DZ) planform area 
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surface area 


SAREA (DZ) 

FN(DZ) 

FA(DZ) 

CN(DZ) 

CA(DZ) 

CL(DZ) 

CD(DZ) 

CM(DZ) 


normal force 
axial force 

normal force coefficient 
axial force coefficient 
lift coefficient 
drag coefficient 
moment coefficient 


Note: The above coefficients are nondimensional ized by the incremental span 

areas. 


Block #2 


AXIAL INTEGRATION OF AREAS AND COEFFICIENTS 


Variable 

Z 

XMIN 

XMAX 

SREF(Z) 

STOT(Z) 

CN(Z) 

CA(Z) 

CL(Z) 

CD(Z) 


Description 
axial station 

location of centerline or trailing edge of 
wi ng 

Maximum spanwise coordinate or leading edge 
of wing 
planform area 

surface area 

normal force coefficient 
axial force coefficient 
lift coefficient 
drag coefficient 


CM{Z) moment coefficient 

Note: The last set of numbers represent the total areas, forces, and moments 

of the configuration. Coefficients are nondimensional ized by the 
computed planform area. 
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Block #2 

The freestream conditions, total lift, drag and moment coefficients, and 
the aerodynamic constants used in this computation are printed. 
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10.3 SAMPLE CASES 


Figure 44 illustrates a sample user supplied CONUSE subroutine for 
generating the cross sectional coordinates of an arrow wing with 33 degrees 
leading edge sweep and 27.214 degrees trailing edge sweep with a symmetrical 
4% thick NACA four digit airfoil specified for the chordwise thickness. The 
wake extension is specified by extending the trailing edge thickness. 

Figure 45 illustrates two sample input data sets for M«, = 1.70 and a = 5° 
with a marching step size OZ = 1.0. Figure 45a illustrates the input data for 
the bow shock capture method (BSC), and Fig. 45b illustrates the equivalent 
data set for the bow shock fit method (BSF). 

Figure 46 illustrates the printed output when NOUT = 1 and JOUT = 0 for 
the BSF method. The first station printed is at R = 0 or the conical 
station. The NG spherical geometry points are printed at R = R + A R. The 
geometry is then assumed to be conical between R = 0 and R = 1. The 
singularity location in both the physical and mapped space are printed. At 
the conical station (R = 0), all the R derivatives of the body (B), bow shock 
or outer boundary (C), metric (H), mesh derivatives (H^^, H 2 ), and mapped space 
singularity location (S^) and its conjugate (S-j^ = S*^ ) as well as (S 2 ) are 
identically zero. 

Following the mapped space metric data, the iterative solution 
commences. For the BSF method, the solution is converged when both the 
maximum residual (RESMX) of the full potential equation and the maximum 
residual of the isentropic shock jump condition (DELS) are reduced to the 
specified tolerances. The tolerance specified in the input was 10“'^ for the 
maximum residual. An internal tolerance of 10“^ has been set on the 
correction for the shock position C represented by DELCMX. DELS and DELSHK 
represent the maximum residual of the shock jump condition which is used at 
nonconical stations. The location (I,J) of the maximum correction and 
residual for the potential are given and the location of the maximum residual 
in the shock jump condition is also specified. For the conical station, the 
internal flow usually converges first followed by the shock jump condition. 

The body and shock solutions are then printed for the J = 2 and J = JC 
rings. With JOUT = 0, the solution on three rays from J - 2 (body) to J = JC 
(shock) are then printed for 1=2, 9, 16 corresponding to the two symmetry 
planes (I = 2,16) and the leading edge (I = 9). 
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This output repeats itself for each radial station. Spherical geometry 
data is always computed one station ahead of the solution. Hence, the 
spherical geometry is printed at R = 2 and the solution is computed for R = 1. 
Figure 47 shows the differing output corresponding to the BSC method. In the 
BSC method, there is no shock jump condition and the parameter DELS is missing 
from the iteration listing. The tolerance is specified in the input as 10“^ 
and applies to the maximum correction to the potential (DELMX). When DELMX 
converges to less than 10"^, the solution has been obtained. The first 
printout yields the ci rcumferential bow shock location relative to the mesh. 
The body data (J = 2) is then printed. The bow shock location is not printed 
since the shock has been captured. Rays from the body to the shock are also 
printed followed by the interpolated bow shock coordinates. The bow shock 
coordinates are interpolated as the outer cross flow sonic line. Pressure 
data is not interpolated for the bow shock because it is considered to be too 
inexact due to the capture method. 

The next set of printed output occurs at the end of the run and gives a 
summary of the number of iterations required at each step and the total number 
of iterations. Figure 48 shows this set of output for the BSF input data. A 
total of 35 stations was computed. This output is designed to give the user a 
quick way to determine if the run was successful. Since the maximum number of 
iterations specified in the input was 200, all the stations have converged to 
the specified tolerance. 

The user has the option of specifying arbitrary spanwise and chordwise 
stations where the pressure coefficient data will be interpolated. Figure 49 
shows an example of this output. 

The last set of printed output is shown in Fig 50 and corresponds to the 
aerodynamic output data for the BSF run. 
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100 

no 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

2246- 


= SUBROUTINE CQNU^<2 j XX » IFLAG? Y j NX » NY ? NZ ) 

= COMMON /TIP/XTIPjYTIP^XSINGjYSING 

= C 

=C ARROW WING WITH SPECIFIED AIRFOIL SECTION 
=C 

= PI=3 *14159265 

= X=XX 

= THLE=33* 

= THTE=62,786 

= TANTHL=TAN(THLE*PI/180* ) 

= TANTHT=TAN(THTE*PI/180. ) 

= XT=Z‘*TANTHL 

= CR=^20* 

= ZA==CR*TANTHT/ ( TANTHT-TANTHL ) 

= XA = ZA!)(TANTHL 

^ IF CX*GT*XT) X=XT 

= IF (IFLAG,EQ*0) GO TO 500 

= ZLE=X/TANTHL 

= ZTE=CR+X/TANTHT 

= C=ZTE-ZLE 

= ZC==(Z“ZLE)/C 

= CALL AIRFOL(ZCjYTHC) 

= YTH=C*YTHC 

= IF (IFLAG*EQ*1) Y=YTH 

= IF (IFLAG*EQ*2) Y=-YTH 

= RETURN 

= 500 YTIP=0* 

= XTIP=Z*TANTHL 

ZC=Z/CR 

= CALL AIRFOL (2C,YTHC) 

= YTH=CRiJ:YTHC 

= XSING=SQRT(XTIP**2-YTH**2) 

= YSING=0* 

= RETURN 

= END 

049W 


CONOOOlO 

CON00020 

CON00030 

C0N00040 

C0N00050 

C0N00060 

C0N00070 

CQN00080 

C0N00090 

CONOOlOO 

CONOOliO 

CON00120 

C0N00130 

CON00140 

C0N00150 

C0N00160 

C0N00170 

CONOOISO 

C0N00190 

C0N00200 

C0N00210 

C0N00220 

C0M00230 

C0N00240 

C0N00250 

CQN00260 

CON00270 

CON00280 

C0N00290 

C0N00300 

C0N00310 

C0N00320 

CQN00330 

C0N00340 

C0N00350 


Figure 44 Sample of CONUSE Subroutines 
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450 = 

460 = 
470= 
480= 

490 = 
500= 

510 = 

520 = 

530 = 
540= 

550 = 
560= 

570 = 
530= 
590= 
600= 
610= 

620 = 
630= 
640= 

650 = 

2246-049W 


SUBROUTINE XTES ( X » ZTE » DUHHY ) 
COMMON/TE/IZWrITE?XTEC»YTEBKjXTIP 
IF (ITEtNEvO) YTEBK=0. 

PI=3*14159265 

CR=20* 

IHTE=62.786 

THLE=33*00 

TANTHL=TAN(THLE*:PI/180* ) 

TANTHT=TAN(THTE*PI/ 180.0) 

IF (ITE.NE.O) XTIP=CR*TANTHT/(TANTHT-TANTHL) 

ZTE=CR+X/TANTHT 

RETURN 

END 

SUBROUTINE AIRFQL (XCtYTH) 

T=.04 ■ 

IF (XC.GT.l.) X=1.0 
IF (XC.LE.l.) X=XC 

YTH=(T/.20)2t( .2969*SQRT(X)-.126-*X-.3516*X>|{X 
l+.28435tJX*X*X-.1015*X*X*X*X) 

RETURN 

END 


Figure 44 Sample of CONUSE Subroutines 


C0N00360 

C0N00370 

C0N00380 

CON00390 

CON00400 

C0N00410 

C0N00420 

C0N00430 

C0N00440 

C0N00450 

C0N00460 

C0N00470 

CON00480 

C0N00490 

CON00500 

CON00510 

C0N00520 

CON00530 

CON00540 

C0N00550 

C0N00560 
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100= 4-III6IT ftIR0IL.» ♦ARROW WING 

110= H=1 ♦ 70 » ALP=5 ♦ 00 » THC=30 ♦ 00 » DELC= 1 ♦ 000 

120= $INPUT IC=16 7 JC=16 7 KMAX ( 1 ) =200 7 KMAX ( 2 ) =200 ? KREF=1 7 IPL0T=2 7 KRR ( 1 ) =0 7 KRR (2) =0 7 

130= KRR(3)=07KRR(4)=07KRR(5)=07nZ=l*7RZNEW<l)=1004 7RZNEW(2)=100f 7RZNEW(3)=900* 7 

140= RZNEW { 4 ) =1000 ♦ 7 RZNEW ( 5 ) =1 000 ♦ 7 EMINF= 1 . 70 7 ALP=5 ♦ 00 7 AZ=33 ♦ 00 7 BZ= 1 ♦ 000 7 

150= GAiiMA=l,47EP=l»407lQPT6=37l0PTS=07lFIT=07lAER0=l7lZW=07lTE=l7lSYH=07 

160= IXWI=107lZWI=57ZTE=20.007CBAR=14*7477YCG=0*7ZC6=16.70l7 

170= RELSHK < 1 ) = ♦ 50 7 RELSHK ( 2 ) = * 25 7 RELSHK ( 3 ) = ♦ 10 ? RELNC=1 . 5 7 KfiAXNC=200 . 

1S0= W(l)=l*57W(2)=l*57DMIN(l)=i*E“67DHIN(2)=l*E-67EST=-l*E-l7NS=47NT=507W(3)=l»5 
190= MQUT=l7 J0UT=07NG=397lRUN=07lPUM=07REND=34*?D«IM(3)=l*E-67KHAX(3)=200 $EM)t 
200=5>0 10»0 15*0 20*0 25,0 

210=0*0 2*0 4,0 6*0 8*0 10.0 12*0 14*0 16*0 18*0 CHORD! 


100= 4-DIGIT AIROIL* * .ARROW WING 

110= N=1 ♦ 70 7 ALP=5 * 00 7 THC=30 * 00 7 DELC=1 * 000 

120= $INPUT IC=167JC=167KMAX(l)=2007KHAX(2)=2007KREF=l7lPL0T=2?KRR(l)=07KRR!2)=07 
130= KRR ( 3 ) =0 7 KRR ( 4 ) =0 7 KRR ( 5 ) =0 7 DZ=1 * 7 RZNEW ( 1 ) =1 00 , ? RZNEW ( 2 ) =1 00 * 7 RZNEW ( 3 ) =900 . 7 
140= RZNEW!4)=1000*7RZNEW(5)=1000*rEMINF=i*707ALP=5.007AZ=33*007BZ=l*0007 
150= GAMfiA=l *4 7EP=1 *07 I0PTG=37 I0PTS=07 IFIT=1 7 IAER0=1 7 IZW=07 ITE=1; ISYN=0? 

160= IXWI=107lZWI=57ZTE=20*007CBAR=14*7477YCG=0.7ZCG=16*70l7 

170= RELSHK! 1)=*507RELSHK(2)=*257RELSHK(3)=*107RELNC=1.57KMAXNC=2007 

130= W(l)=l*57W(2)=l*57DNIN(l)=l*E-37DHIN(2)=l*E-37EST=-l*E-l7NS=07NT=507W(3)=l,5 

1 90= N0UT=1 7 J0UT=0 7 NG=39 7 IRUN=0 7 1 PUN=0 7 REND=34 * 7 DMIN ! 3 ) =1 * E-3 7 KMAX ( 3 ) =200 $END 

200=5*0 10*0 15*0 20*0 25*0 

210=0*0 2*0 4*0 6*0 8*0 10*0 12*0 14*0 16*0 18*0 CHORDI 

f I- 

COMMAND- 

2246-048W 

Figure 45 Sample Input Data Sets 
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NON-CONICAL RELAXATION SOLUTION 


r-j 

o 


li 13 Cr 
M X O 
<x »“• a: 
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RZNEW(1)= 100.00000 RZNEW(2)= 100.00000 RZNEU(3)= 900.00000 RZNEU(4)=1000. 00000 RZNEW(5)=1000. 00000 



LEADING EDGE DR MAXIMUM HALF WIDTH COORDINATES 


I 


r-4 r-4 c-j CM c-4 r-4 cm r>4 r*4 cm cm r^i r-4 cm r-4 cm r-4 r-j c-4 r-4 r-j cm cm C“-4 C4 rg r>i r>j C4 c>i r-j C4 c-4 C4 C4 C4 r«j 
ooooooooooooooooooooooooooooooooooooooo 


LiJ LU JjJ LjJ UJ Ld UJ 

-.-H rs r4 o -o c^j 

^ b"J -O CO ^ 

♦— I UT ro 'O ^ CO ^ 

cnrororO'^rb'^MDCO 


LU Ld 
r J r-4 
|j*:i 0^4 CD CO 


UJUJLiJLiJLLJLLlUJUJUJLJLlJUJUJLLJLiJLiJLUUJLLJUJLdlii 

CO ii“D ■«r rs ^ o o c> o o M3 •c' LTi CO <r CN t-i r-j r-i 

r-4 O r-4 ''O M3 O MD -.-1 NO fJ O C'4 ^ 00 CD CJ bT 

ro On o r-4 ro r-4 o L"> CD o CO m o r-j ro r-j 


o-f-irO’^r'NorNvONO 


r-4 r-4 r-4 ro r-4 CM r-4 -r-H o On rv M3 ■«- r3 


CM r-4 r-4 rNi r-4 r-4 r-4 ro ro ro ro ro ro ro ro ro ro CO C4 r-4 r-4 r*4 r-4 r-4 r-4 


00 -rH ^-1 

^ 00 LT -«r 
O 00 NO LT 


M3 LO ■•3“ 
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Os 


tn 

m 

IK 

M3 

O' 

TK 


eg 


M5 

n 

■«r 

o 

O 

ro 

CK 

ro 

tn 

O' 

IK 

o 

rg 

Oh 

MD 

ro 

CK 

o 

rg 

CO 

04 

FK 

FK 

MD 

ro 

tn 

MD 

7—1 


CK 

3: 


CO 

tn 

Os 

Os 

tn 



rg 


04 

IK 

o 

rg 

ro 



ro 

Os 

rg 

CO 

Oi 

tn 

o 

o 

MD 

Os 

IK 

«T 

ro 

ro 

■«r 

tM 

MD 

tH 

IK 

t'i 

7-1 

CK 






O 

rg 

tH 

rg 

o 

ro 

CO 

r-4 

in 

o 


tn 

M3 

o 

CM 

ro 

IK 

M3 


o 

Mi 

tn 

rg 


O' 

IK 

IT) 


ro 

rg 

7-1 

o 

un 

FK 

O 

ro 

FK 

04 

MD 





UJ 

CD 

I 

rg 

1 

04 

1 

7—1 

CO 

1 

IK 

[ 

M- 

1 

M3 

t 

ro 

MD 

1 


tn 

1 


in 


1 

ro 

ro 

1 

ro 

ro 

1 

rg 

CM 

1 

CM 

eg 

I 

7-1 

1 

7—1 

tH 


t 

7-M 

1 

CK 

CO 

CO 

FK 

M3 

M3 

tn 



1 


Lc: 


CM 

to 


un 

M3 

IK 

00 

CK 

O 

7— f 

r-4 

ro 

‘•T 

tn 

o 

IK 

GO 

CK 

o 


04 

ro 

<T 

in 

MD 

IK 

00 

CK 

o 

tH 

CM 

ro 

•«r 

in 

MD 

rs 

CO 

CK 

O 

7-1 

CM 

ro 

Ul 

H- 

1— i 












7-H 

7-^ 



7—1 




04 

rg 

04 

r-4 

04 

04 

04 

rg 

eg 

rg 

ro 

ro 

ro 

ro 

ro 

r*3 

ro 

ro 

ro 

ro 






124 
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3721E-03 2 10 .3939E-04 .6296E-01 2 2 .B090E-02 0 .llE-02 


C>4 rj CN CN r-4 CM C4 CX CN C4 C4 C4 C4 C'i C4 C4 rj C4 C4 04 04 C4 C4 04 04 04 04 C-4 04 CM 04 04 C-4 04 04 C4 C4 C4 04 04 C4 C4 04 04 04 04 04 04 04 

o4roi^mroMroroh'>rorotororororororoMrohorofOhorohororo^^^^^^-^^^<r^^^^'^^«o^'^^^^ 
oooooooooooooooooooooooooooooooooooooooooooooooooo 

I I i I I ! I I J J I I I I I I I I 1 I I I 4 t I < I I I t I t t I I I I I I j j I I t r I 1 I I I 

U4 IxJ LaJ LlI lij LaJ LlJ LlJ LlJ LlJ LU LjJ LlJ UJ LlJ LtJ LiJ LU LlJ tiJ LU IxJ LiJ tiJ t.iJ LU UJ LU UJ tiJ tiJ t.tJ lii lil lil i.i i.i ||^ ||> iii i.i ||} i.i ni i«i 

ojj‘iooo^oorooxMooin»-*co^ 04 0'-rs'^ro*Ho^rv.-<im‘^04^-iO'Ooo-i-«U“5co?oo3fot>'in^ooirt04Ch>-OsiofO‘i-iOHoa*»o 
■T-io>oocDrs‘OHomm^^'^rororoo4 04 04 04 04 T-i*-i^^ ‘»h »-4^«-<o^oooorN.^'Oinm-^^^rofotoo4C4 04 04 04 TH^T-f 


oooooooooooooooooooooooooooooooooooooooooooooooooo 


04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 

ooooooooooooooooo 

J I I I I ( f I I I J I I I 1 I I 

LlJ III i|i lit ffi 1^1 III III ^|| I,i ^.1 111 111 111 111 111 ifi 

ooo.or^CTN^O'^o^rxCN'r-i^toO'T-* 
04 04 b") b"> rv ^ L*5 CD bT <* Si O b*:i li“> o> o 
‘Q*cnc4r>.04C0‘^orv.<rTsCKSD'«T04ocs 
Os nq ^ ifj iTj ^ ro ro ro 04 04 04 04 04 ts 


04 040404 04 04 04 roroMrooiroroporo 
oooooooooooooooo 

I i I I ( r ( I 1 I j 1 I 4 I I 
LiJUJUJU4UJUJU4UJUJUJUJLjJLUUJUJUJ 
-orof0 04 04*»HaD'«rD'ONrso«.oov<iro 
<*orsin«T^^04ro*»-ib’i^CNOsor*s 
0HSD^r004-r-iOSDCDT-^^CD04 0sr0C0 
■rM^,-l^i-H^»HOvGOCOOsS5**OU’5b“D^ 


rorororororomrororororot^howroro 

ooooooooooooooooo 

I I 1 I I I ( i I i 1 I I J t I I 

iii iii til 111 III III III 111 |«i iii III \ti i»i III j.i 111 III 

so-»-«s3CDUTtnrNCOs>04r4sDrO'.s-^ocD 
0HT-H{S.sD00 04 C0'<»OC0-rHm*.-i00-«0b*i-<r 
'^-i-iOs-^rrHCNsi’^r 04 ocsOs<i<rho 04 -rH 
•«r ro ro ho 04 04 04 04 04 -rs th *.M -.-i th -rH 


04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 


04 04 04 04 04 0 4 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 0*4 04 04 


•H tH %-C -.-I ---I -r-f ^ ^ ^ ^ T-l ^ -r-< ^ ^ ^ OJ 04 04 04 

ooooooooooooooooooooooooo 

f J I I I I I I I I I I I I I I I I t t I I I I I 

LU LU L»J LlI III tjJ lil III til 111 til l»l tjJ lil LlJ tiJ lil III III ill III til 111 fil |il 

r4«5--i-H^rOCKrOPsCNCDOi>'O^Ob">r^^COT-i'«T^LTO-.bT 
COOOIs.b“)siOJO.CDr^^bT‘4bO‘^0*Ob‘irO'«r^OJi-HTHKi 
^rocD^c>Ov'«ri-icci*-o«3-o4 0 0Nrv.'*o-^ror4'f-*0'^coi-<^ 

UT bT ■«T ^ ■«r ro m ro 04 04 04 04 04 -.-H ^ tH -r-H -rS -rH tH -r-H Ch> CD CO fv. 


I t 1 I I kf I I 1 1 


04 04 04 04 CS 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 


04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 04 0 4 0 4 0 4 04 04 04 


04 04 CN 04 04 04 04 04 04 04 04 04 04 0-4 04 04 04 04 04 04 04 04 04 hi ro 
OOOOOOOOOOOOOOOOOOOOOOOOO 
I r I I 1 j I I I t I I I I I j j I t I I I I I I 

UJ LlJ LxJ UJ UJ UJ UJ LlI LlI LlI LlI UJ LlI UJ LlI UJ UJ UJ UJ UJ UJ UJ UJ LU UJ 

S3 ro -«r b~> o sD CO Oh Si o ro 04 04 CO CO ^ o o CO CO Si m bT uT 

’«r0H0-0-0xSDOSJ^r^ 04 03-^**00.OU^-*-*C0b'3<4"^b0*OC0 
03 04 0.04CD^t-40v^^Ch*‘^’«r 04 00NO».SD<rr<JOJi-iO-<ICCi 
'O S3 bT b"D ^ ^ ^ ro ro ro 04 04 04 0 4 04 -.-H T-I ^ -t-i ^ 7-1 T-I -»-i Ch CO 


I - I I I r I 1 I ] I I I 


OOOOOOOOOO 
I I I I I I I I I I 
LUUJUJlJJUJUJLUUJUJUJ 
ISO sDCNcsor^ococN 
tS 04 ^ CN S3 S3 S3 CK 04 07 
''Ohioo^LThi'SCNCDo 
hi ro m 04 04 04 04 7-1 ^ ts 


«r ^ ^ ^ ^ ^ b“i bO b*i tn 
OOOOOOOOOO 
f I I I I I J I I I 
UJUJLULULULULUUJLULlJ 
04^7-*<rO.Ov’«5-COS3Cr' 
^»HOCNCNOinCr**0'0 
UT'«rhi7soooj<rco>i-i 

7-H 7—1 ^ 7-H 7—1 CN 00 Os 


in m m m bT b*i b*i U”i u*i m 
OOOOOOOOOO 
r I J I I I I J I I 

LU tu lU 111 UJ ill III 111 l»i III 

b^ Os b*i T-i S 3 hi ro o ■«T 04 

co^b“3ocoob'ihihis3 

L’iOb'iTssihiCNsDhiO 

s3S5b’ib'i-«r^hihihihi 


m b*i b") b") b“i b *5 uT in in \n 
OOOOOOOOOO 
I I i I I I I I I I 
UJUJLlJUJUJUJUJIiJUJLiJ 
hihimos7-HCOOD7-i<rco 

ts CO Os rs o hi 00 b*i 04 o 
coiDhiT-iocasib’i^Tro 
04 0-4 04 CN 04 tS tS 7-1 IS 


inb’iLnsDsisi*^s3>os3 
OOOOOOOOOO 
I I I I I I 1 I I I 
UJUJUJUiUiUJLUUJUJLU 
•r-i oi hi CD CN 04 0 4 hi hi rN 
OOTSO^U'i7-i04 03CD 
04 T-< O hi m 00 04 •<» O bO 
^tSt-iOsCOOsOsnOsOLiT 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

^7-1 7— trH^THTHT-HTSTS-iST-HT— I t-Ht-HtHt— It-HtSt— I l<rM7 — It— t7"H7 — It— (7-4t— I t— It-It-It— It-It-I 

r-4 04 0-4 04 0-4 04 04 C'J 04 04 04 04.04 04 04 0-4 0-4 04 04 04 04 04 04 04 04 0-4 04 04 04 04 CN 04 04 0 4 04 04 04 04 04 0 4 0 4 04 04 04 04 04 04 04 04 04 


hihihihihihihihihihihihihi'«i''^<r*^^^^^^^^^^^^^^^«<T<r^^‘<r’^^'^^bTb‘ib'ibTbTbTbTb'i 

oooooooooooooooooooooooooooooooooooooooooooooooo 


UJ UJ Ui LU LU LU UJ LU UJ LU LU LU UJ LU LU LU LU LU LU LU LU LU LU LU LJ LU UJ LU LU LU LlI LcJ LU LU LU LU UJ LU LaJ LU LU LU LU UJ LU LU LU LU 
O'OCN»-^C^*04C0^^00CNOCN0shs0ss0O0''0C0O4•^O0*•hs 04 04 U‘i0v04hiO 04 0□'*0'0^*s0s’70hi047-*hi<"CNhi04hi 
C0<j"04r0'<rc0 04 C0inh0hi04hi04^hi000sr4 7-iin04his0hih3LncNb'i^<T<iCNhi0s'C^hihi^00oc0 04 7-4s0bTC0 
oos37fojooohs.m‘<rhi 04 7-ioir>OsOhihs04hsC4co^ors.^Tsoo>o^ 04 0cohsin^hi04 7-HOb*icoo^ooo4rs,r4 


LU LU LU LU LU 


LU LU LU LU LU 


■ if III I »l It I III 1 I.I 


oo'O'^cNooohstn'^hioj 


omos0hihs04hs04co 


04O4 04 04 O4-rSTH7ST-iT-i7S7-H7SCNC000OssOsObrib'i^^^hihihiO4 04 CN04O4T-<T-iT-iTHT-i-T-i7-iTsCNC0COrssOsOmb"5 


■^lO'^rNOOosOT-iojhi 


■^b*i'OrscoCNo-.-io4hi-^ 

(OLnb^Lnb'il/isQsosiTOso 


tn-^hsCDOsOTHCNhi 

sO'«o*^**0'^hshshsOs 


-^LnsJ07a30sOTH04 

OsOHps.rshNrs.ooco(DOOCo 


m>0hsC0CNO7-i04hi^ 

CDOOCDCDOOCNCNO^OsOs 
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BODY COORDINATES AND SURFACE DATA 


o o o o 

o o o c> 


o o o o o o c> 
o o c> o o o <? 


^cM-«rvnc>0'*r'Oiii'0^2;?:52D; 

CDCao-.'»-i«rO'irtOOOcor^|^<^(^ 

rjroro»r'«r«rLnmMrJi--'’-‘'^*^'^ 


■«r-rHoooo— »oooo 
ooooooooooo 

UJUJUILJUJUJLULUULU^ 

ooor-^co<icr-»rir2«^o 

cKt*3«c>5ra;Menjoo-c^ 

r-.'Ot»-r^ooifJ^oogw 

^OO.^CMWtNOr-J'OSS'O 


UJ Ui 
M2 

o o 
if> ro 


UI UI 

in CD 
O' in 


o o o o o o 


OOOOOOOOO 


UlWUlUlUlUIUIUJUjyUJUJUIJiJ}^ 
§§■<*■ 111 

oo-fv<---<oo'orMroO'«oinoo2; 


Soooooooooooooo 

lit ^.1 1.1 111 111 111 LU UI UI UI Ul UJ Lil ^ 

oO'Cor'i)">yM'COT'''Oco2;S:l^S!r 

m c-4 c-4 r -1 Cl c''i c*'< ^ 2 ^ "i ^ 


Oooooooo 
oooooooo 

LUUJLUUiUJUJLUUJijjLUUiLLjWyy 
ococ30fvcKcrir^oijT^ro^2"Ogg 
so oo Ch- ^ CNi ro ^ CO ^ o ^ o 

CM oj rw r*j ro ^o hO Ki 04 CJ r-i ^ ‘r<H ^ 


OOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOO 

ujSluujujluluujujuji^^ujwi^ 
^ *f+ b*i ^ ^ rvr ^ M uo o c-4 ^ ^ 

rvic-i'M'MOO'O'Orvjmt^o-O'M^ 

oacccocDCDrMrM.cDcoaacocDO'O'O' 


OOOOOOOOOOOOOOO 

^ -f- *4- -+- ~C '+' *+' 

UIUJUJUJUJUJUJUIUIUJUIUJUJIUW 


OOOOOOOOOOOOOOO 

iLClJJUJUAUJUiUiUJUJUJ^^UW 

C*4 >0 b*^ ^ TS, t<5 >C Cr-' ^ C'J ly^ b^ 

r<i rtv CK in CO '*0 -43 CO r*-" o*» Cr- *-o 

«3-Kot*5mr*3mf*5Mm-«r<"ii^minm 


UI 

o 


o o o o o 


o o o 


oooooooooo 


*-( 


'4' -+- -t- -+■ 


OOOOOimOO 


I I I I 


UiUJUtUJUJUJUlUlUlUlUJLUUJ 

^ T-1 vO -C 111 O' CM 'S- T 4 5^ lO ^ 

CO UI CO r-i r-J *r CM n 2 JM. ^ CM CN CO 
ot'viMUTrDOacr'C'JO'OOfnincM 

£M r-J '-H O' 


UI 


un f*i 'O O' 


lift 


« I I 


«9-OOOOOOOOOOOOOM 

i.mOOOOOOOOOOOOOi-1 

LiJujtijuiUfUJU/UJUJUIUJI^UiWy 

fnrM'OU’^'OOfO’-iMO'Oin^r^'O 

<r T'' m r-J CD r-4 o o n CD « in CM ^ 

r-.j I'' ^ O' p' O' p'l O' ^ o ' * rn m 

CJ ^ C'4 <■ b“l 'O ^ **0 li^ ^ ^ ^ 

I 


C=i 

1Z. 

<t 


OOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOO 

o o o o O O O o o o o o o o 


•4- -+' *4“ . ^ 


UJUIUIUlUIUJUJUJUJUiyWyi^W 

^ CM O'! o C4 ro UI fv •<r -«r u^ I-' CP C;4 

fM. rn r-J fM. pM rsi ^ CM 

CP'OCDthUTOMOO 
r-4orofnt>'<‘rM.'pco 


CD 


CM PO »r CD 
. u"> oJ -o UI O' 

QQfM^VpO'-in'-'MmjM'CncpO' 


I I I I 


I I I 


^ o O O O ^ ''M ^ 

tMOOOOOOOO 


O O C> O PO 

0 0 0 0''* 


UlUJUJUJUJUJUlWWyyWt^W 

cp'ircpo-Hror-i'Ci'O'O’TrjinrM 

PM n <■ D4 uT o un (M, M 


o PD o 


M T'4 PO CO 

c> C“ P'' 'O 

-fM '"I O' 02 


CO 

•X. 

•— I 
DC 


O f— ( T— # W 


ooooooooo 


1 I I t I I 


I I I 


UJ UI U4 UI UJ UJ Uf 
(jT O CP Pp U"> 'O 

-1-1 o 

C& VH 

CD {*i 
"I CP CD PP pp 


CO ^ Pm 


i.i 111 Lij UI Ui UI UI 

to CM MJ 00 «r 'O 

O' m CD M Pm I-I 

(S' '•T rj tm r J 'C CP 

O O' CP O v"i ''I O' 

PM 'O nO Pp pp PP 'D 


I I I • 1 I 


I I 


00 CM 

o o 

c-4 CM 


PP 00 CP Ii1 


CP Pp 'O pp 
CP CP CP CP 


po 

c> o 

CJ CM 


c> o o o o o o 
o O O Oi o o o 

U1 Ui UI UJ u> Ul UI 
in pp PD to o ■P4 tp 
Cp 00 fp- oo ^ 


o Pp M3 pp 


c-4 00 
O' o o 
•rM CM C-J 


r-4 ro <r n •<) Cp CO CP o CM po < r i n 


OOOOOOOOOOOOOOO 


ricriOOOOOOOOOOOOO 

OOOOOOOOOOOOOOO 

UlUIUIUlUiUlUIUIUJUIUJUJWUJUJ 
cpor-4-^c'-4PonKi»MCO^ro^rp^ 
^moPDCa-M-'J-CP'OMO o ^ ^ 

pp.tnoc'4PO'OCMt*i»M<'0^oopfn 

^^iTPCpG3Pp.rpfP030'-i-4C-4^^in 


c-jPO'^m'OppcDD'O 


cj po •»• m 'O 
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Fiaure 46 Sample Printed Output for BSF Method (Sheet 7 of 14) 



RADIAL OUTPUT DATA IN CROSSFLOW PLANE 


LULUUJLiJLUUJLiJUJUJLxiUJUIUJLiJLLJ 

CO MD o r-j MD CO o o- o r-j 00 O' r*:i -.-I 

r-j ^ ro ro rj cj cn ^ rj Ch. Ij"D o 

•oiJ*:ir*jairorxor’Jvu">'«rrooLncN 

OD CO CO rs . r-v Mi Hi bri -«r m c * j o CO -o 

hi ro ro tn Ki ro hi ro K) ho ro ro r*4 rj 


UJUjLlJLiJLLJUJLUUJLiJLiJLULLlLiiUJLU 
in rv hN U'i Oh Ch- CO hi ^ GO C4 ^ o 
'<iOT-ioor’4Hors.fv.fMr-Hij'io^r<50 
OHHQ^>oojHOO'^coc-4Hoor*ib'irM 
o-i-HOOHCOHOb^ro-rHocohNinhi-i-i 
b“i b“i bT ■•r «r ^ ^ TO hi hi hi hi 


ill til 111 111 III III Ifl lij III ttJ Irl l.tJ IjJ III 111 

o UT T-H -»-H o> T-« o o *»—i ch o hv ^ 

04 Oh CD ho r-4 h'j hi hn, hn. >o |v. hi th u“i h-. 

OHCObTTH'«r<rr-jcahir-.ohiHOOHhi 
r^b*i<T'iMrH^hiOH<romi-»*HO'rH%oc4 
•<1 hv -rH c-4 C'j hi hi ^ UT uri •<) <1 hs hs CO 


UJLULxJLlJLiJLLiUJLiJLiJLUUJLiJLlJ 

hir^^hHCO^m»-*orOL’ihira 

HOii‘ir-4COhiorHshicO‘^^-r-ioo 

co‘^r-j«r»HCDhicO'i-<^<»oooH 

•HOii*iCDrJhH^HS30b‘iOHhihs-r-4 

Csi r 4 04 hi hi ^ «r irt bT iiT ho ho 


ooooooooooooooo 

-+“ “f“ -h -H -+" -f- -+** “+" -4“ “f“ ■+" ■+“ *+■ -+- 
LULULiJLlJUJUJUJUiUJLULUUiUJllJLlJ 
orM«^^ocor4r'4 0H^OHhiCNCoch 

o CO u*i CK r«4 c-4 c'4 o *hq r-4 -hq o r>i ^ in 

oOHOHcncohH<»b’ihir*4 00HrMb"5hi 
hi r-4 rj r-i r*4 o4 r-j r j 04 oj r-4 -f-i th -ci -th 


ooooooooooooooo 

"+• "+H ^ -4- -+“ -h “h ^ ^ *+■ -4“ •+“ “+“ 
UJUJLiJUJLiJLiJLaJUiUJLLiUJLLlUJLLlliJ 
hi b"i ^ th o o th r-4 ^ 00 hi 1-1 -rH iiT r«4 
«f-»CN00hH.b‘iCMC0h>hH.ohim'*0'H0'0 
rvHC'O'OHOHOinb'i^^hiCM^^ODH 


Ijl Ifl 111 li) ^II l.i III tii Ltt III lt> itl IjJ Itl 

ooHhihsoocohiU ' iHOHorMhirsC '4 
o o hi NO ■.-! hO -i-H CO b“i hi 04 04 ’•r hi 

O O O O -r-H t-h 04 04 hi «r b“i -<1 Os CO O 
hihihihihihihihihihihihihihi^ 


ooooooooooooooo 

LlI LU LlJ UJ LU L^J LU LiJ tj_l 111 liJ t.l lil l.i lit 

o,co»-ibTcooNb*ihiC'4'«ro3'0-i-iU‘i*^ 
-rHHOo4'r-io*4^hsOhiHocNhicohiaH 
ocKOTHO-4hi^HOhxcociN^c-4*^in 
04 ■»-» 0*4 04 04 04 04 0-4 04 0-4 0 4 hi hi hi hi 


ooooooooooooooo 


ooooooooooooooo 


00000000 0*0000 
ooooooooooooo 

“4“ -f- ‘I- H- -h -4- H“ “fn *4- -4“ "4- -4~ 

iit III III 111 III ill III III Lil t«J III LiJ tiJ 

^ ^ r-H. C ‘4 ^ LT T -4 b"J hi hO o rv. o 

CO hO Mi T 4 in hN 1 -H pN . CD - 3 “ CD Oh f-H 

o r -4 hi ^ U“i U"i Mi oi o hi CD 

CD C' 4 Mi o ^ CO r -4 Mi o b“i Oh hi 

r-4 r-4 hi hi -«a“ U"i L"i Mi ho Mi Pm 




^ ''T tt -«r <r 


LlJ UJ LU LlJ LlI LU LlJ LU LiJ LU LU I i 1 1 1 1 I 1 1 lil 

■«r HQ CD b"i o bT o UT b*i n bT m ^ pN. 

HQMif-ihihiH-»oaHO^*«-ic4coOHHO 

r-4 O CO bT r 4 Oh Mi r-4 O hs bT PO -»-H o o 

rN.or4b"iCOOhiHocK-r-i^TrHvOhiMD 

t-h C‘ 4 r-4 c* 4 c*-4 hi hi hi hi ^ ^ b*i in LT 


«j- *«r bn b"i bi b"i b*i in U“i ^ ^ 


lit ill III 111 lil 111 » 1 1 1 1 j ii \ III 111 

hi -r-t -T CO CO o C'4 hi CO ho r4 

<rhi*p-HMirN,coopN.-^hiMD 
h“4 ch b*i hi r-4 r-4 -r-n r-4 
H-I ^ LO C4 C|H Mi hi o Pm C3 -i-H 

r-4 r 4 b“5 hq vO PM CO Cth Oh 1-1 i-H 


LU LlJ LiJ LjJ LjJ LlJ LiJ LU LiJ LtJ Lii lit Irl iii itl 

T-I Pm rv hO hi CD b"i C4 PM pN. Li r-4 -•-H hC r-4 
•^•so-'Or-4hirNvinHoc»Hn'^HOT-icKfi 
Ch Ps pH, Oh C4 MD r*4 Ch pv. pH^ CD o <• CO hi 
^Micoohib'iCDOhi'QOHhiHoaHO 

HONOHOpH - PN . rHvpvOOCOCDCDOHOHOHH -1 


LULULUUJLlJUJUiUJUJUJIJ-lUJ 
04 pN. O U“3 On PM 11 i-» L"i in C'4 hi 

•«r Mi r -4 hi bT b“5 o On u"i -«r hi 

■fh b*i ho r-4 hi Oh o hi no hi b*i 

<J" hi hi b”i hO On T-1 -«g- pH^ ii 

Oh CO pN. Mi bT "T ro hi 04 *1 th b"i 


UJLUUJLjJUJLiJUJUJUiUJUJLiJ 
LI pn. Oh pH^ Qh. <T b*i hi o r-4 ^ r-4 
L1COO-T-*OCO^ChhiHO-^HO 

hi r-J ^ rv r*4 CO pH ph^ o -«r Ch. hi 

0‘«rair-4pN.iiHOTipHwr-4 0 0H 

hq b"i -«r -*er hi ro r -4 r -4 i -< ii 00 hi 


LiJ LiJ LlJ LU UJ UJ’ L&JLiJLiJLiJLiJLiJUJUJUJ 
On b“i O NO 1P4 CO hi Ch b’i O S3 r-4 PN. hi Oh 
o>i-ihi^HOpN.OHor-4*«rbTrNCDOT-i 
cohirHv-nb*Di>hiCDr4HOO'^CDhirM. 
orghiL'iHOPMOHorghiii’iHor'HONO 
r-4 rg 04 rg r-4 r-4 r 4 hi hi hi hi hi hi hi -«r 


LUUIUJUJLLlljJUJUJUJLiJLiJLiJLiJUJUJ 
rgONMi-^-*HcoLnrgopN.-^-riOH-»ohi 
b*iNocoorghimrN.OHor'i'5“mhH.OH 
hiU‘iPH^or4^NocoohibTrH^Ch-^hi 
o-nr-j^rb’iHorM.coOTir-Jhi'^NorM 
rg C'4 r-4 eg eg r-4 eg rg ro hi hi m po hi hi 


C'jro^iTiHorHsOoOHOi-grghi^biHO 


r-4 hi ^ LT ho Qo Ch o -t-H r-4 hi ^ b*i -o 
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12645E+00 -.62S55E-01 .11V59E-13 .29226E+00 0. .14968E+01 .14946E+01 .8035SE-01 .17343E+00 


ooooooooooooo 

ooooooooooooo 


UJUJUJUJUJUJLiJliJLULLlLLlUJLjJ 

c-4LjT'^*or^h'ar^coorv.‘^’^i-« 

^rO'<»^r^L“3'00sicr‘r^o>co 

r-4i-HCr'rv.<rTH|v,Kjr^oc4^co 

rvr>.'0'0>o*obrrij“D^‘^rocMO 



o o 
o o 


ooooooooooo 

ooooooooooo 


LjJ LlS 
<1 l> 


LlI LU UJ LU 
th ro 

o ^ -.P-* 

GO o ro 
CN MD ro CK 
rj ro ^ 


Ld UJ UJ LU 

MD rx o 
«•-* no -.-I r«4 
u~> m u’i ^ 

th K) 

to o <1 


CN 

o 


ooooooooooooo 


U1 LU 
U"5 ^ 

y~^ 
CD 03 


LU UJ LU UJ 
CO O O CM 
o CO ro S3 
m <r cj 
^ ^ ^ ^ 


UJ LU LU LU 
fs OJ S3 

fs rv <p ro 

O CO S5 ^ 

■<r ro ro ro 


LU Ul LU 
CN L*3 r-4 
Ch hD O 

cn r*^ 
r-i r-j c-4 


LU 

o 

O 

o 

ro 

ro 


M 

o 


o o o o o 


*rH tH < 1 -^ 

oooooooo 


LU LU 
S3 

03 
cs o 
^ <- 



LU LU lU UJ lU 

rs o CO OJ ro 

o ro LT ro 

o o o o -r-» 

L“) L“i ti"j LT LT 



Cl. -ih -p-i -fH r-j r»4 r-4 c*^ cm cm r-j C4 r^j r^i ra Cvi r-j r-j r-4 r-j r-j r-j c-4 cj rj 

ooooooooooooooooooooooooo 

•+“ -f- H“ S" -+- ■+• ^ Hh ■+• •+' ^ ■+■ *"h “H ^ “H •+* S- -h ■+■ S** S“ “H 
LULUUJLULUUJIUUJLULLJUJLULULULUUJLUiUUJLULULUlUUJLU 
<T r-j fs CJ fs b"3 03 00 o c*^ o O th <r r-j o S3 -i-i o th M3 o rj ’•T 

OFvi-iooL*3ot"3r>4rMrsrs>Oi-tOi-iro*iT'r*ja30CDrJ<’ro*i-* 
i~<’«rL“3s3*H00rvS3>0 1J*3^r^CNf3Ch--r-f«rHCh-U‘3C0O03U’3CN^*'-* 

to CN o 00 n-h r>4 tr S3 CD o rj ^ LT 00 o --H ^ r-j rj ro CJ rj *r-i th o 

u_ 00 CK ^ ^ ^ 1-1 r-4 r-J f '4 CM r-4 ro M r<3 ro ro M ro ro to Ki ro 


OOOOOOOOOOOOO 


C‘4 CM rj rj rg rj c-4 cj r-j r-j cj C 4 -r-i th ^ -.h t-i ro ^ -th ,h ^ 

OOOOOOOOOOOOOOOOOOOOOOOOO 


lit 111 111 III m III lit til III I rl UJ LlJ 

oo-CNr^c4r^-r-iixi^r'3r^S3 
o r-3 ^ ^ CO u“3 ii") T'- o* ro 
CDoro-^^ro'^CM-r-iODiHCOCh* 
Luoroo3s3rv.oin'.-Hr^b")C-4 
cs rs u*3 <■ ro m c-4 C J 


fM. O O LT 


LU LU 
^ S3 


LU LU LU 
Kj cj C 3 
r-0 S3 


^ b'^ O O 


^■i-itMO^rs‘<i<"roC'j^c^j 


r*3‘^-lOS30»-^C30K^ 

>o c^ c^i <r S3 C'j 


C4 C3 ■<T S3 


I I 


LU LU LU LU III til I J i Li J LU LU LU LU LU 

^ Ck O i> C3 CJ Cr- Cf- CO O C> CD C4 

S3 ro o pN. o o o- ^ cci C3 h’j M3 

C4 C'J C4 f-H C-4 C3 -C" CO C4 C^' C- CO -i-* 

^ ^ D-^ ^ CP- ^ CP- bO O -Si C4 O-- 

f-o ro ^ "C" Lt"3 ii"3 '-o '-o rM CO CO cp- o> 


ooooooooooooooooooooooooo« 

-♦-*+- S- -H -+- -+• H“ S- -H H- H- H- H- *+■ -t" •+“ H- -t- -+• -t- -t“ -J- -H- -C 

lululululululuujujujlululululuujluujluluujujlululu 

pM. fH o CP- b“3 -•-» C-J C-J c-4 c-4 CN Cp- O CN r-4 PO -O C"J -O PO c-4 Cp- O 

UT LO o C-J O-' S3 b"3 C-J 00 C3 CD P3 O O S3 CK Pm |> m3 Pm i-« O 

pv,p>,rMSD'«rro-r-rcNr*^<-cjcp-rMb*3P0-»-CKC0pMrs.rMO3CP*-r-»pci 

PMCP-0-I>'0-0--’CP-Ch-COlXi03CDr-vrvi\rv.r^MDM3S3S3-<3‘<lM3rMp*^ 


poromropopopiiroropor^roho 


j I 


} I I I 


lululululululululululjlulu 

^cob"j'M"s3rv.uDs3rMK»corMrM 
COS3^POP-NfOC-4^0M30 0'PO 
b*3 C4 T-H <r CO c-4 m3 T-I CD O ^ 
-»-<po<5“PO'*rb’jr^C30b’3r^030 


OOOOOOOOOO 00000**-H*^-*MT-<b“3-r-l-t-<-rHi-H0 

oooooooooo ooooooooooooooo 

I I I I I I * !-♦- 

LULULULULULULULULULULUUJLULUUJUJUiUJLULJJ 


UJ LU LU 
Pm M3 LD 
CO C J 00 


-*J CD M3 pN -i-H 

SI n bD r*-. O 


'f-(or-s-<robDoc3s3r^cocDroor4coccicoc-JOP''5CD 


c-4 C' J C-J 04 04 


o CD M3 Csl O 00 bD c-4 SJ c-4 bD rv. CO O 


,-H-rM-rHCDS3'^C-4'«r04^'CCD 


i 1 1*1 I I 


I I I 1 I I 


oooooooooooo 


I I 


LULULUUJLULULULULU 
Cp-OOC- LDCO-f-«bDS3 
o 04 LO pM CO S3 Csi «r 

04 O- CP 04 03 r*M CP C3 O 

S3 CP ro 00 04 Pm C 4 CO -^r 

LD «r r*3 C3 C'4 c-4 -rM -rH 


} I 


1 I 


I t I o 


o o o o o o o 
o o o o o o o 


ooooooooooooo 


-4“ *+- H“ 


lulululululululululuujlululu 
CJ CP O CP- S3 «r 03 LD LD 
UD C4 ■'T S3 UD ■<r pM C> rp 
cpcprpporpcobDCprp''-HCP^r^o- 
CPrPS3UDPO-i-4CP-SirOOLDi-HLDCP 
CO ■•-4 C J C3 "'T bD UD -S3 pN CO CO CP CP CP 


UD CP Cp bD O UD ^ UD 
* ■ PO LD rp 00 


UJ LU LU 
PD S3 <T 
CP LD sr 
CO 03 


OC5 0000000CP 


o o o 


o o o o o 
o o o o o 


o o o o o 
o o o o o 


LULULuLULULULULULULU 
OS 3---epN.04COPOCP<rO 
CP PD CD C-J P^ I S3 O LD O 
rO-r-iCDS3p3-rMCDS3pD-T— » 
■^-OPMCPTHfO^-OCOO 
c^^ C-4 OJ C-4 04 ro ro P3 ro ro -«r 


ooooooooooooooooooooooooo 


S- -t- -P -4“ 

LU LU LU LU LU LU UJ 


lulululuujlululuujlululu 


000000 0-000 oooooooo c>c*ooooo 

CCOOOOOOOOOO OOC^OOOOOOCjOOOOO 

oooooooooo ooooooooooooooo 
0000000000000 0-00000000000 
C>J C J c-4 C-J C-J C'J C'J r-J C-J C'J c-4 CM CJ C'J c-4 04 04 C-4 C-4 OJ C-4 C-4 04 04 04 


PO<*bDS3P^COCPO 


04 ro -«r LD S3 


CD -r-c C J ro -O- LD -o p^ CO CP o 


CvirO'^TbDSJpMOOCPO 
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:OOOOE:.+Oi ,9J/T7tU)0 .12ti/EjE+00 . I/Si'JVE'tOI ,7654^E+01 .28901E+02 


UJUJUJLiiUJLiJLUUJliJUJLiJLULiJ 

o o c-4 o CO CD u*5 r>^ r-4 CN rs. c j «T 

T-# NO C'4 u"5 o u*i 00 th rv. o 

li") CN C 4 ^ UT <1 'O CO O UT 

^ r .1 o CO '<1 ■'T c-j T-i CD o Os 

fs| 04 f%: C-4 0-1 1-H -*-1 ^ O' O' CD 


LUUJbJUJLiJLiJUiLiJLjJUJLUUJUJ 

r^<jr^rO‘*3"CDs-tr^r-4rN.o-oo 

r'Nr-oO'rN.r'.bTU">co*3“’^h'jCKo 

oi>-cOi^CD-r-ic;-«j"ro'*r-«rp'3o 

r-j o r-4 UT pN t-h ii"3 o ^N so CO ro o 
Os ^ ^ T-i ^ c*-j c-j ro ro UT P's O' 


ooooooooooooo 

H- ■+■ -H "+■ *+* -+- -f— ■+■ ■+- »+- -f- 

LiiLULULULiJLULUUJUJLULLJLjJUJ 
O' r-4 r j r-J rj -rn u“j O' o ^ t-h rs. 
PI CO r-o CQ r-4 ut so •«r o>* r-4 o m 
r'^0'rj'*^’r'sO'‘.'Hp'j'?rNor'NPspN. 

OOCOCDCOO'O'O'O'O. O'O' 


iiJ III tjl LlI LlJ LU Lj-I I_l> nl 111 III 111 iji 

Cvl O' -rH rx O CD C4 -F-I >0 fsH T-l ^ 

CO O b“l <;> pN. Lin L"5 NO CO r-4 CO pN CD 

pN'OrooL'no<rr'^0'T-Hi-iT-io 
•«r'OCDOT-ip"j-crvjn'Ococr'0-r-i 

f-H ,-H th r *4 r -4 r -4 c j oj r -4 r-j cn ro ro 


\t\ Iti Ui LU ti! til l>l LiJ \jl 

O' T-» in tn CD ^ NO cr* o O' ci 

n ^ p's o pN ^ un sO <r C'4 un 

o>^iN^i>LncopNp'jPNO'CK 
L"5 O P‘5 -o CN ^ hO bn <i pN O' 

CO CD Pn so on un '•T p:i r-4 CO 


ooooooooooooo 

“p *+* *+" "p “p “p "p -p "p ^ 

UJUJUiLiJLiJUJLLJLLiLULULLJLULiJ 

OOOOOOOOOOOOO 

ooooooooooooo 

OOOOOOOOOOOOO 

ooooooooooooo 

:**4 r-4 r*4 r-j r-J r-4 r-4 r-4 r-4 c-4 r-4 r-4 r-4 r-4 


r-4 C4 r-4 p:i r-4 r-4 CP 

o o o o o o o 

J I i i 4 1 I 

UJ LU LtJ Ul LU LU LlI 

fp. r-4 O PD ^ CO *0 

iin C 4 in r-j o pN 'O 

c-4 Ps CD CP pN CP L"J 

ro 1-1 Pi o Pn C4 •'O 

bn p:i 1 -+ un c*4 p‘i C4 


o o o o o o o 

r j I I I I 1 

UJ LU LU LU LU LU LU 

in ^ 00 1 -H bn pN un 

r-4 CO pH ^ M C4 CP 

r-4 pn pn -r-i 00 M5 pn 

CO un m un Ps 00 


lULULULULUlUUilU 

C 4 pN O' o Pi bn bn 'C 

CO r-4 1-1 CP o P-s ro pN 

un i-« O' 00 Ph 1-1 1 -+ O' 

c>i“4i-i-sObncpoopi 

i-hCpOOPnP-.PsCDO' 


LU r-4 

O' o 


O Q. K 
CD r-4 Z II 
tn q: z 
Ci_ 

Z o Cu 

O Z CD r-4 
»— I O CD 


<c I— c I— 
<L <C 
CD CD 
LU LU iU LU 

z> o- z> => 

P- I — f— H- 
c <C -<C 
r>o^i> 

♦-H I »—* 

iIe: cc: £t: cc: 

LU Ui LU LU 
d P4 CCi Cu 


“ S? “ 55 

<L <c <j: <c 

I —I I _J 

=1 rD =3 
ta ca ta CD 
3: z 2: = 
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r j r-4 r<i r-j a c-4 

o o o o o o 

I I J I i I 

LU LjJ UJ UJ LU LU 

oci --I ro o r-J 

o r-j r-4 in 

CK CK CD 

04 O Kl «»-i hO 

ro TnJ Ijt ■•“I ho un 


o o o o o o o 

I I 1 I J I I 

l.t lU 111 Iti l«l LU 

u*^ 00 un ‘•T 

oj ro ^ CD r-4 <3 

so CO ■*-! ro ro r-j -rH 

U*D L“5 <» CO ^ 


ooooooooooooooo 


ro 04 r-4 04 04 r-j in oj c-a 04 o4 c- 4 m *>0 

thOOOOOO^OOOOOO-ip-i 

I f { ) ; J i I ) I t I I I I 

i.t 1 I I I II I tJ l_Li LU LU UJ LiJ LU m LU UJ LlJ UJ 

ir-l *<1 O ^ O MD ^ 

04 ^ -r-t bn <■ CO 

Ijt Os o CO 04 CO m 

CO 04 «rH CD CO CO Ov 

• 04 -r-l CD On 


<rH ^ in -rH Os -C- r -4 -i-H 

Os CD 04 CO O CN bn o- 


socof“i 04 ’<rm*v^^ 


CN» r-j 04 r-4 04 m 04 ca ra rrt 
0000000000 
I I I I I 4 I I i I 
UJUJUILiJLLlUiUJUJLiJLiJ 
bncrsCNinoosom 

sO^'0-r-HLn*^0040T-i 
04 CO *0 Crs CTs LT OD bn CK 

OD Os un o un 00 'O 

CM CN C-4 04 


I I I 


CK ^ Ov ^ Co -i-H 


UJ UJ UJ UJ LiJ 

CO Os bn >0 

bn ^ <1 

T-I CO sO CO CvJ 

o bn Os CO 

04 04 C'i 04 


ooooooooooooooo 

ooooooooooooooo 

Ll) LU UJ UJ LU 


ooooooooooooooo 


uiuj UJUJUJUJUJUJUJUJ 
o m tn un bn CD r-4 r-4 

^CD^TCr-MDUnW^OO 
co'Oi-^soosunO'-unr^o 
‘•T'^roor^unrom^r^ _ 

OsO^OsCh-COCOCDCOOOCOO^O^CsCKOs 


Os <r sM Os 

m un <r bn Os 

m m m 00 

KI sO CD CD 


r-4 CA r-4 r 4 c^j r -4 ^ 
0000000 

1 I I I I I I 

UJ UJ UJ UJ UJ UJ UJ 

un bn m o rs r*4 

hi O Ch- ssH C4 C0‘ 

^ i-H rs, CD Os ^ un 

CD Ch- Un -sO t-H -r-l o 

CO rv. rs r-s CO Os t-h 

I 1 I I 1 i 4 


ooooooooooooooo 


ooooooooooooooo 

ooooooooooooooo 


UJ UJ Ui U! UJ UJ 
CN CO *4D bn -T-H o 
Os CO Ck ^ bn SD 
bn r-4 CO 

- 1-1 o o r-4 bn r-4 
m m m kj ho 


•+■ H- -+- ■+* 

UIUJUJUIUJUJUJUIUJ 
bn bn r-4 CO Ch -rH r*4 

hi bn bn m m <r h'. m bn 

'*0 r-4 bn so ^ hs un 00 

i-» -«r r-4 CD o un m hi ^ 

bn bn -cr r -4 r -4 t-h w 


I I 


i 4 I ( 


I 1 


0000000 


UJ Ul Ul UJ UJ UJ 

bn hi Os rv un o 
so o o un bn hs. 
o o m hi bn o 

C>CD>.MCh UnOs 
i-H hi hi r-4 hi ^ 

I 


T-ii— t-^-rH-f-H-r-tOOOi-ir-J-r-fT-rOO 

ooooooooooooooo 


I I 


i t 


UlUJUIUIUJUJUJUJUIUJUIUtUJ 

• 5 - hi Ln CD hi Os C'‘J CD CD hs o <1 

o hi hi C 4 o Os ■>-< Os Os ^ un 1 -H 

crs un *«r hs rn r*s CD hi o o hs rs. o 

^ r>. s<) rv. o r-4 o r-. r-4 ck so 

hs. -<1 ^ r-4 sQ ^ T-t th un CD hi CD 


bn o o o o o o 


o o o o o o hi 


«pHOOOOOOOOOOOOO-tH 
-h -h \ "+■ ^ “1“ I 

UIUIUJUJUIUJUIUJUJUJUIUJ 
CDOsh 1 ‘« 5 -sQ-C-sOOQ^'*-^!^sO 

<r un r-4 Ch in bn hi bn o Ch hi CO 

-t-h r-j so Os - 1-1 CD - 1 -H bi r-4 Os Crs o 

fs^i-ii>sor^’i-<ocDr' 4 rsj-i-co 
r -4 hi r -4 -nH cci hi -«r hi hi r-4 T-f -r-i 


t j 


000000 


CO ■«r -sO bn o 

hi un r-4 r-4 r-j 

Os bn -r-H r-j 


0000000 

0000000 

“h “h ^ ^ "h 
UJ Ul UJ UJ UJ Ui Ul 

un hi CO o ck 

c»s CO hs -.-I CO Os 

o hs -sO i-v. CN r-4 CD 

O O' Ch Ch o o 

Cl -r- -r-c -.-I -rH r-J rj 


000 Ci 000000000 


I : 


M I 


I I 


UJ Ui Ul UJ Ul Ul 


Ul Ul UJ Uf ul 

o bn r-4 un -r-f SCI rn 

CO CO un <?“ so r -4 “O o hi o Ch 

cDO’C-hi'so^unhiso-i-icrs 

CD -sQ un un bn o un hi o rv. r-j 

rj hi hi r-4 bn r-4 rv so un 


i 1 


ooooooooooooooo 
ooooooooooooooo 
-h “h ■+“ “h H- -+“ -h H- -h *h *4“ 
UIUIUIUIUJUIUJUJUIUIUJUJUJUIUJ 
o o r J ‘«T r 4 hi un hi -t-H CD r j CQ ^ rv so 
-^hiOhico<“^OssOhio Oi*-icDhi 
hs un o ri hi so r-4 hi -i-H o so o CN hi 
Ci o o ch CD hs Ts. hs CO Os — 4 r-4 ^ un 

-^•«r-c-hihihirohihihi^^<r<r'«r 


ooooooooooooooo 

UlUJUIUJUIUJUJUJUJUJUiUJUJUJUJ 
oo-^'OOsps.cor-joor'sDsso-r-Hoo 
^ o 04 hs. r-j r-4 r-j r-4 r-4 rj Os r-j <r 
r -4 -r -1 hi 1 -H T-r CO hi o hi CO -•-* -rH hi rH r 4 

C 0 C 0 ‘sQh 1 Ch<"‘»^O-»-f^ 0 sh 1 '-ClCDC 0 

r-4 r -4 rj r-4 ^ T-i th -iM t-< r-4 r-4 oj rj 


o o o o o 
o o o o o o o 


T-i i-H o o o o o 

o o o o o o o 

I I -+- H- -h *h ^ 

UJUJUJUJLUUJUJUJUIUIUIUJUJUlUJ 

^ -r-H so so bn Os r-4 ^ 04 Os un sO so -*-H 

00 in CO 04 r-j rs. sO hs <r r-4 o4 co un co 

o ps. 04 un hi CO Os r-4 Os CD hi bn 04 rv. o 

<*f-HCh<»hibn-*-thst-iunhisOCh-i-H<r 

f.j r -4 1 -H ^ -I— 4 Os un un Os -r-4 1 -H -I— I 04 04 


i I r } 


{ I 


o O O O r-l 

0000000 

H h -+- H 1 H-h 

LU UJ UJ UJ Ul UJ UJ 

O O O O O ^ CO 

rr CO r-4 -o r-4 -sO o 

<r- 00 hi hs 04 r-s 

r-4 hs Os -rH hi un 

-<r sQ CP ^ tH -rH 


i-htH^OOOO^OOOO-t-*-^^ 

0000000-^0000000 


I 

lil 111 l»l III III III l;i UJ UJ UJ LtJ 

CD<r 00000 « 3-000 
o so r-4 *so r 4 CO ^ un -«3- CO r-4 
rss-^r-4rs.hico<-o^coro 
un hi --H Cfs hs -«r r 4 -f-t 04 ^ ps. 

^ ^ ^ 00 so ^ r -4 p^ r -4 ^ so 


o o ^ CO 
so r-J *<1 o 
hs r 4 p^ 


1 i 


\ I 


-^OOOOOOOOOOOOOhl 


0000000000000 


UIUIUIUJUIUJUIUIUJUJUIUIUI 

hi ps. so un so o ro ^ hi o so un so 

ps* un r -4 CG r-J o o r -4 co r -4 un 

r- 4 hi-»-ro<"Osrs^osrssDs^o-»-t 

-,-1 hi *<1 CO 00 o r-j ■«r r -4 '*0 00 00 ^ 

r -4 -r-i C 4 hi un sc sQ so un ^ ro r-J 

I 


r -4 r -4 c-j r -4 cm r -4 r-j r -4 r-j r -4 r-j r -4 r -4 r -4 r-J 


o -rH r-J hi un so 


r- 4 hi-«runsops 00 i>-o 


r-j hi ‘<r un so 


c-4 hi un so hs. CO CN o 


r-j hi <r un >o 


5 

(0 

o 

VO 

CM 

CM 
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5973E+00 -.97070E-01 0. -.1510BE-03 .27701E-02 0, .19344E-18 .6B435E-17 


I 

I 


o o o o o 
o o o o o 


m r* 4 ' r -4 oj n bi n c-i a c-j C 4 ho a 
O OOOOO'r-iOOOOOO-r -4 
I I I i I I t I r I I j I j 

m ill 1 1 1 111 111 < 1 1 < t i lit til til til 1,1 1 1 1 [1 1 

ro CK rv c j r '4 MD <r ^ ch t-* i-i m 

so sD<i-r-iroLTCCiro*»'-»-irx[^rNT-H 
rv. h"jiv.Ma-'*oo-*o<>-rorxMr^r-4 

mooot>-^uT-.-Ha*-ocoocio^co 

a 3 ^rj'^<r^TH-<r'<T<rr' 4 ^oor^ 


•«r ^^roroc-Jin oi m ^ 

O 00000*^000000 t-h 
I I I I J (( I I I I I I I 

1 1 1 ilf iti ti« Ul itl Iti 111 III itl iir iti 

tn o3CCi^rv<j^cDu"5 0 '^rv.cor^ 
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Figure 46 Sample Printed Output for BSF Method {Sheet 14 of 14) 
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KRR(1)= 0 KRR(2)= 0 KRR<3)= 0 KRR(4)= 0 KRR(5)= 0 

RZNEU(1)= 100.00000 RZNEU(2)= 100.00000 RZNEU(3)= 900.00000 RZNEW(4)=1000. 00000 RZNEU(5)=1000.00000 
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-r-TT-fOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO-rSi-^ 

ooooooooooooooooooooooooooooooooooooo 

LiJ LjJ LU LU LlJ lU LiJ LU UJ UJ UJ LiJ LU LlI LU LlJ LU LjJ m UJ LtJ LiJ LiJ 111 l.l lit III LlJ UJ LiJ LiJ IjJ 14 ! lil I |J 

NOL-jo-^cor4 0Nr-JrNOOboorNrorN..,-ecc-<‘CD-.-Hr'xrors.ob’5oorN.r-4CrNr-4co-«robTNO 
X pN.-<rrN.cDrN.r- 4 CDU' 3 CciTNvcoch-<i>pNsi~sOvrjrN.Norvr -4 0 NT-irN.oc>NCorN%cob" 5 cor- 4 iNs.corN.-« 5 j“rN, 
0 NS 3 r 0 -o 03 ONrNs«TCDOCKb ’3 0 NC 0 U“jrN.rN.ri^r* 4 rN.rsb '3 000 NbT 0 NOC 0 ^rs, 0 NC 0 --or 0 N 00 N 

On ro pN. 1-s b“) On ro S 3 o r 4 bT rs On -rs r -4 ro ^ ^ -«r ro r -4 -r-< CK p'v b *3 r -4 o NO ro On b^ T-i r-. ro On ^ 

CD is 1-f 04 C ’-4 r -4 ro ro ^ ^ «r ^ ^ b "3 in b "3 bo bo bo bo bO m -<r ^ ^ ro ro r -4 r -4 r -4 1-s is 00 ^ 


o o 


000000000000000000000000000000000000000 


UJ LU LlJ LlJ LiJ UJ UJ UJ LiJ LiJ LjJ 1 _U UJ UJ LiJ LlJ LlI UJ LU UJ UJ 111 LjJ US It! <1 1 111 l.l Ut »«t Lil >4 1 UJ Uf l« > It? U 1 lil 
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ooooooooooooooooooooooooooooooooooooooo 
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5 

in 

o 
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SINGULARITY LOCATION IN PHYSICAL SPACE AT R= 1.000 

XSIN/Z= .59424 YSIN/Z= 0.00000 -YSIN/Z= 0.00000 

SINGULARITY LOCATION IN MAPPED SPACE AT CURRENT STATION 

S0= .2747E+00 0. Sl= .2747E+00 0. S2= 0, 0. 

SINGULARITY LOCATION IN MAPPED SPACE AT NEW STATION 

SON= .2747ET00 0. S1N= .2747E+00 0. S2N= 0. 0. 


XTIP= 1.00928 YTIP= 0-00000 ZTIP= 1.67734 OMTIP= 0.00000 PSID= 33.00000 RTIP= 2.00000 


SINGULARITY LOCATION IN PHYSICAL SPACE AT R= 2.000 

XSIN/Z= .62557 YSIN/Z= 0.00000 -YSIN/Z- 0.00000 

SINGULARITY DERIVATIVES AT CURRENT STATION BOPR= 0. 
SINGULARITY DERIVATIVE AT CURRENT STATION S2PR= 0. 
SINGULARITY DERIVATIVES AT NEW STATION S0PRN= .6159E-02 

SINGULARITY DERIVATIVE AT NEW STATION S2PRN= 0. 

PSIMU= 36.031879 
SHOCK I 


R= 0.000 


I 

X 

B 

BPR 

BSEC 

BR 

BRR 

BRTH 

7 

-.15708E+01 

.20899E+00 

0. 

-.24553E+00 

0. 

0. 

0. 

3 

-.13464E+01 

.20281E+00 

-.22043E-01 

.49070E-01 

0. 

0. 

0. 

4 

-.11220E+01 

,19910E+00 

-.12548E-01 

.35555E-01 

0. 

0. 

0. 

5 

-.89760E+00 

.19718E+00 

-.52650E-02 

.29357E-01 

0. 

0. 

0. 

6 

-.67320E+00 

.19673E+00 

.15416E-02 

.31309E-01 

0. 

0. 

0. 

7 

-.44880E+00 

.19787E+00 

.93857E-02 

.38603E-01 

0. 

0. 

0. 

8 

-.22440E+00 

.20095E+00 

.125B8E-01 

-.10065E-01 

0. 

0. 

0. 

9 

-.71054E-14 

.20352E+00 

-.17316E-11 

-.10213E+00 

0. 

0. 

0. 

10 

.22440E+00 

.20095E+00 

-.12588E-01 

-.10065E-01 

0. 

0. 

0. 

11 

.44880E+00 

.19787E+00 

-.93B57E-02 

.38603E-01 

0. 

0. 

0. 

12 

.67320E+00 

.19673E+00 

-.15416E-02 

.31309E-01 

0. 

0. 

0. 

13 

.89760E+00 

.19718E+00 

.52650E-02 

.29357E-01 

0. 

0. 

0. 

14 

.11220E+01 

.19910E+00 

.12548E-01 

.35555E-01 

0. 

0. 

0, 

15 

.13464E+01 

.20281E+00 

.22043E-01 

.49070E-01 

0. 

0. 

0. 

16 

.15708E+01 

.20899E+00 

0. 

-.24553E+00 

0. 

0. 

0. 


R= 

0.000 







I 

X 

C 

CPR 

CSEC 

CR 

CRR 

CRTH 

2 

-.15708E+01 

.49170E+00 

0. 

-.64481E-01 

0, 

0. 

0. 
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3 

-.13464E+01 

.49008E+00 

-.26059E-01 

-.16777E+00 

0. 

0. 

0. 

4 

-.11220E+01 

,48001E-f00 

-.5B190E-01 

-.11860E+00 

0. 

0. 

0. 

5 

-.89760E+00 

.46397E+00 

-,763B1E-01 

-.43527E-01 

0. 

0, 

0. 

6 

-.67320E+00 

.44573E400 

-.75266E-01 

.53456E-01 

0. 

0. 

0. 

7 

-.44880E+00 

.43019E+00 

-,51272E-^01 

.16040E+00 

0. 

0. 

0. 

8 

-.22440E+00 

.42272E+00 

-.27591E-02 

.27197E+00 

0. 

0. 

0. 

9 

-.71054E-14 

.42895E+00 

.66030E-01 

.34112E+00 

0. 

0. 

0. 

10 

.22440E+00 

.45235E+00 

.13259E+00 

.25214E+00 

0. 

0. 

0. 

11 

.44880E+00 

.48846E+00 

.17371E+00 

.11434E+00 

0. 

0, 

0. 

12 

.67320E+00 

.53031E+00 

.18394E+00 

-.23211E-01 

0. 

0. 

0. 

13 

.B9760E+00 

.57101E+00 

.16464E+00 

-.14880E+00 

0. 

0. 

0. 

14 

,11220E+01 

.60420E+00 

.12036E+00 

-.24586E400 

0. 

0. 

0. 

15 

.13464E+01 

.62502E+00 

.58153E-01 

".30853E+00 

0. 

0. 

0. 

16 

.15708E+01 

♦63030E+00 

0. 

-.20976E+00 

0. 

0. 

0. 


R= 0.000 


4 


I 


J 

x/z 

Y/Z 

H 

UI 

VI 

UI HI 

H2 

2 


2 

-.21243E-14 

-.24084E+00 

.28240E+01 

.84953E-15 

-.31799E+00 

.94810E+000. 

0, 

3 


2 

•13377E+00 

-.21751E+00 

.28100E+01 

-.10254E+00 

-.30488E+00 

.94686E+000. 

0. 

4 


2 

.26156E+00 

-.19286E+00 

.26321E+01 

-.19701E+00 

-.30596E+00 

.93144E+000. 

0. 

5 


2 

.38025E+00 

-.16526E+00 

.23176E+01 

-.27148E+00 

“.32215E+00 

.90693E+000. 

0. 

6 


2 

.48486E+00 

-.13325E+00 

.19129E401 

-.31259E+00 

-.35851E+00 

.87963E4000. 

0. 

7 


2 

.56920E+00 

-.95849E-01 

.14827E+01 

-.29623E+00 

-.42460E+00 

.85555E+000. 

0, 

8 


2 

.62703E+00 

-.51972E-01 

.11328E+01 

-.16994E+00 

-.51635E+00 

.83935E+000, 

0. 

9 


2 

.64941E+00 

-.17264E-14 

.10012E+01 

.87156E-01 

-.54257E+00 

.83548E+000, 

0. 

10 


2 

.62703E+00 

,51972E~01 

.11328E+01 

.31854E+00 

-.42555E+00 

.84702E+000, 

0. 

11 


n 

.56920E+00 

.95849E-01 

.14827E+01 

.40136E+00 

-.28632E+00 

.87002E+000. 

Or 

12 


2 

.48486E+00 

.13325E+00 

.19129E+01 

.3855BE+00 

-.2015BE+00 

.90039E+000. 

0. 

13 


2 

.38025E+00 

.16526E+00 

♦23176E+01 

.32204E+00 

-,i5747E+00 

.93354E+000. 

0. 

14 


2 

.26156E+00 

.192B6E+00 

.26321E+01 

.22991E+00 

-.13779E+00 

.96341E+000. 

0. 

15 


n 

.13377E+00 

.21751E+00 

.28100E+01 

,11933E+00 

-.13531E+00 

.98359E+000. 

0. 

16 


2 

.11566E-13 

.24084E+00 

.28240E401 

.10086E-13 

-.14852E+00 

.98891E+000. 

0. 


J 


HX HY 

HR 

HIX 

HIY 

HIR H2X 

H2Y 

H2R 

2 

2 

0. 

-.37363E+00 0. 

0, 

0. 

0. 0. 

0, 

0. 

3 

2 


19198E+00 -.3919 

2E+00 0. 

0, 

0, 

0. 0, 

0. 

0. 

4 

2 


49237E+00 -.34407E+00 0. 

0. 

0. 

0. 0. 

0. 

0. 

5 


♦ 

71920E+00 -.24384E+00 0. 

0. 

0., 

0. 0. 

0. 

0. 


6 2 -.83490E+00 -.11168E+00 0. 0. 0, 0, 0. 0. 0. 

7 2 -.78004E+00 .36485E-01 0. 0. 0. 0. 0. 0. 0. 

8 2 -.48149E+00 .18161E+00 0. 0. 0. 0. 0. 0. 0. 

9 2 -,43627E-11 .25891E+00 0. 0. 0. 0. 0. 0. 0. 

10 2 ,48149E+00 .20573E+00 0. 0. 0. 0, 0. 0. 0. 

11 2 .78004E+00 .43174E-01 0. 0. 0. 0. 0. 0. 0. 

12 2 .83490E+00 -.15495E+00 0. 0. 0. 0. 0. 0. 0. 

13 2 .71920E+00 -.35036E+00 0. 0. 0. 0. 0. 0. 0. 

14 2 .49237E+00 -.50791E+00 0. 0. 0. 0. 0, 0. 0. 
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19198E+00 -.58903E+00 0 
-.56B34E+00 0 


Oooooooooooooooooooooooooooooooooooooooooo 

-f- H- -t* “H *+• “t- - 4 “ ^ - 4 " ^ ■+■ -t- lllllllllllltfllflllltll 

111 111 III III III Itl 111 III 111 III Itl lil 111 Ijl LiJ 111 LiJ LiJ LiJ lil lil 111 1 1 1 1 1 1 1 1 ] I J I LiJ 1 iJ LU LU LiJ Ll] LlI LiJ UJ LiJ LlJ LiJ LxJ LlJ LlJ LlJ 

Mo~.r^T-^oL"5-^o-^-*rmc>'i>'T-iT— i'orj*riXiC)OMrric'jcDr4ror<>cncioroc<'0-'i>-'0 r-j rj o kd lt 

o •«T- o I rj Pv. m <j- C 4 T-i O' fo 'T <r T 4 M r--. m 'O o CD -<3 T'J 1-1 -o T 4 ~o C J o ro ~o T'- LT O' M 'O 

0 'OC 0 t'ii' 0 'O 0 ''«i»H 00 rii 0 'Cci' 0 ' 3 -f*ir-ji-ir- 4 iiT-«rc 0 rop'i( 33 i-irvij'ir 0 T-tr'..«a-t' 0 <r' 0 C 0 i-i<‘r 00'^'0 

iH th «o i 4 T -V w Ki ro r -4 T '4 1-1 1-1 1-1 »-( iH T-i 1-1 O' rv io iiT UT •«T po ro rj fii r -4 r -4 1-1 1-1 T-i 1-4 1-1 i-< r -4 r -4 r -4 1-t 1-1 O' 


CM PO LO iO rv 03 O' O' O O O 1-1 1-1 1-1 1-1 1-4 Til C-4 Cil fil (N CM Oil r-4 f^4 CM 04 C-J C-4 r-4 0-4 r>4 fil Cil r-4 CM r-4 04 CJ r-4 0-4 04 

i»-H T-^ ^—4 tH tH t— 4 ^4 ^—4 i-H t-H t—4 

<roDcococoaDr^oorvCDi>^cot>.in’«rc- 4 CDC>**<i>o<iu‘ 5 ininiJ' 3 Ln^<’^'« 3 “^‘<'rororor-or-jr- 4 rsii>JC‘ir -4 


LU LjJ LJ 
O CO <1 
C -4 CK t> 
b "3 CD 

ro CN 


III 111 IjJ IjJ LU LiJ III l«l III Itl Itl Itl lit |J> liJ IjI LU LtJ LU LiJ LU LU LjJ UJ LU LiJ UJ LU LlJ LU 

CN cs L“5 fo u“i ^ cs i> ^ ro cs i> -wH rv. u“5 o> ‘O', CD ro rj LT o ro iiT 0'4 cr. ro bT 

o CK r-4 C J CO r-j CO o ^ uri T-i c^j u"> o 1 -H CO CD r-4 *rH -T o CD c^- r-’j rv. LT UT rv 

•«5- r-4 o !>• r-j r- r-4 ii"5 sD (V. r-j o r-j r-4 ro *«3“ UT ro "Q* CN rv Ch- o CD rv <“ o 

i-iif-Hor^coo^csr^coLiT^*'0<i^r-4ro^'«3‘h'Dr-4^*^r-4roro<r<iNO'C'r-4'*-H 


r-4 04 04 04 IP-4 04 04 04 C-4 04 O'! 04 PO r<3 ro PO m PO PO PO P*3 PO PO PO PO PO PQ PO PO PO «r IT IS- IT P<3 m r*3 PO PO PO 45- ■<r 

oooooooooooooooooooooooooooooooooooooooooo 
t I I 1 I I I I I 1 I I I I I I I I t I I I I I I I I I 1 I I I 1 f I I I I I I I I 
LLlLiJLiJLUIXlLlJUJlJjUJIjJUJlJJLlJLlJUJUJUJUJUJLlJUJLiJlxIUJLlJUJLlJUJLiJLiJLLjLiJUJUJLiJUJLLJUJLiJLiJ UJ ' UJ 
Oi' 000 O'CDr''ij' 3 L "3 0 '''Tp 0 i-i'>Tr- 4 C 0 C 0 r 00 '^p>icrii-ii- 4 r* 3 ' 0 a 3 -rt 04 r'-cauTpnpnc>O 0 'ri 4 ii’Ji-iP 0 04 iT 
0'iTOorsi04r^'0aD0'Ki'00'i-ir^0'0'0'0'04'0 04 rM0'i-iir^r0 04 '0-«Tp0'0'^ 04 04 u-3 0'0'OOC0 
aD' 43 i- 4 ' 0 P 0 i-iC 0 ' 0 mMi-io 0 '' 0 aD-«TOi-iO 04 ' 0 P 0 Ob' 3 i-iP'.ii -3 4 TP 0 i-i' 00 '' 0 r^Oi-ip 0 '^p 0 i-ia 3 i-i 
poror*ir -4 04 r'ii-ii-ii-ii-ii-iiHO'CDr''r'ir'''Oin'rpopopocM 04 i-i'-Hi- 4 i-ii-iO'r''r'. 03 i-ii-ii- 4 i-ii-i,-irMLT 


O' O' P'. oo 1 -I 1 -i 


04 04 04 04 04 04 04 04 .04 04 C-4 04 O' 04 04 rj 04 Til 04 04 04 C-4 04 Til O'! 04 0-4 04 04 C4 0-4 0-4 0-4 O-J 0-4 04 0-4 Osl 04 C-3 C-J 0-4 

T— 1 ▼— ! n-H i I if-H t*H t-H i— I i-H t— H ^—4 i—i t— 4 t-H 

TO r-4 r-j r-4 r-j r-4 r4 r-4 ur» ^ r*4 r-4 >o b") <“ ro r-4 r-4 ^ <i o ut l*:i ut in ^ ^ ^ ■«r -«r ro lo lo ro ro r-4 c-* r-4 


LU LU LU 

in b*:i <■ 

-rH CO CN 
<1 o ro 
ro ra -.-4 


LULULUlULULULULULULULULULULUUJlULULULUlUiULULULULU LU LU UJ 
r^'«r-.-Hb")Ln'«rr -4 0 -»“*CKbTrvCDor' 4 oro'or- 4 ^i-iCsi>-ro<r-ocn'^ 
■<r'^r^‘<rr>vCM‘C-b')cor -4 0 ororN.ioocDb' 3 CKi-ico-i-i-«rcKb“ 3 rvrar- 4 r^ 
ch‘ 00 *oinro'<jfsi^-i-«tiT>o*«rT-ioo-^c- 4 r%JcooTor- 4 r^<roKiro 


rsiro^rm'OP^coD^o 
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1321E-03 2 2 ,3490E-04 -.4541E-01 11 12 .6445E-02 


o o o 

-H 

LU UJ LU 
irt m 03 
CO CO o 

o ro ■«t 
o o ^ 


o o o 
r i -+■ 

U 4 Ui UJ 
'O 

^ *<5 CS 
m \n -c 

^ CQ 

^ ^ 


o o o 
■+■ •+- 
Ui UJ 1 x 1 

o cs 
rs CN ro 
00 00 o 

CM CM eg 


O 

<o 


I I 


I i i 


I 1 I I 


_lxlUJliJLiJUj ____ 
ror^OMDc-4corx-.-iro'Ooo 
U’3*<JljTOO'«3“CM-«rr0-0'0 
<rrjir3^aDL*>rn-.-H03 
^ TO r-4 CM i-» ’ 


r I 


Li I 111 liJ LiJ liA 


I I I I 


o o o o O o 


r Lxl UJ Lxl LU Ul UJ Ld 
i>o"sOi-HO<-o-io-r-icDO^rv.o 
'O i-t<;>’»“iC4‘^^OJ‘OiiTCNr^oro»-«orx 
car*4ToincocNr^!oo^^ro’i-i&vtooch.o^ 
c^oofv.^oLl^«r^oroc^Jr^<-r-4‘l-l^THa^oD^^m•^ 


•o 

o 


U 

CO 

CD 


T'4 C'J CN r-j r*j r-4 w 


«»-iiOToroioro roToiororororororoTororoiorororoToro 


^ ^ ^ ^ ^ oi C'4 r4 r-4 0-4 c-J CM CM r-4 r*J r-4 c-4 r-4 cj TO to ro ro ro TO TO TO TO 

ooooooooooooooooooooooooooooo 


Ld Lxl U4 

o m cs 

CO b“J -r-i 


I I 


i I 


I I 


TO r -4 

CK 

CO <■ 


TO TO r-J t-H 


TO m CK 
r -4 r j 
-f-i 00 -c 
T-l CK CO 


Ld Ld Ld Lxl 


O' r-4 CK O^ 


rv MD uo TO ro CM r-4 

r r I* .* 


Ld Lxl Ld Ld Lxl 

i^^i-iCKC-j-rHor^o-roTHCKcorvsoci 
-rror>-.ooror-4r*4CM'Cooir>r‘4CKLr3 
<3r-4G3U“DrO-^^CK-^*-OrKCKTOr-sTO 
. MD UT TO TO CN CM 


I I 


I 1 


o o o ^ 

3 

o o o O 

o o o ^ 

Ld Ld lxl Si 

O- tH MD 3 

00 bO 00 

o TK r -4 w 


I I I 1 I 
Ld lxl Ld Ld Lxl 
aOMOCM^CK 
GO ^ Tv ^ -*0 
CK b^ TO 


0‘irtminb^inbTirtboir3bTb">mb“DLn-^'-o>o 

oooooooooooooooooooooooo 

I I I f I I J I I I I I I t I I I I I I I I I I 
LxIlxJLULdLULdLdLdLdLdUJLdLdLxILdLdLxILdLiJLdlxILLlLdUJ 
ooTo-MDrocKT-4 cM^-i^i»-Hro»-irv^^r-4rKro-^o“rKfn-<r'<j 


_ - CK*o>ocKrKCOTK...^-^-.--w-w 

,CKo-or>^r-4ocsmr0T0in0DCM<jrsj 

oorKMDin^r^rTor-jr-JCM-t-i-f-i-i-HTHCDrv.sob'D^-TorotMCN 


^-lO-'O^TO^^OO 


rvi ^ 


TO r4 

CM r>i 

-•-t ^ 

TO 
O TO 
CK TO 
00 TO 


Ld Id Id 
TO FK -<J 

^ fK in 
r-j to -i-H 


01 


TOTOToroToroTororoToroTOTorororororoTororoTororoTOTOTOTO 


I I 
lxl Ld Ld 
PK fK TO 
00 S 3 CQ 
Tk CK ^ 

>o m m 


^<-^^*^<^<'^-^mminmu^ir3b"3iriin 

oooooooooooooooooo 


III] 


I I 


Ld Ld LU Ld Ld Ld lxl 

TO CO CM CO CK CO rs 

^ ^ r -4 m CO •<) 

04 PK TO CK ^ TO 

TO C ^4 r-J ^ ii-i 


I I 


I I 


IdLdlxlLdLdldLdLdLxl 

cKr-4bOCDrKbooorK.ro 

o«a‘CMi-<ooo*^CKro 

OOCMCKODOOi-t^OO^ 

O^00<Jb0’^‘^TOCK4CN 


I I I J I 


^ rj o in 


I I I 


I I 


pK ro 

so ro 


TO TO 
r si TO 


•-I CM TO ^ w 
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19718E+00 -.11182E+00 .38025E+00 -.16526E+00 0. ,12625E+01 .12357E+01 *25876E+00 .43205E+00 


o o o 


o o o 


oooooooo-»-<roro 

oooooooooo^ 


lit i_> J UJUiLlJ LU LU UJ UJ LxJ UJ 
roro-.-HOQr^'Oocs-i-Hr^'O 
•*ocKc>'«rrv.rv.cNCKinr-j'»-H 

o in>Of-^rHco-*oooo'Oin 
o o <jc;C'J-.-ic>'Cnoqodco 

inininMCMTH-rH-i-Hf-t-H 


O O *r-« 


O O 


UJ L±i LU UJ UJ UJ UJ UJ UJ UJ LlI 

oo*oQO*«ot^oo^&^roco 

CNJoin-'O^^ro^oor*'--^ 

CM CM^CMin-O'O^ro^rs 


ooooooooooo 
^ ^ ^ "4“ ^ ^ ^ H“ -4- 
ii I UJ UJUJlllUJUJUJLlJUJUJ 
CD'-HCKr^coi>''00»-Hor^ 

orv-mMSCJCoroch-^rvoo 


ooooooooooo 

^ ■+“ -+- •+“ •+“ H- -H ^ “K 

iti l>l 111 111 jil til Ltl UJ liJ LU LU 

rKincor»s^-r-iij*5 03ros3rv 
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Figure 47 Sample Printed Output for BSC Method (Sheet 14 of 14) 
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INTERPOLATED Z STATIONS SPECIFIED AS INPUT 
21 -^ 5.0000 


I XI YI CPI 


1 

-.00000 

-.39603 

.07728 

*;) 

.66285 

-.37143 

.07778 

3 

1.30799 

-.33800 

.08309 

4 

1.91370 

-.29407 

.09831 

5 

2.44957 

-.23809 

.12895 

6 

2.87706 

-.16922 

.1855? 

7 

3.15500 

- . 08784 

.30383 

8 

3.24704 

- . 00000 

•7 ‘=■,'70'=; 

« .w tJ U vJ 

9 

3.15500 

.08784 

-.10914 

10 

2.87706 

. 16922 

-.13759 

11 

2,44957 

.23809 

-.12083 

12 

1.91370 

.29407 

-.10324 

13 

1.30799 

. 33800 

-.09035 

14 

.66235 

.37143 

-.08015 

15 

o 

o 

o 

o 

o 

.39603 

-.07548 



ZI = 

10,0000 


I 

XI 

YI 

CPI 

1 

-.00000 

--.35294 

♦05442 

O 

1.32423 

-.35352 

,05339 

1 

iJ 

2.61471 

-.34361 

.05582 

4 

3.82715 

-.31877 

.06462 

5 

4.90017 

-.27327 

.08654 

6 

5 . 75615 

-.20285 

.13149 

7 

6.31274 

-.10753 

.23903 

B 

6.49408 

“.00000 

-.12945 

9 

6.31274 

.10753 

-.21065 

10 

5.75615 

.20285 

-.19996 

11 

4.90017 

.27327 

-.16707 

12 

3,82715 

,31877 

-.13311 

13 

2.61471 

.34361 

-.11060 

14 

1.32423 

.35352 

-.09730 

15 

.00000 

♦35294 

-♦ 09219 


Figure 49 Spanwise and Chordwise Interpolation Output (Sheet 1 of 6) 
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CPI 


zr-" 15.0000 


XI YI 


1 

"t 00000 

-.21070 

.04085 

1 

1.98525 

-.23922 

.04068 

3 

3.91938 

-.25920 

.04149 

4 

5.73572 

-.26562 

.04637 

5 

7.34247 

-.24928 

.06147 

6 

8.62395 

-.20021 

.09557 

7 

9.45795 

-.11299 

.16495 

8 

9.74111 

-.00000 

-.32318 

9 

9.45795 

.11299 

-.29193 

10 

8.62395 

.20021 

-.2357/ 

11 

7.34247 

.24928 

-.19094 

12 

5.73572 

.26562 

-.14698 

13 

3.91938 

.25920 

-.11968 

14 

1.98525 

.23922 

-.10640 

15 

*00000 

.21070 

-.10221 



Z3> 

20 . 0000 


I 

XI 

YI 

CPI 

1 

-.00000 

-.00340 

.04819 

n 

2.64587 

-.06772 

.03052 

3 

5.22303 

-.11878 

♦03064 

4 

7.64268 

-.15732 

.03403 

5 

9.78240 

-.17672 

.04287 

6 

11.48870 

-.16412 

.06683 

7 

12.59977 

-.10371 

.09909 

8 

12.98815 

-.00000 

-.34928 

9 

12.59977 

.10371 

-.34131 

10 

11.48370 

.16412 

-.26006 

11 

9.78240 

.17672 

-.20577 

12 

7.64268 

.15732 

-.15503 

13 

5.22308 

.11878 

-.12645 

14 

2.64587 

♦06772 

-.11334 

15 

.00000 

.00840 

-.09520 


2246-043W 


Figure 49 Spanwise and Chordwise Interpolation Output (Sheet 2 of 6) 
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TRAILING EDGE LOCATION » XTE=- 


9.72310 


ZTE= 


25.00G00 



ZI = 

25.0000 


I 

XI 

YI 

CPI 

1 

-cOOOOO 

--.00840 

.08614 

2 

3.30727 

-.00697 

.10399 

3 

6.52S33 

-.00554 

. 10487 


9.55182 

- . 00722 

.06062 

cr 

J 

12.22453 

-.05786 

.03136 

6 

14.35473 

-.03715 

.04167 

7 

15.7415? 

- . 07489 

.04227 

8 

16.23519 

-.00000 

-.34990 

n 

7 

15.74159 

.07489 

-.37036 

10 

14.35473 

.08715 

-.27660 

11 

12.22453 

. 05786 

-.21016 

'I n 

X ^ 

9.55182 

. 00722 

-.13860 

13 

6.52838 

.00554 

-.07414 

M 

3.30727 

.00697 

-.05247 

15 

.00000 

. 00840 

-.06078 


2246-043W 


Figure 49 Spanwise and Chordwise Interpolation Output (Sheet 3 of 6) 
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CHORDUISE INTERPOLATION PROGRAM 


XSPAN= 

UPPER SURFACE 


ZUP YUP CPU 


0*00000 

0*00000 

*17370 

1*97570 

*31081 

-*04595 

2*97885 

*35563 

-*06174 

3*98170 

*38216 

-*06950 

4*98430 

*39596 

-*07541 

5*98664 

*40013 

-*07993 

6*98875 

*39667 

-*08361 

7*99063 

*38700 

-.08678 

8*99230 

*37220 

-*08962 

9*99376 

*35308 

-*09218 

10*99504 

* 33030 

-*09451 

11*99614 

*30435 

-*09665 

12*99703 

*27560 

-*09872 

13*99787 

.24434 

-*10048 

14*99852 

*21075 

-*10220 

15*99904 

*17492 

-*10400 

16*99945 

*13687 

-*10584 

17*99974 

*09653 

-*10779 

18*99992 

.05378 

-.10991 

20.00000 

*00840 

-*09520 


0*00000 

LOWER SURFACE 


ZLO YLO CPL 


0*00000 

0*00000 

*38624 

1*97570 

-*31081 

*09929 

2*97885 

-*35563 

*09352 

3*98170 

-*38216 

*08517 

4*93430 

-*39596 

* 0773 P . 

5 * 98664 

-*40013 

*07116 

6*98875 

-*39667 

*06607 

7*99063 

-*38700 

.06173 

8*99230 

-*37220 

*05789 

9*99376 

-*35308 

*05444 

10*99504 

-*33030 

*05131 

11*99614 

-*30435 

*04845 

12*99708 

-*27560 

*04585 

13*99787 

-*24434 

*04330 

14*99852 

-*21075 

*04085 

15*99904 

-*17492 

*03847 

16*99945 

-*13687 

*03602 

17*99974 

-*09653 

*03343 

18*99992 

-*05378 

*03061 

20*00000 

“*00840 

*04819 


2246-043W 


Figure 49 Spanwise and Chordwise Interpolation Output (Sheet 4 of 6) 
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XSPAN= 2,00000 


ITEU= 26 ITEL= 26 
ZTE= 21,028^*8 


UPPER SURFACE 


LOWER SURFACE 


ZUP 

YUP 

CPU 

ZLO 

YLO 

CPL 

3,07973 

“,00000 

,33169 

3,07973 

-.00000 

,3316"^’ 

3,17301 

,07431 

“,03599 

3,17301 

-.07431 

,34399 

3.46086 

,14408 

-♦09853 

3.46086 

-.14408 

.22219 

3,47890 

,14825 

-.10180 

3.47890 

-.14825 

.21518 

4,08379 

,22248 

-.10538 

4.08379 

-.22248 

,14268 

4,56970 

,25828 

-.10689 

4.56970 

-.25828 

.11944 

5,22535 

,29749 

-,10540 

5,22535 

-.29749 

,09610 

5.64143 

,31330 

-,10583 

5,64143 

-.31330 

.08898 

6,68999 

*34133 

-.10489 

6.68999 

-.34133 

.07542 

7,64750 

.35593 

-.10331 

7,64750 

-.35593 

,06621 

7,73635 

,35654 

-.10338 

7,73635 

-,35654 

.06566 

8,75433 

,35710 

-,10383 

8,75433 

-.35710 

.06027 

9*77240 

,35048 

-,10426 

9,77240 

-.35048 

.05581 

10,79047 

,33787 

-,10469 

10.79047 

-.33787 

.05203 

11,80849 

,32021 

-.10510 

11,80849 

-.32021 

.04878 

12.82645 

,29824 

-.10560 

12,82645 

-.29824 

.04598 

13.84431 

,27256 

-.10597 

13.84431 

-.27256 

.04338 

14,86206 

,24363 

-.10643 

14,86206 

-.24363 

,04099 

15,11175 

,23589 

-.10657 

15,11175 

-.23589 

.04044 

15,86684 

,21173 

-.10754 

15,86684 

-.21173 

,03381 

16,86733 

,17727 

-.10889 

16.86733 

-.17727 

,03661 

17,86771 

,14037 

-.11039 

17.86771 

-.14037 

,03428 

18.86799 

,10099 

-.11208 

18.86799 

-.10099 

,03172 

19.86817 

,05900 

-.10969 

19.36817 

-.05900 

,03385 

20,86826 

,02727 

-.10037 

20.86826 

-,02727 

.05001 

21,02848 

,02407 

-.09760 

21.02848 

-.02407 

.05389 


XSPAN= 

4.00000 

ITEU= 22 

ITEL= 2 

ZTE= 

22,05696 



UPPER SURFACE 


LOWER 

SURFACE 


ZUP 

YUP 

CPU 

ZLO 

YLO 

CPL 

6,15946 

-.00000 

♦19495 

6,15946 

-.00000 

,19495 

6,33707 

,09505 

-.14306 

6,33707 

-.09505 

*28337 

6,92546 

,18444 

-.16764 

6,92546 

-.18444 

.16423 

6,94887 

,18780 

-.16839 

6,94387 

-.18780 

,15984 


2246-043W 

Figure 49 Spanwise and Chordwise Interpolation Output (Sheet 5 of 6) 
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8*05636 

,26235 

-,15583 

8,05636 

-.26235 

♦10335 

8,16225 

,26849 

-,15452 

8,16225 

-,26849 

,09875 

9,14929 

,29552 

-,14607 

9,14929 

-,29552 

,08085 

10,23871 

,31443 

-,13663 

10,23871 

“,31443 

,06517 

10,45198 

,31673 

-,13476 

10*45198 

-*31673 

,06253 

11.29573 

,-31493 

-,13202 

11,29578 

-,31493 

,05720 

12,34501 

,30629 

-,12871 

12,34501 

-,30629 

,05177 

13,39427 

,29170 

-,12547 

13,39427 

-,29170 

,04723 

14,44349 

,27193 

-,12244 

14,44349 

-,27193 

,04345 

15,30941 

,25205 

-,12011 

15,30941 

-,25205 

,04081 

15,48715 

,24736 

-,12005 

15,48715 

-,24736 

,04041 

16,5047? 

,21780 

-,11983 

16,50477 

-,21780 

,03819 

17,52230 

,18526 

-,11982 

17,52230 

-,18526 

,03597 

18,53974 

,14988 

-,12010 

18,53974 

-,14988 

,03361 

19,55709 

,11169 

-,12055 

19,55709 

-,11169 

,03104 

20*57435 

,07052 

-,11790 

20,57435 

-,07052 

,03293 

21,59151 

,03523 

-,10704 

21,59151 

-,03523 

,04904 

22,05696 

*02743 

-,09984 

22,05696 

-,02743 

,06060 


7UP 

XSPAN= 
ITF.U= 19 
ZTE= 

UPPER SURFACE 
rUF CPU 

6,00000 
ITEL^ 19 
23,08543 

LOWER 

ZLO 

SURFACE 

YLO 

CPL 

9,23919 

-,00000 

-,04739 

9,23919 

-.00000 

-.04739 

9,50436 

,10638 

-,20300 

9,50436 

-.10638 

,24458 

10,38997 

,20062 

-,20370 

10,38997 

-,20062 

,13207 
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